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PREFACE. 


The  object  of  the  following  treatise  is  to  exhibit,  in  a  concise 
form,  the  elementary  properties  of  the  expressions  known  by 
the  name  of  Laplace's  Functions,  or  Spherical  Harmonics, 
and  of  some  other  expressions  of  a  similar  nature.  I  do  not, 
of-  course,  profess  to  have  produced  a  complete  treatise  on 
these  functions,  but  merely  to  have  given  such  an  introduc- 
tory sketch  as  may  facilitate  the  study  of  the  numerous 
works  and  memoirs  in  which  they  are  employed.  As 
Spherical  Harmonics  derive  their  chief  interest  and  utility 
from  their  physical  applications,  I  have  endeavoured  from 
the  outset  to  keep  these  applications  in  view. 

I  must  express  my  acknowledgments  to  the  Rev.  C.  H. 
Prior,  Fellow  of  Pembroke  College,  for  his  kind  revision  of 
the  proof-sheets  as  they  passed  through  the  press. 

N.  M.  FERRERS. 

GONVILLE  AND  CaicS   COLLEGE, 
August,  1877. 
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.  CHAPTER  I.  ,       , 

INTRODUCTORY.      DEFINITION  OF   SPHERICAL  HARMONICS. 

1.  If  V  be  the  potential  of  an  attracting  mass,  at  any 
point  X,  y,  z,  not  forming  a  part  of  the  mass  itself,  it  is 
known  that  V  must  satisfy  the  differential  equation 

fJ^V    (J^V    fPV 

or,  as  we  shall  write  it  for  shortness,  V'F=  0. 

The  general  solution  of  this  equation  cannot  be  obtained 

in  finite  terms.    We  can,  however,  determine  an  expression 

which  we  shall  call   F<,  an  homogeneous  function  oix,y,  z 

of  the  degree  i,  i  being  any  positive  integer,  which  will 

satisfy  the  equation ;  and  we  may  prove  that  to  every  such 

solution  F,  there  corresponds  another,  of  the  degree  —  (i  +  1), 

V 
expressed  by  -^^ ,  where  ?-'  =  a:*  +  y'^  +  z^. 

For  the  equation  (1)  when  transformed  to  polar  co-ordi- 
nates by  writing  x  =  r  sin  Q cos ^^  y  =r  sin  Q  sin  (p,  z  =  r  cos  0, 
becomes 

d'{rV)  ,     1      d  (^._a^V\  ,      1     d^V 


d 


sm  a  da  \         do  J      sm  0  d(p 

And  since  V  satisfies  this  equation,  and  is  an  homo- 
geneous function  of  the  degree  i,  F|  must  satisfy  the  equa- 
tion 

/  ^         '     *     am  6  dd  \  dd  J      sm^  d(p 

F.  H.  1 
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since  this  is  the  form  which  equation  (2)  assumes  when  V 
is  an  homogeneous  function  of  the  degree  «*. 

Now,  put  Fj.  =  r**"^*  Uf,  and  this  becomes 


or 


'  ('  +  ^)  ^  +  sfiT^  31  (^'°  *  W)  +  S?^  5^  =  « (2)- 

Now,  since  Z7^.  is  a  homogeneous  function  of  the  degree 

dr^  dr 

or  r— ^  =  *(*+l)C7.: 

therefore  equation  (2)  becomes 

d'irU)         1      J    /  .       Jf/A  1     d'U,_ 

''     ar     +  sin  e  dd  V""  ^  de)^  sin^  6  df  ~   ' 

shewing  that   CT.  is  an  admissible  value  of  V,  as  satisfying 
equation  (2). 

It  appears  therefore  that  every  form  of  U^  can  be  ob- 
tained from  V^,  by  dividing  by  r"*'^\  and  conversely,  that 
every  form  of  V^  can  be  obtained  from  11^  by  multiplying 
by  r^''*'^.  Such  an  expression  as  F^  we  shall  call  a  Solid 
Spherical  Harmonic  of  the  degree  i.  The  result  obtained 
by  dividing  F^  by  r*,  which  will  be  a  function  of  two  inde- 
pendent variables  ^  and  ^  only,  we  shall  call  a  Surface 
Spherical  Harmonic  of  the  same  degree.  A  very  important 
form  of  spherical  harmonics  is  that  which  is  independent 


DEFINITION  OF  SPHERICAL  HARMONICS.  3 

of  (f).  The  solid  harmonics  of  this  form  will  involve  two  of 
the  variables,  x  and  y,  only  in  the  form  x^  +  y^,  or  will  be 
functions  of  o?  +  y^  and  z.  Harmonics  independent  of  ^  are 
called  Zonal  Harmonics,  and  are  distingnished,  like  spherical 
harmonics  in  general,  into  Solid  and  Surface  Harmonics. 
The  investigation  of  their  properties  will  be  the  subject  of 
the  following  chapter. 

The  name  of  Spherical  Harmonics  was  first  applied  to 
these  functions  by  Sir  W.  Thomson  and  Professor  Tait,  in 
their  Treatise  on  Natural  Philosophy.  The  name  "  Laplace's 
Coefficients"  was  employed  by  Whewell,  on  account  of 
Laplace  having  discussed  their  properties,  and  employed 
them  largely  in  the  Mecanique  Celeste.  Pratt,  in  his 
Treatise  on  the  Figure  of  the  Earth,  limits  the  name  of 
Laplace's  Coefficients  to  Zonal  Harmonics,  and  designates 
all  other  spherical  harmonics  by  the  name  of  Laplace's 
Functions.  The  Zonal  Harmonic  in  the  case  which  we  shall 
consider  in  the  following  chapter,  i.e.,  in  which  the  system 
is  symmetrical  about  the  line  from  which  0  is  measured, 
was  really,  however,  first  introduced  by  Legendre,  although 
the  properties  of  spherical  harmonics  in  general  were  first 
discussed  by  Laplace;  and  Mr  Todhunter,  in  his  Treatise, 
on  this  account  calls  them  by  the  name  of  "Legendre'jj 
Coefficients,"  applying  the  name  of  "  Laplace's  Coefficients " 
to  the  form  which  the  Zonal  Harmonic  assumes  when  in 
place  of  cos  ^,  we  write  cos  ^cos^' +  sin  ^sin^'cos(0  — ^'). 
The  name  "  Kugelfunctionen "  is  employed  by  Heine, 
in  his  standard  treatise  on  these  functions,  to  designate 
Spherical  HarmoDics  in  general. 


1—3 


CHAPTER  11. 

ZOXAL   HARMONICS. 

1.  We  shall  in  this  chapter  regard  a  Zonal  Solid  Har- 
monic, of  the  degree  i,  as  a  homogeneous  function  of 
{x^  +y^)^,  and  z,  of  the  degree  i,  which  satisfies  the  equation 

dW     d'V     d'V_ 
'  dx'  "^  df  "^  dz' 

Now,  if  this  be  transformed  to  polar  co-ordinates,  by 
writing  r  sin  6  cos  ^  for  x,  r  sin  0  sin  (f>  for  i/,  r  cos  0  for  z,  we 
observe,  in  the  first  place,  that  x^  +  y^  —  r^  sin^  0.  Hence 
V  will  be  independent  of  <^,  or  will  be  a  function  of  r 
and  6  only.  The  differential  equation  between  r  and  0 
which  it  must  therefore  satisfy  Avill  be 


r 


-^^+^6de[''''^d0)==^- 


Now  V,  being  a  function  of  r  of  the  degree  i,  may  be 
expressed  in  the  form  r^F^,  where  P^  is  a  function  of  0  only. 
Hence  this  equation  becomes 

or,  putting  cos  B  =  fi, 

In  accordance  with  our  definition  of  spherical  surface 
harmonics,  P^  will  be  the  zonal  surface   harmonic   of  the 
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degree  i.    When  it  is  necessary  to  particularise  the  variable 
involved  in  it,  "we  shall  write  it  P,  (/i). 

The  line  from  which  6  is  "measured,  or  in  other  words 
for  which  /m  =  1,  is  called  the  Axis  of  the  system  of  Zonal 
Harmonics;  and  the  point  in  which  the  positive  direction 
of  the  axis  meets  a  sphere  whose  centre  is  the  origin  of 
co-ordinates,  and  radius  unity,  is  called  the  Pole  of  the 
system. 

Any  constant  multiple  of  a  zonal  harmonic  (solid  or 
surface)  is  itself  a  zonal  harmonic  of  the  same  order. 

2.  The  zonal  harmonic  of  the  degree  i,  of  which  the 
line  //,  =  1  is  the  axis,  is  a  perfectly  determinate  function  of 
fjb,  having  nothing  arbitrary  but  this  constant.  For  the 
expression  r'P,  may  be  expressed  as  a  rational  integral 
homogeneous  function  of  r  and  z,  and  therefore  F^  will  be 
a  rational  integral  function  of  cos  d,  that  is  of  fi,  of  the 
degree  i,  and  will  involve  none  but  positive,  integral  powers 

of  fJi. 

But  Pj  is  a  particular  integral  of  the  equation 

i\i^-'^'^^}  +  ^^'  +  ^^fM  =  0 (3). 

and  the  most  general  form  of  f(fi)  must  involve  two  ar- 
bitrary constants.      ISuppose  then   that  the   most  general 

form  oifiji)  is  represented  by  P^  Ivdfi.    We  then  have 

^  (1  -  ^•)  ^^ = (1  -  ^)  g  .jvdf. + (1  -  /)  p,v. 

d  r,_     ,,  <?./(/*))     d  f„     ,^  dP)  f  , 


+  2(l-^')g.  +  P.||(l-^').| 


Hence,  adding  these   two  equations  together,   and    ob- 
serving that,  since  I\  satisfies  the  equation  (3),  the  coefficient 
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of  Ivclfju  will  be  identically  eq[ual  to  0,  we  obtain,  for  the  de- 
termination of  V,  the  equation 

whence      -P.  (1 -/>'') ^  + 2 1(1 -/='')  ^-/^^.j'-^O.   . 

dv 
or 


the  integral  of  which  is 

log  t?  +  log  Pi  (1  -  /w.")  =  log  (7^  =  a  constant ; 

a 


v  = 


Hence  jvdfi  =  C+C^j  p,  ,i^_-^.  5 

and  we  obtain,  for  the  most  general  form  oif{fi), 

Now,  Pi  being  a  rational  integral  function  of  //.  of  i 
dimensions,  it  may  be  seen  that  i,   _   -2]  p-^  "^^^^^  assume  the 

form  of  the  sum  of  i  +  2  logarithms   and   i  fractions,  and 

therefore  cannot  be  expressed  as  a  rational  integral  function 

c  J 

of  /A.    Expressions  of  the  form  P^  I  j- .    p.^  are  called  Kugel- 

functionen  der  zweiter  Art  by  Heine,  who  has  investigated 
their  properties  at  great  length.  They  have,  as  will  hereafter 
be  seen,  interesting  applications  to  the  attraction  of  a  sphe- 
roid on  an  external  point.  We  shall  discuss  their  properties 
more  fully  hereafter. 

3.     We  have  thus  shewn  that  the  most  general  solution 
of  equation  (2)  of  the  form  of  a  rational  integral  function  of  a 
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involves  but  one  arbitrary  constant,  and  that  as  a  factor. 
We  shall  henceforth  denote  by  Pj,  or  P<  {fju),  that  particular 
form  of  the  integral  which  assumes  the  value  unity  when  fi 
is  put  equal  to  unity. 

We  shall  next  prove  the  following  important  proposition. 

Ifh.be  less  than  uniti/,  and  if  (1  —  2/-th  +  h*)~*  he  expanded 
in  a  series  'proceeding  hy  ascending  powers  of  h,  the  coefficient 
of  h}  will  be  Pp 

Or,  (l-2/ih  +  hT^  =  P,  +  P,h  +  ...  +  P,h'  +  .,. 

We  shall  prove  this  by  shewing  that,  if  H  be  written  for 
(1  —  2/j.h  +  h^)~^,  H  will  satisfy  the  differential  equation 

For,  since  H=  (1  -  2/ti/i  +  h^^, 

.:  ^-,  =  l-2fjLh  +  h'; 

XdH 
'''  ti' df4,~' 

h  dfi  * 

=  -2fjLip  +  s{i-fjL')hrr. 

,    ,  I  dH  J 

^"^  H^dh^f"-^'' 

=  H'[l-2fih  +  h''  +  h(iJ.-h)] 
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.•:|,(A^)=^i^'(i-M)}   '■::'.:■; 

=  3  (1  -  till)  II'  ifi  -  h)  -fiTP. 

=-  -  S/xH'  +  3  {(1  -  lii')  h+(l-/JLh)(jj,-  h)]  IP 
=  -  ti^lIP  +  3  {/a  (1  +  h')  -  2fj.'h]  H' 
=  -Sfi{IP-il-2fxh  +  h')H'} 
=  0,  since  1  -  2fxh  +  Ji'  =  H^. 

Therefore  J  jd  -  Z'')  f } +/' |.  (^ff)  =  0- 

This  may  also  be  shewn  as  follows. 

If  X,  y,  z  be  the  co-ordinates  of  any  point,  /  the  distance 
of  a  fixed  point,  situated  on  the  axis  of  ^s-,  from  the  origin, 
and  R  be  the  distance  between  these  points,  we  know  that, 


and  that 


©-• 


Now,  transform  these  expressions  to  polar  co-ordinates, 
by  Avriting 

x  =  r  sin  &  cos  ^,     y  —  r  sin  6  sin  ^,     z  =  r  cos  9, 

and  we  get 

R^  =  r^-2zrcosd  +  z'\ 

and  the  differential  equation  becomes 

d^  (r\        1      d  {  .    .d  (1\\      ^ 

'-dr\R)-^^eder^d6[li)r^^ 


ZONAL  HARMONICS, 
or,  putting  cos  6  =  im, 

'■|'(a)+i{(^-''''i(i 

Now,  putting  r  =  zh,  we  see  that  , 

R^  1 

^,  =  ^'-2^  +  1  =  7^., 

1     H 

or  P  =  "7  > 


=  0. 


.*.  the  above  equation  becomes 
z  dk 


or 


4.     Having  established. this  proposition,  we  may  proceed 
as  follows: 

If  p,  be  the  coefficient  of  A'  in  the  expansion  of  //, 
//=  1 +pj/i +p,/t' +  ... +_?Vt' +  ... 
/.  hH=  h  4-  iV^'  +i3,/t''  +  . . .  +  p,h"'  + . . . 

Also,  the  coefficient  of  7i*  in  the  expansion  of 
d   {  i.dH)    .      d   (,^        ^^dp) 

Hence  equating  to  zero  the  coefficient  of  A', 
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Also  p^  is  a  rational  integral  function  of  /i. 
And,  when  /x  =  1,  //=  (1  -  2A  +  k^^ 

=  l+h  +  h^  +  ...+U  +  ... 

Or  when  /i  =  1,  p,  =  1. 
Therefore  p,  is  what  we  have  already  denoted  by  T^^. 

We  have  thus  shewn  that,  if  h  be  less  than  1, 

If  h  be  greater  than  1,  this  series  becomes  divergent. 
But  we  may  write 

since  ,  is  less  than  1, 

P      P  P. 

Hence  P^  is  also  the  coefficient  of  ^-(*+i>  in  the  expan- 
sion of  (1  —  2/A/t  +  /i")   '^  in  ascending  powers  of  y-  when  h  is 

greater  than  1.     We  may  express  this  in  a  notation  which  is 
strictly  continuous,  by  saying  that 

Pi  —  P-[i+\]' 

This  might  have  been  anticipated,  from  the  fact  that  the 
fundamental  ditferential  equation  for  P;  is  unaltered  if 
—  (i+l)  be  written  in  place  of  ^;  for  the  only  way  in 
which  i  appears  in  that  equation  is  in  the  coefficient  of  P^ 
which  is  i  (i-\- 1).  Writing  —  {i  +  1)  in  place  of  i,  this  be- 
comes —  {i  +  Vj  [—  [i  +  1).  +  1]  or  (i  +  1)  i,  and  is  therefore 
unaltered. 
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5.     "VVe  shall  next  prove  that 

^*     ^     ^  1.2...idz'[rJ' 
where  7'^  =  x^  +  y*  +  z^. 

Let  ^=(a:'  +  2/'  +  ^r'=/(^), 

and  let  k  be  any  quantity  less  than  r.       . 

Then  {a?-\-f  +  {z-Tcf]-^=f[z-k), 

and,  developing  by  Taylor  s  Theorem,  the  coefficient  of  k'  is 

Also      {a;'  +  3/"  +  (^  -  ^)T-  =  (^'  -  2^2  +  k')'^ 

1/,    ^  ^•   ^•v^ 


1/,     ^   k     k'Y^ 


since  z  =  fir, 

in  the  expansion  of  which,  the  coefficient  of  /c'  is 

Equating  these  results,  we  get 

^*        ^  ^    1.2...td3'W 

The  value  of  Pj  might  be  calculated,  either  by  expanding 
(1  —  2fili  +  h^y^  by  the  Binomial  Theorem,  or  by  effecting  the 

differentiations  in  the  expression  (— 1)' ,— o   « .-t-J-    1 

1.2.3...  taz\rj 

and  in  the  result  putting  -  =  /i.  Both  these  methods  how- 
ever would  be  somewhat  laborious  ;  we  proceed  therefore  to 
investigate  more  convenient  expressions. 
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6.     The  first   process   shews,  by  the  aid  of  Laigi-ange's 
Theorem,  that 

1         ,  ^*      , 

■^' "  2M.2.3...^  ^-  ^^'  ~  ^^'* 
Let  y  denote  a  quantity,  such  that 

-I'V-     - 

h  being  less  than  1. 

Then 

1 
dy  _  h 

dpi 


3^  =  x-li-X+AV' 


Also  [y-^^l^^^^,:,  ■ 

'-y   h      k ' 

.         /.  3/  =  ^  +  a(^-1).         .  . 

Hence,  by  Lagrange's  Theorem, 

therefore,   differentiating  with  respect  to  /a  and  observing 
that 

i  ,*2  ...  id/jb  \    I    J 
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7.  From  this  form  of  P^  it  may  be  readily  shewn  that 
the  values  of  fx,,  which  satisfy  the  equation  P^  =  0,  are  all  real, 
and  all  lie  between  —  1  and  1. 

For  the  equation  "  ' 

(/t*  —  1)*  =  0  has  %  roots  =  1,  and  i  roots  =  —  1, 

.'.  -7-  ifJi'^  —  1)*  =  0  has  i  —  1  roots  =  1,  {i  —  1)  roots  =  —  1,  and 
one  root  =  0,  ■ 

-T-2  (/«.'  — 1)*  =  0  has  (t  — 2)  roots  =  1,   one   root  between   1 

and  0,  one   between  0  and  =  —  1,  and   {i  —  2)    roots  =  —1, 
and  so  on.     Hence  it  follows  that 

d*  i  i 

T~i  i}^  —  1)*  =  0  has  Q  roots  between  1  and  0,  and  -  roots  be- 

tween  0  and  —  I,  if  i  be  even, 

i—1  i—l 

and  roots  between  1  and  0,  — ^—  roots  between  0  and 

1,  and  one  root  =  0,  if  i  be  odd. 

It  is  hardly  necessary  to  observe  that  the  positive  roots  of 
each  of  these  equations  are  severally  equal  in  absolute  mag- 
nitude to  the  negative  roots. 

8.  We  may  take  this  opportunity  of  introducing  an  im- 
portant theorem,  due  to  Rodrigues,  properly  belonging  to 
the  Differential  Calculus,  but  which  is  of  great  use  in  this 
subject. 

The  theorem  in  question  is  as  follows: 
If  rale  any  integer  less  than  i, 

dx'--  ^*      -^^  -  1 .  2  ...  (i  +  m)  ^"^      ^^    duT-  ^"^      ^^- 
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It  may  be  proved  in  the  following  manner. 

If  (a;-—!)*  be  differentiated  i-m  times,  then,  since  the 
equation 

has  «  roots  each  equal  to  1,  and  t  roots  each  equal  =  —  1,  it 
follows  that  the  equation 

Ji-m 


dx 

has  i  —  (i  —  m)  roots  (i.  e.  m)  roots  each  =  1,  and  m  roots 
each  =  —  1,  in  other  words  that  (x^  —  1  j"'  is  a  factor  of 

We  proceed  to  calculate  the  other  factor. 

For  this  purpose  consider  the  expression 

{x  +  aj  (x+a;)  ...  (a;  +  a,)  (a;  +  y8J  (x  +  yS,) ...  (^  +  /3,). 

Conceive  this  differentiated  (I)  i  —  m  times,  (II)  i  +  m 
times.  The  two  expressions  thus  obtained  will  consist  of  an 
equal  number  of  terms,  and  to  any  term  in  (I)  will  corre- 
spond one  term  in  (II),  such  that  their  product  will  be 
(aj  +  otj)  (^4-aJ...(a;  +  «.)  (a;  +  /3J  («  +  /3,)  ...  (a;  + /3,),  i.e.  the 
term  in  (II)  is  the  product  of  all  the  factors  omitted  from 
the  corresponding  term  in  (I)  and  of  those  factors  only. 
Two  such  terms  may  be  said  to  be  complementary  to  each 
other. 

Now,  conceive  a  term  in  (II)  the  product  of  p  factors  of 
the  form  x  +  a,  say  x+a.',  x  +  a'  ...x  +  a'^',  and  of  q  factors 
of  the  form  x+/3,  say  x  +  ^^,  x  +  /3^^...x  +  /S^^y  We  must 
have  2)  +  q  =  i  —  m. 

The  complementary  term  in  (I)  will  involve 

p  factors  x  +  ^,  x  +  jS"  ...  a  +  yS''", 

q  factors  x  +  a^,  x  +  a^, ...  x  +  Ui^^ 
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Now,  every  term  in  (I)  is  of  i  +  m  dimensions.  We  have 
accounted  for  p+  q  (or  i  —  m)  factors  in  the  particular  term 
we  are  considering.  There  remain  therefore  2m  factors  to 
be  accounted  for.     None  of  the  letters 

can  appear  there.     Hence  the  remaining  factor  must  involve 
m  a's  and  m  ^S's, — say, 

There  will  be  another  term  in  (II)  containing 
{cc  +  ^')  {x  +  /3")  ...{x  +  ^<-')  {x  +  a)  (x  +  a  J  ...(x  +  a  J. 
The  corresponding  terra  in  (I)  will  be,  as  shewn  above, 
(x  +  a')  {x  +  a") ...  {x+o!")  (x  +  ^)  {x  +  ^^)  ...  (x  +  ^J 

ix  +  ,a)  (x  +  .-x)  ...  {x  +  ^i)  {x  +  A  i^  +  A) '"  {^  +  ..^)- 

Hence,  the  sum  of  these  two  terms  of  (I)  divided  by  the 
sum  of  the  complementary  two  terms  of  (II)  is 

{x  +  ,0)  (x  +  ^a)  ...  {x  +  ^oi)  {x  +  ^0)  (x  +  A)  -  i^+.^)- 

Now,  let  each  of  the  a's  be  equal  to  1,  and  each  of  the  ^S's 
equal  to  —  1,  then  this  becomes  (a^  —  1)'".  The  same  factor 
eaters  into  every  such  pair  of  the  terms  of  (I).     Hence 

dx'-^        ^  ^     "  ^   da^^"'       '        ^  numerical 

factor  prh. 

The  factor  may  easily  be  calculated,  by  considering  that 

the  coefficient  of  a;'^  in       J  .^„  ^  is2{(2i-l)...(t+7?i+l), 

and  that  the  coefficient  of  ic*"^  in -T~+r. — -  is 

dx' 

2i(2*-l) ...  (t  +  7n  +  l)  (i  +  w)...  ({-w  +  1). 
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Hence  the  factor  is  .      '' 

1  1.2... (t-m) 

or ^^ - 

(t  +  w)  (i  +  m-1)  ...  (i-m  +  1)'       1.2...  (t  +  m)*..    .U 

9.    This  theorem  affords  a  direct  proof  that  C-7— ^  {j^ —  \y, 

C  being  any  constant,  is  a  value  of  /  (/a)  which  satisfies  the 
equation 


rf/i 


d     d^ 


d' 


(W 


// 


\  ^'\     ^-^dy.-dfi'^-^=^-'^dy^^-^ 


from  above, 

••■i[f'^'-i)il|.('^'-i)' 


=;(i+i)|;.(/-i)' 


=i(i+i){|.(/''-i)'}. 


or 

d_ 
dfjL 


Hence,  the  given  differential  equation  is  satisfied  by  put- 
tmg/Oa)  =  a|^.(/x^-iy. 

Introducing  the  condition  that  P^  is  that  value  of  /  (//) 
which  is  equal  to  1,  when  /a  =  1,  we  get 


P.= 


1  d' 


2M.2.../  d/jt,' 


10.  We  shall  now  establish  two  very  important  proper- 
ties of  the  function  P. ;  and  apply  them  to  obtain  the  develop- 
ment of  Pf  in  a  series* 
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The  properties  in  question  are  as  follows  :• 
If  i  and  m  ha  unequal  positive  integers, 

And  f_p,',,  =  ^^. 

The  following  is  a  proof  of  the  first  property^ 
We  have 

I  {(!-''')  fh '■('+!)  ^'=«' 

Multiplying  the  first  of  these  equations  by  P„,  the  second 
by  P„  subtracting  and  integrating,  we  get 

+  [i  {i  +  1)  -  m  (m  +  1)|  jP.PJfi  =  0. 

Hence,  transforming  the  first  two  integrals  by  integration 
by  parts,  and  remarking  that 

i  {i  +  1)  —m  {m  +  1)  =  (i  —  m)  {i+  m  + 1), 
we  get 

(-')^f-^.t)-/(-''')(tf-fS>'' 

+  (i  -  m)  {{  +  m  +  l)jP,PJfji  =  0, 


or 


(1  -  /.')  (P„  ^^'  -  P,  ^-)  +  {{  -m){i  +  m  + 1)  jP,PJf.=0, 

since  the  second  term  vanishes  identically. 

r.  H.  2 
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Hence,  taking  the  integral  between  the  limits  —  1  and 
+  1,  we  remark  that  the  factor  1  —  //.^  vanishes  at  both  limits, 
and  therefore,  except  when  i  —  m,  or  i  +  m+l  =  i), 

-1 

We  may  remark  also  that  wo  have  in  general 


/: 


/, 


pdP,_p  dP„ 


^  {i—m)(i+m  +  l)' 

a  result  which  will  be  useful  hereafter. 

11.     We  will  now  consider  the  cases  in  which 

i  —  m,  or  i  +  m+  1  —  0. 

We  see  that  i-f  m  + 1  cannot  be  equal  to  0,  if  *  and  m  are 
both  positive  integers.  Hence  we  need  only  discuss  the 
case  in  which  7n  =  i.     We  may  remark,  however,  that  since 

P^  =  P_(i+i),  the  determination  of  the  value  of  I    P]  dfx,  will  also 

n 
give  the  value  of  /     P^  P_^^.^)  dfi. 

The  value  of  1     P/J/t  may  be  calculated  as  follows; 
(l-2/x7t+;0-^  =  Po  +  P,^  +  ...  +  P/i*  +  ...; 

=  Po^  +  p^'h'  +  ...  +  p;h''+... 

+  2P,PJi  +  2P,PJi'  +...+2P^PJi'  +  ... 
Integrate  both  sides  with  respect  to  //. ;  then  since 


/ 


(1  -  2/j,h  +  hY  diJ^  =  -li  log  (1  -  2/i^^  +  ^0> 


2^ 
we  get,  taking  this  integral  between  the  limits  —  1  and  +  1, 

all  the  other  terms  vanishing,  by  the  theorem  just  proved. 
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Hence,  equating  coefficients  of  A", 


/ 


>''^''  =  2.-  +  l- 


12.     From  tlie  equation  I     PjP„,tZyLt  =  0,  combined  with 

the  fact  that,  when  fi  =  l,  P^  =  1,  and  that  Pj  is  a  rational 
integral  function  of  /j,,  of  the  degree  i,  P,  may  be  expressed 
in  a  series  by  the  following  method. 

We  may  observe   in  the  first  place  that,  if  m  be  any 

integer  less  than  i,  I     jj^Tidfi  =  0. 

For  as  P„,  P^_^ . . .  may  all  be  expressed  as  rational  in- 
tegral functions  of  fi,  of  the  degrees  vi,  m  —  1  ...  respectively, 
it  follows  that  //.'"  will  be  a  linear  function  of  P^  and  zonal 
harmonics  of  lower  orders,  /t'""*  of  P^^  and  zonal  harmonics  of 

lower  orders,  and  so  on.     Hence  jfi^P^dfji,  will  be  the  sum  of 

a  series  of   multiples  of  quantities  of  the  form   1  PJP.djx, 

m  being  less  than  i,  and  therefore  I     fJ^P^dfx,  =  0,  if  ??i  be  any 
integer  less  than  i. 
Again,  since 

(l-2/.^  +  /0"^  =  Po  +  P./H-...+P,/V+... 
it  follows,  writing  —  h  for  h,  that 

{l  +  2/Mh  +  h')-'-  =  P,-PJi+...  +  i-iyPJU... 

2—2 
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And  writing  —  /^  for  //,  in  the  first  equation, 

P„',  P/...P/...    denoting  the    values  which  P„,  P,...Pi, 

respectively  assume,  when  —/a  is  written  for  jx.  Hence 
P^=P^  or  —Pi,  according  as  i  is  even  or  odd.  That  is, 
P.  involves  only  odd,  or  only  even,  powers  of  i,  according 
as  f  is  odd  or  even*. 

Assume  then 

P,  =  A,t.^  +  A,_,fj}^+... 

Our  object  is  to  determine  A^,  Ai_^.... 

Then,  multiplying  successively  by  ya'~^,  /u,*"*,  ...  and  inte- 
grating from  —  1  to  + 1,  we  obtain  the  following  system  of 
equations : 

'■     ..   •^'-2    1.       I :^»-2» 1       —  f> 


2i  -  1     2{  -  3  2i  -  2s  -  1 


J !=i-    4-      .     J '-2' L  =0 


2i-S     2i-5  2i-2s-3 


^-     +^^^-^+...+^r"— r  +  ..•  =  0. 


2z-2s-l  '  2i-2s-3  2t-4s-l 


And  lastly,  since  P^  =  1,  when  fi  =  l, 

^i  +  A_2  +  -.-  +  ^i-2.+ •.•  =  !; 
the  last  terms  of  the  first  members  of  these  several  equa- 
tions being 

AAA 

-^\,   ^-\...-:.^,A^,iUheexen, 
t— 1      i—o       1 

^^1-       ^1        _J    ^     if?;  be  odd 

13,  The  mode  of  solving  the  class  of  systems  of  equa- 
tions to  which  this  system  belongs  will  be  best  seen  by 
considering  a  particular  example. 

*  Tliis  is  also  evident,  from  the  fact  that  P,  is  a  constant  multiple  of 
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Suppose  then  that  we  have 


X  If  z  ^ 


a+a      6+a     c+a 
X  y  z        ^ 


a  +  /3     i  +  /3     c  +  y3 

jg  +  y  ^  g  ^  1 

a  +  <u      i+w      c  +  a)      o)  * 

From  this  system  of  equations  we  deduce  the  following, 
6  being  any  quantity  whatever, 

X         j^     _s_  ^  \  {0-a)  (0-/3)  (a +  0))  {b  +  Qj)  (c+(o) 
a+0     b  +  d'^c  +  9~  0)  {(o-a){(o-l3){a+e)  {b+d)  (c+O) ' 

For  this  expression  is  of  —  1  dimension  in  a,  b,  c,  a,  /S,  y, 
0,  w;  it  vanishes  when  0  =  a,  or  0  =  ^,  and  for  no  other 

finite  value  of  0,  and  it  becomes  =  — ,  when  0  =  «. 

a 

We  hence  obtain 

'^'^^''^^\b+0'^ 0+0)      o>  («-a)(<y-)8)       (b  +  0){c+0)       ' 

and  therefore,  putting  0  =  —a, 

_  1  (g  +  a)  (g  +  /3)  (a  +  <»)(&+  &>)  (c  +  «) 
~  ci)  (a  —  i)  (a  —  c)         (w  —  a)  (&)  —  /3)       ' 

with  similar  values  for  y  and  s. 

And,   if  ft)  be  infinitely  great,  in  which  case  the   last 
equation  assumes  the  form  a;  +  y  +  ^^  =  1,  we  have 

^~\a-b){a-c)' 
with  similar  values  for  y  and  z. 

14.     Now  consider  the  general  system 
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rjn  /v. 


^'     4-     ^'-8— 4....  4.     ^*-^»      4....=!  . 

the  number  of  equations,  and   therefore  of  letters  of  the 

■i  + 1  .  i 

forms  X  and  a,  being  if  i  be  odd,  ^  +  1  if  i  be  even ;  and 

i—\ 

the  number  of  letters  of  the  form  a  being  — - —  '\ii  be  odd, 

and  ^  —  1  if  I  be  even. 

A 
We  obtain,  as  before, 

. » 1 ir9 L  4.  •-?^ L 

^  1^  (^,)  {e-a,^ . . .  (^-a,. J^  (a.^o,)  (a,.,+  e.)^ (a,_,,+a,) . . . 
0)  (a,-a,)(«-a,_ J . . . (a,-a,_ J . . .  (a,+  ^)  (a,_,  +  ^) . . .  («,_,.+  i?) . . . 

and,  multiplying  by  a^_^,  +  ^,  and  then  putting  6  =  —  a,._2„ 

a-  =  1  («f-2..  +  «»•)  (o^i-2.  +  «i- J  '  •  •  («i-2»  +  «.--2»)  •  •  • 

*-2'     o)        (to  -  a .)  (ft)  -  a,.  J  . . .  (g>  -  a._  J . . . 


{a,  +  0))  (a,_j  +&))...  (a,_2,  +  «) . 


(«i-2.  -  «i)  (^.-2.  -  «i-2)  •  •  •  («.-2.  -  «i  or  «o)  ' 

15.  To  apply  this  to  the  case  of  zonal  harmonics,  we  see, 
by  comparing  the  equations  for  x  with  the  equations  for  Ay 
that  we  must  suppose  o)  =  00  ;  and 

Oi  =  i  ai_2  =  »  -  2,...a^2.  =  «  -  25  ... 

Hence 

.      _  (2^•  -  2.9  - 1)  (2t  -  2g  -  3). ..{2  {i  -  2g)  - 1|... 
•-2'  -  (_  25)  {-  (25  -  2)}...{(i  -  2s  -  1)  or  {i  -  2s)  j 

(2i-2s-l)  (2i-2s-3)...{2(t-2s)-l|... 


=  (-!)« 


2s(2s-2)...2x2.4...(i-2s-l)  or  (t-2s)* 
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Or,  generally,  if  i  be  odd, 

_  (2^•-l)(2^-3)...(^  +  2) 
'~  2.4...(i-l) 

A      -  _  (2^'-3)(2^-5)...^ 
*-=^~       2.4...(i-3)x2"' 

^(2*-^5)_{2t-7),^,(t^) 
^-*         2.4...(t-5)x2.4     ' 


'~^     ^     2.4...(i-l)* 
And,  if  i  be  even, 

_^2^-lK2t-^)^^/^+l) 
'"  2.4...i 

_     (2^-3)(2^-5)...(^-l) 


A-.= 


2.4...(t-2)x2 
(2i-5)(2i-7)...(t-3) 
2.4...(*-4)  X  2.4      ' 

We  give  the  values  of  the  several  zonal  harmonics,  from 
P,  to  P,o  inclusive,  calculated  by  this  formula, 

3    ,     1 


i". 

"       2 

=  2^- 

"2' 

3/i'- 

2 

■1 
> 

p. 

5       3 

=  2^" 

3 

"2^ 

5/i,'- 

-3^ 
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p_7.5    ,     5.3   „     3.1 

8 
p_9.7   5     7.5    3_^5.3 

~  8  ' 

'       2.4.6^      2.4  x:i^  "^"2x2.4^      2.4.0 

_  231/x^  -  315;/,^  +  105yL6^  -  5 
16 

P  _  13. 11 .9   ,     11.9.7   5     9.7.5    ^     7.5.3 
^~    2.4.(5   '^~2.4x2^''"2x2.4^~274Tg'^ 
_  429yu,^  -  G9Sf/  +  315/^^  -  35/^ 
16 

_  15. 13^  11.^   8     13.11.9.7   g     11.9.7.5    , 
■^'  ~     2.4.678    ^       2.4.6x2^  ■*'2. 4x2. 4^ 

9.7.5.3    ,  \  7.5.3.1 
2  X  2 . 4 . 6  ^  "*■  2 . 4 . 6 . 8 

_  6435/  -  12012/x'^  +  QdSOfi*  -  1260yu,''  +  35 

128  '       . 

17.15.13.11  3  15.13.11.9  7  ,13.11.9^7  , 
'~     2.4.6.8  ^       2.4.6x2  '^  "*■  2.4x2.4  ^ 

11.9.7.5  3  9.7.5.3 
2x2.4.6^  "^2.4.6.8^ 

_  12155ya"-  25740/^^  +  18018a<.^-  4620//,"+  315/^ 

128 

^19^7.15.13.11  ,0  17.15.13.11.9  g  15.13.11.9.7  , 
"~  2.4.6.8.10  ^-  2.4.G.8"x2^'^2.4.6x2.4^ 
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13.11.9.7.5   ,     11.9.7.5.3   ^      9.7.5.3.1 


/*  +Os,o   A  a   ft^ 


2.4x2.4.6^      2x2.4.6.8^      2.4.6.8.10 

,  _  46189/a"-  109395/+  90090/^^-  30030yLt*+  3465/^^-63 
~  256 

It  will  be  observed  that,  wben  these  fractions  are  reduced 
to  their  lowest  terms,  the  denominators  are  in  all  cases 
powers  of  2,  the  other  factors  being  cancelled  by  correspond- 
ing factors  in  the  numerator.  The  power  of  2,  in  the 
denominator  of  P„  is  that  which  enters  as  a  factor  into  the 
continued  product  1 .  2...i. 

16.  We  have  seen  that  /  /jT F..dfjL  =  0,  if  m  be  any 
integer  less  than  i. 

It  will  easily  be  seen  that  ii  m  +  i  be  an  odd  number,  the 
values  of  If/J^F^.dfi  are  the  same,  whether  fju  be  put  =  l  or 

—  1 ;  but  if  m  + 1  be  an  even  number,  the  values  of  I/a'"  P,  .  dfi 

corresponding  to  these  limits  are  equal  and  opposite.   Hence, 
(m  +  i  being  even) 

\  fjrP^.d/jL  =  0,  if  m  =  i-2,  i-4> 

Jo 


and  then 


We  may  proceed  to  investigate  the  value  of  /    fj'"' Pi- dfj, 

Jo 
if  m  have  any  other  value.     For  this  purpose,  resuming  the 
notation  of  the  equations  of  Art.  13,  we  see  that,  putting 
6  =  171  +  1,  and  o)  =  oo  ,  we  have 

+  ^— --''— -T  + + ^ 7  + 


tti+m  +  l     a,_2  +  TO  +  1      '  *"       aj_2,  +  m  +  l 

^  (m  +  1  -  g^)(w  +  1  -  or,J  ...  (m  +  1  -  Q  ...  . 
\a^  +  TO  + 1)  (a^2+  TO  +  1)  ...  {ai_^+m-\-l)  ... ' 
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and    therefore,    putting    x^  =  A^...,    a,  =  » . . .,    a^  =  i  —  1 . . . , 

we  get 

1  A  A  A 

o'^      *     '^     *+m+l"^i-2+m+l'^***"^i-2s+m+l'^'" 

(m -t  +  2)(m-i  +  4)  ...(w-1)         .„  ..        ,. 
—  if  *  be  odd. 


/, 


I 


{m->ri+V)(m  +  i-l)  ...  (7^  +  4)  (w  +  2) 

-                    (m  — «  +  2)  (wi  — 4+4) ..,  m  ._., 

and    =  ;^ .     ■,,  , .     -, ,       ,    — jiT-. z^if  I  be  even. 

(m  4-  *  +  1)  (w  +  «  -  1)  ...  (m  +  3)  (m  +  1) 

In  the  particular  case  in  which  m  =  «*,  we  get 

iD^  2.4  ...  (t-1)  ..     ,,v 

'^    '  '^      (2i+ 1)  (2i— 1)  ...  (i  +  4)  (i  +  2)  ^         ' 

A                                         2.4...Z  ,.         . 

a-nd  =  T^r-; — ...  ,^. — r- -T-. — — — ; — -7-  n even). 

(2i+l)(2i-l)  ...  (t  +  3)(i  +  l)  ^  -^ 

17.  "We  may  apply  these  formtilge  to  develope  any  positive 
integral  power  of  /z-  in  a  series  of  zonal  harmonics,  as  we 
proceed  to  shew. 

Suppose  that  m  is  a  positive  integer,  and  that  /i."*  is  de- 
veloped in  such  a  series,  the  coefficient  of  F^  being  0^,  so 
that 

then,  multiplying  both  sides  of  this  equation  by  P,.  and  inte- 
grating between  the  limits  —  1  and  1,  all  the  terms  on  the 

right-hand  side  will  disappear  except       Cj  P^  d/x,  which  will 

2 
become  equal  to  ^. — ^  C,. 

Hence  C,  =  ?^  j\''F,dfi, 

which  is  equal  to  0,  if  m  + 1  be  odd.     Hence  no  terms  appear 
unless  m  4-  «'  be  even.     In  this  case  we  have 


•I  n 
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Hence  the  formula  just  investigated  gives 
r  -  W4-  ..         (m-t  +  2)(m-i+4)...(m-l) 
*~^     "^    ^(m  +  i+l)(m  +  «-l)...(m  +  4)(m  +  2) 
if  %  be  odd,  and 

C<  _  i^,'  ,  ^^  (m-t  +  2)(7n-^  +  4)  ...  m 

'  -  ^^*  f  ^;  (w+  4+  1)  (^4-  j_i) ...  (m  +  3)  (m+  1) 

if  i  be  even. 

Therefore  if  m  be  odd, 

^_  2.4.6...Cw-l) 

/t   -  (^»»+  1)  (-2;^+  1^  ^2m  -  1)  ...  (m  +  4j  (m+  2)  ^"^  -* 

(m  +  4)(7/i  +  2j^^m  +  2    ^* 

If  m  be  even, 

m_/9        i\  2.4.G...m  p 

/i    -(-^  +  l;(2m+l)  (2w-l)  ...  (to +  3)  (m  +  1)     "•■^•" 

15  ^  P   I      ^      P 

^       (7/1  +  3)  (to  +  1)      '^  ^  7H  +  1      »• 

Hence,  putting  for  to  successively  0,  1,  2  ...  10,  we  get 
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13.11.9.7^11.9.7^9.7'^"^7 

16   „      24  „      10  „      1 
—  P-l —  P  +      P+- 

231     «^77    '     21^^7 


^_.,.      2.4.6       p      .,■,      4.6     p      ^    6  Sp 

'^  ~ ^^  15. 13. 11.  9 ^^"^-^^3. 11. 9^^"^^  11:9 ^« "^9^' 

=  iip+lp+lip+lp 
429^     39    ^33^     3^' 

2.4.6.8  4.6.8 

'^  ~  ^^  17.15.13.11.9^^      15. 13.11. 9^« 

4.  9  _Al^  p    .  5  _^  P  +  ^  P 
^13.11.9^^    11.9^  -9     °' 

_  128^  64^  J.8  40  1 

~  6435    '  '^  495     '  "^  143    '  "^  99     =^  "^  9     «' 

.0=  19  ^-^-g-S p  +15       4.6.8 

^  19.17.15.13.11^^      17.15.13.11^ 

4.11 ^?__p   I  7      ^     p  +s  — P 

^15.13.11^^    13.11^^    11^' 

^128  192  16  56  3 

12155  ^2431  ^  "^  65  ^  "^  143  ^  "^  11  ^' 

.0  _  21 2^4 .  6. 8 J0__  4.6.8. 10_ 

'^   ~       21.19.17.15.13.11     ^o"*"      19.17.15.13.11^ 

,10     6.8.10  _8_^10_p  4-5_l^P  +  J_P 

^17.15.13.11^'^    15.13.11^'*'    13. 11^"^  11     " 

_  _156  m  82  4^      +i^p  +lp 

"46189    ^"^2717    '187    '143    *^  143    =^^11    ° 
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18.  Any  zonal  harmonic  P^  may  be  expressed  in  a  finite 
series  of  cosines  of  multiples  of  $,  these  multiples  being 
ie,  {1-2)6....     Thus 

{l-2fih  +  h;r^  =  P,+PJi+...  +  PJi'+...; 

therefore,  writing  cos  6  for  /jl,  and  observing  that 

l-2coseh  +  h'=(l-  he'^~^ «)  (1  -  ^e-^^«), 

we  obtain 

(1  -  Ae^^^)-*  (1  -  he-  ^^«)-^  =  P^  +  PJi+,„  +  pje+,,^- 

or 

\       2  2.4 

2  .4  ...  2i  / 

=  P,  +  PJi+...  +  PJi,'+.., 

whence,  equating  coefficients  of  h\ 

„      1.3...(2i'-l)„        ..     1.3...(2i-3)l„        ,.     „,  ^ 
*         2.4...2i  2 .4 ...  (2i  — 2j  2  ^        ' 

.  1.3...(2*-5)1.3^        ,.      ,,  . 
+  „- , )^. — -.(  „    .  2  cos  (i  -  4)  (9+  .... 

2 . 4  . . .  (2i  —  4)  2 . 4  ^         ^ 

the  last  term  being  \  —^—z-^^ — -. —  [  if  i  be  even,  and 
°   {     z  .4)  ...I     )  ' 

1.3...(i+])  l^^(-2)_         .....        ,^ 
-      2  .^4T:7(fqrr)  2.4...(^•-l)  2  '"'  ^'  '^  '  ^^  "^^- 

19.     Let  us  next  proceed  to  investigate  the  value  of 

I     P^  cos  md  sin  ^  dd. 
Jo 
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This  might  be  done,  by  direct  integration,  from  the  above 
expression.     Or  we  may  proceed  as  follows. 

The  above  value  of  P^  when  multiplied  by  cos  m6  sin  0 

(that  is  by  ^  (sin  (m+ 1)  ^—  sin  {m  —  1)  6])  will  consist  of  a 

series  of  sines  of  angles  of  the  form  [i  —  2n  +  (m  +  1)}  ^,  that 
is  of  even  or  odd  multiples  of  6,  as  i  +  m  is  odd  or  even. 
Therefore,  when  integrated  between  the  limits  0  and  tt  it 
will  vanish,  if  z'  +  m  be  odd.  We  may  therefore  limit  our- 
selves to  the  case  in  which  t  +  w  is  even. 

Again,  since  cos  mO  can  be  expressed  in  a  series  of  powers 
of  cos  6,  and  the  highest  power  involved  in  such  an  expression 
is  cos"'^,  it  follows  that  the  highest  zonal  harmonic  in  the 

development  of  cos  mO  will  be  P„.    Hence  /    P^  cos  mO  sin  6  dd 

Jo 
will  be  =  0,  if  m  be  less  than  i. 

Now,  writing 

P,  =  Ci cos ie  +  Ci_^ cos  (i-2)e  +  ... 

we  see  that  i^.  cos  md  sin  6  dO  will  consist  of  a  series  of  sines 
of  angles  of  the  forms  {m  +  i+l)6,  (m  +  i  —  l)  0  ...  down  to 
(t/j  _  I*  _  1)  0^  there  being  no  term  involving  m0,  since  the 
coefficient  of  such  a  term  must  be  zero.     Hence 


Jo 


Pi  cos  m0  sin  0  d0, 

0 

will  consist  of  a  series  of  fractions  whose  denominators  in- 
volve the  factors  m  +  i'+l,  m  +  i  —  l...m  —  i  —  1  respectively. 
Therefore  when  reduced  to  a  common  denominator,  the  result 
will  involve  in  its  denominator  the  factor 

(»z  +  /+l)  (m+i-1)  ...  (wi  +  1)  (m-1)  ...  {m-i-l) 
if  m  be  even,  and 

(m  +  «  +  l)(m  +  t-l) ...  (m  +  2)  (771 -2)  ...  {m-i~l) 
if  m  be  odd. 

For  the  numerator  we  may  observe  that  since 

r 

Pj  cos  m  0  sin  0dd 


Jo 
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vanishes  if  m  be  less  than  ?',  it  must  involve  the  factors 
»n  — (i  — 2),  TO— (i  — 4)  ...  7M  +  (i  — 2),  and  that  it  does  not 
change  sign  with  m.    Hence  it  will  involve  the  factor 

[>,n -  {i- 2)}  (m -  {i -  4)}  ...  (w  -  2)  m''  (w  +  2)  ...  (m  +  i  -  2) 

if  TO  be  even,  and 

{w-(i-2)}[to-(;-4)}  ...  (w-1)(to  +  1)  ...  (to +  1-2) 

if  TO  be  odd. 

To  determine  the  factor  independent  of  ??i,  we  may  pro- 
ceed as  follows : 


P,  =  (7,  cos  ie  +  (7,_j  cos  (t  -  2)  (9  +  . . . ; 
1 

2 


.*.   P^  cos  TO^  —  a  ^i  [^OS  (w  +  t)  ^  +  COS  (m  —  l)  $} 


+  2  ^'-a  1^°^  ('^  +  t-  2)  ^  +  COS  (to  - 1  +  2)  ^}  +  ... ; 

.*.  P^  COS  TO^  sin  ^  =  J  C^<  (sin  (m  +  i  +  1)  6  —  sin  (to  +  *  —  1)  ^ 

+  sin  [m  —  i  +  V)9  —  sin  (to  —  i  —  1)  6] 
+  T  CLs  (sin  {m  +  i-\)6-  sin  (ni  +  t  -  3)  ^ 

+  sin {m-i+^)d-  sin  (m  -  i  +  1)  ^|  +  . .,, ; 
.*.  I   Pj  cos  mO  sin  ^  c?^ 

I^.fi^ ^+-1 ^4 

2  (m  +  z  +  1      TO  + 1  —  1      TO  —  t  +  1      111  —  i—1) 

+  ^M-1 1— +— 1. i.-l  +  ... 

2     (to +  1  —  1        7/1  +  1—3        TO  — t  +  3       TO  — l  +  lj 


TO''-(i+i/  '  to''-(z-1)' 
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Now,  when  m  is  very  large  as  compared  with  %,  this  be- 
comes 


=  -2 


a+c:._2+...     2 


m  m 

since  C^  +  (7,._2  +  ...  =  1,  as  may  be  seen  by  putting  ^  =  0. 

/•t                   ^                                    .    .         2        , 
Hence   /    P,  cos  m^  sin  6  dd  tends  to  the  limit o ,  as  m 

is  indefinitely  increased. 

The  value  of  the  factor  involving  m   has   been  shewn 
above  to  be 

{m  -  (i  -  2)]  {m  -  (t-  4)} ...  (m  -  2)  m'  (m  +  2)  ...  (?n  +  i-  2) 
[m  -  [i  +  1)]  {m  -  (i  -  1)}  ...  (w  -  1)  (m  +  1)  ...  (m  +  1+1) 

if  m  be  even,  and 
{m  -  (/-  2)1  {m  -  (i-  4)1  . . .  (m  -  1)  (m  +  1)  . ..  (m  +  i-  2) 
"{m  -  (iTl) f  [m  -  (t  -  1)1  ...  (m"^2y  (?n  +  2)  ...  (m  +  i+  1) 

if  m  be  odd. 

Each  of  these  factors  contains  in  its  numerator  two  factors 
less  than  in  its  denominator.     It  approaches,  therefore,  when 

m  is  indefinitely  increased,  to  the  value  — ^  .     Hence 


/, 


m 
Pi  cos  m6  sin  6  dO 


0 

(^-(^-2)]{m-(z-4)]...(m-2)m^(m  +  2)...{m+'(/-2)} 
~~    [m-{i+l)][m-{i-\)\...{m-l){m  +  l)...{m+  («  +  l)| 
if  m  and  *  be  even,  and 
_  _    \in-{i-2)]  |m-(i-4)]...(m-l)(m  +  l)...{m+(/-2)} 

[m-  {i+  1)]  [m-  (t-  1)} ...  (m -2) (m+2) ...  [7/l+(^  +  l)j 

if  w  and  i  be  odd. 

In  each  of  these  expressions  {  may  be  any  integer  such 
that  m  —  i  is  even,  *  being  not  greater  than  m.  Hence  they 
will  always  be  negative,  except  when  i  is  e^ual  to  m* 
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20.     We  may  apply  these  expressions  to  develop  cos  mO 
in  a  series  of  zonal  harmonics.. 

Assume 

Multiply  by  Pi  sin  6,  and  integrate  between  the  limits  0 
and  IT,  and  we  get 

|^_(,•-2)}{m-(^-4)}...{m+(^-2)|^      2 
{m-{i+l)][m-{i-\)]...[m-{-{i+l)]     2i  +  l    '' 

Hence 

B=.-(^i^^^  {^^-  (^-  2)}  \m-{{-^)]  ...  {m  +  (/-2)} 
^i         ^  '^^^  i^,n-(i-vl)\[m-{i-l)]...[m  +  {i^l)] 

Hence,  putting  m  successively  =  0,  1,  2,  ...  10, 
cosO^  =  P,; 
cos6  =  Pp 

=^P -^p- 


2 .  4^  6  „  4^^ 

-1.1.3.5.7.9^        +1.3.5.7 


cos  4^  =  -  9 -p^,,^,^^^  P,  -  5 -;-^-,,-^  P, 


-1-^P 

3.5    " 


^64  16  J^ 

35     *     21^      15     *" 


F.  H. 
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Kz,    T-,    2.4.6.8    „   ^   4.G 
cos  5^  =  -  11   -.  -.  o  ^— ,r-T-^  -P«  -  7 


-1.1.3.7.9.11  °   1.3.7.9  " 


-3-^P 


nQ    5    0   8    7   1' 


COS  G^  =  -  13   -.  .  '^  •.  :  ;  ..  .» Pe 


G3^  9  «  7 
2.4. G^ 8. 10 
-1.1.3.5.7.9.11.13 

4.6-. 8  6-       1 

1.3.5.7.9.11  '   3.5.7.9  2  5.7  » 
^512    884    £    J^ 
231  "^   385^  21^  35  »' 

^a-     ^K        2.4.6.8.10.12    ^ 
'''^'^^-"^  -1.1.3.5.9.11.13.15^' 

_     4.6.8.10      y  •  6.8      J^ 
1.8.5.9.11.13^   3.5.9.11  ^^  5.9^ 
1024    l'>8    112     1 
429  '      117^  495^  15^' 

^.    ^^     2.4.G.8M0.12.14 
^^^^^  =  -^^-1.1.3.5.7.9.11.13.15.17^- 

_-,. 4J^8^10J.2 6.8^JlO 

1 .  8 .  5  .  7 .  9  .  11 .  13  .  15  «   3 .  5  . 7 .  9 .  11 .  13  ^ 

5.7.9.11^  7.9  " 

^  16384    4096  p 256     64  p  _  JL  p 

~  6435  «  3465  «  1001^   693  '^  63^' 
cos9^  =  -19    ^•4.6.8.10.12.14.16__^ 
cos  jc;    -^^-1.1.3.5.7.11.13.15.17.19'' 

„ 4  J3^8^  10  .  12  .  14   p  _  -,  -,  __6  .  8  .  10  .  12_  ^, 

1.3.5.7.11 .13 .15  .  17  '  3  .  5.  7 .  llTl37l"5  ^  ' 

_    8^10__   _^^ 
5.7.11.13^   7.11  ' 


COS  10^  =  -21 
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32768   _  3072   _  128   _  16  „   3 
12155  "  2431^  455  '^  143  »  77^* 
2.4.6.8.10M2.14.16.18 


-  1 . 1 .  3 .  5 . 7 .  9  .  11 .  13  .  15 .  17 .  19 .  21 
4.6.8.10M2.14.16 


-17 


1.3.5.7.9.11.13.15.17.19^ 


6.8.10M2.14  8.10M2 

3.5.7.9.11.13.15.17    '        5.7.9.11.13.15^ 

5         10'        p 1     p 

7.9.11.13    '^     9.11    <> 

_  131072   _  32768   ^^12  p  _  128   ^  ^0  „ 
~  46189  ^^  24453^  1683  «  1001  *  9009  '^ 

99  " 

21.  The  present  will  be  a  convenient  opportunity  for 
investigating  the  development  of  sin^  in  a  series  of  zonal 
harmonics.  Since  sin  ^  =  (1  —  /x'^ji,  it  will  be  seen  that  the 
series  must  be  infinite,  and  that  no  zonal  harmonic  of  an  odd 
order  can  enter.     Assume  then 

sin^=  C,P,+  C,P,  +  ...  +  C,P,+  ... 

i  being  any  even  integer. 

Multiplying  by  P^,  and  integrating  with  respect  to  yx 
between  the  limits  —  1  and  +  1,  we  get 

ri  2 

J     PiSmed/j,  =  ^j^Ci; 

2/4-1   fi 

=  ^—j    Pi  sin' Odd, 

supposing  Pi  expressed  in  terms  of  the  cosines  of  6  and  its 
multiples 

2^  + 1  f'^ 
=  rfJlA      F,{1- 00^29)  dd. 

■*     J  ft 
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Hence,  putting  i  =  0, 


4' 

1  .  8 


Putting  i  =  2,  and  observing  tliat  Pg  =  7  +  t  cos  2^, 


a= 


5  f^  (1  +  3  cos  26)  (1  -  cos  2^)  ,- 
4Jo 4 ^^ 

^  f"  |l  +  2  cos  2^  - 1  (1  +  cos  4^)1  d?^ 


32 


For  values  of  *  exceeding  2,  we  observe,  that  if  we  write 
for  Pi  the  expression  investigated  in  Art.  18,  the  only  part 

of  the  expression  |  Pj  (1  —  cos  26)  d6  which  does  not  vanish 

Jo 
will  arise  either  from  the  terms  in  Pj  which  involve  cos  26,  or 
from  those  which  are  independent  of  6.     We  have  therefore 


a  = 


2i  +  l  T'^ 


Jo 


l,3^.(z  +  l)1^3^(z-8) 
2.4...  (z  + 2)  2.4...  (i- 2)  ^^"^-'^ 
,1.3... (t-1 


[2.4...      * 
=  ?i±i    l•3...(^-l)  l■3...(^•-3) 
4      '2.4...      t~  2.4...  (i-2) 


(1  -  cos  26)  d6 


L 


- .—  +  l-'^-l  2  cos  2^)  (1  -  cos  26)  dd 

I         t+2  J^  ^ 


=  2£+l  1.3...  (t-1)  1.3...  (^•-  3)      //-I  _  ^+1 
4      2.4...      i       2.4...  (z- 2)      V    i        t  +  2 


=  —  TT 


2^  +  l  1.3...  (i-l)   1 .3...  (t-3) 

2      2.4...^(^•+2)  2.4...(^-2)^• 


Hence     sin  ^  =  -  P. 
4     " 
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P„-.. 


_  (2t  +  l)7r  1.3...  (/-I)  1 .  3...   (/-3)   p  _ 

2         2.4...f(i  +  2)2.4...(z-2}*     *      *" 
i  being  any  even  integer. 
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dP 

22.     It  will  be  seen  that  -^' ,  being  a  rational  and  integral 

function   of  /a*~S   y"-'*^...,   must   be  expressible  in  terms  of 
■^i_ij  ^i_3«'«     To  determine  this  expression,  assume^ 

dP. 

then  multiplying  by  P„,  and  integrating  with  respect  to  /x 
from  —  1  to  +  1, 

And  jpJ£d>.  =  P^P.~jP,^-^-dj.. 

Now,  since  i  >  m, 

1'''^''-°' 


■■■/: 


Pj^^^d;.  =  [P„AY-[P„Ar=^, 


since  either  m  or  i  must  be  odd,  and  therefore  either  P,„  or 
Pj  =  —  1,  when  /^  =  —  1 ; 

dP 

<ZP      fZP 
Hence  ^p  -  ^'^'  =  (2i  - 1)  P^,. 

dfj.        dfji       ^  ^ 

23.    From  this  equation  we  deduce 

J  ij. 

the  limits  /a  and  1  being  taken,  in  order  that  P^  —  Pi.,  may 
bq  equal  to  0  at  the  superior  limit. 
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Now,  recurring  to  the  fundamental  equation  for  a  zonal 
harmonic,  we  see  that 

24.  We  have  already  seen  that  I  JP^P^  dfi^O,  i  and  m 
being  different  positive  integers.  Suppose  now  that  it  is 
required  to  find  the  value  of  1    P^P^  dfj,. 

"We  have  already  seen  (Art.  10)  that 

j^    *    "•   '^  (^■-m)(^■  +  m^-l) 

And,  from  above, 

n      „2^^-^"»_     wjw  +  l),p     _p    N 
^^"''^"^-       2m +1    ^   '"-^^        "*-^^' 

Jm  (»  — m)  (*+m+l)  [2m  +  l 

_lii±I)p  rP    -P    )1 


1  r^ ( 

IT  Jo  a  — 
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25.  "We  "will  next  proceed  to  give  two  modes  of  ex- 
pressing Zonal  Harmonics,  by  means  of  Definite  Integrals. 
The  two  expressions  are  as  follows : 

P^=-  r{fi±{H>*- 1)^  cos  fy  df. 

These  we  proceed  to  establish. 
Consider  the  equation 

^ 1 

The  only  limitation  upon  the  quantities  denoted  by  a 
and  h  in  this  equation  is  that  6*  should  not  be  greater  than 
a^    For,  if  Jf  be  not  greater  than  a^,  cos  ^  cannot  become 

equal  to  j-  while  ^  increases  from  0  to  tt,  and  therefore  the 

expression  under  the  integral  sign  cannot  become  infinite. 

Supposing  then  that  we  write  z  for  a,  and  V—  1  p  for  6, 
we  get 

•^Jo  3-V~lpcos^     («'  +  p')^* 

We  may  remark,  in  passing,  that 

f "  <^^  _  ['  d^ 

Jo  2  — V— Ipcos^     ^0  ^r  +  V  — Ipcos^ 


Jo 


'o  2'' +  p* cos''' ^  ' 
and  is  therefore  wholly  real. 

Supposing  that  p'  =  £c'  +  2/',  and  that  x^+y^  +  z'  =  r^,v:Q 
thus  obtain 

1  /•'  J^  ^  1 

^io  ^-V-^pcos^     ^  * 
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Differentiate  i  times  with  respect  to  2;,  and  there  results 


*      TT    1 .  2  .3...id^' Jo  5^  — V- Ipcos^ 


TT  Jo    (z-^/'^l 


lj+l 


In  this,  write  fir  for  5r,  and  (1  —  /r)-^  r  for  p,  and  we  get 


TTJo  {/a-OLt'-lpCOS^}' 

which,  writing  tt  —  ^  for  ^,  gives 
1  ^  d'^r 


R  = 


•*       ttJo  {/i+(/i'-l)*COS^f^' 

26.    Again,  we  have 

1  1  l""  __d±__ 

In  this  write  1  —  fxh  for  a,  and  +  (//.^  —  1)^  h  for  &,  which  is 
admissible  for  all  values  of  h  from  0  up  to  fju—  {/jJ'  —  1)^  and 
we  obtain,  since  a^  —  h^  becomes  1  —  2/i^  ^-  K\ 


^ 1   /•"•  cZi/r 


(1  -  2fJLh  +  h')^'      TT  Jo  1  -  /x/t  +  (/a^-  1)*  /i  COSl/r 


^  1  ^ 

TT  Jo   1  —  1/^  +  (/"-^  —  1)^  COS  -v|r}  A  ' 
:  ~  r  cZi|r  [1  +  {/A  ±  (yx'  -  1)^  COsf  I  Ji+... 

+  [fi  ±  [p,'' -1)^  COS  f}'h' +...]. 
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Hence,  equating  coefficients  of  ^*, 

Pt  =  ^j    {fji  ±  {fjr -1)^  cos  fYdylr. 

The  equality  of  the  two  expressions  thus  obtained  for  P^  is 
in  harmony  with  the  fact  to  which  attention  has  already 
been  directed,  that  the  value  of  Pj  is  unaltered  i£  —  (i  +  1) 
be  written  for  i. 


27.  The  equality  of  the  two  definite  integrals  which 
thus  present  themselves  may  be  illustrated  by  the  following 
geometrical  considerations. 


Let  0  be  the  centre  of  a  circle,  radius  a,  C  any  point 
within  the  circle,  PC^  any  chord  drawn  through  C,  and  let 
0  C=  h,  GOP  =  %  COQ  =  ^.  Then  CF'  =  a^  +  b''-  2ab  cos  ^, 
C  Q^  =  a^  +  l/-2ab  cos  y}r.     Hence 

(cr  +  6'  -  2ah  cos  ^)  (a'  +  6'  -  2ah  cos  f)  =  (a"  -  by ; 

sin  ^d^  sin  •x/r  dyjr 


+ 


*  *  a'  +  6'  -  2ab  cos  '^     a^  +  b^  —  2ab  cos  i/r 


=  0. 


42  ZONAL  HARMONICS. 

Again,  since  the  angles    OPG,   OQG  are  equal  to  one 
another, 


ry/i    > 


sin^^ ^  sin  0P(7_  sin  0^(7_  sini/r 
CF  00  00      "TTQ  ' 

sin  ^  sin  yjr 


■  ■  (a-  +  6'  -  2ah  cos^)^     (a'  +  6'  -  2ab  cos  f)^  ' 

whence  — ,  .  ^^ r  =  0. 

(a'  +¥-  2ab  cos  ^)i  ^  {d'  +  b'  -  2ab  cos  f)^ 

In  this,  write  a^  +  Jf  =  /x,  2a&=  +  (/a^— 1)^  which  gives 
a^  —  b^  =  1,  and  we  get 

{/x  ±  0.^-1)4  cos  ^r = -  i^  ±  (^'  - 1)*^  ^^^  ^1^^^- 

We  also  see,  by  reference  to  the  figure,  that  as  ^  in- 
creases from  0  to  TT,  -v^r  diminishes  from  tt  to  0.     Hence 

r ^5_ =  f  {/^  +  ({m'  -  if  cos  -^IrVdylr.  - 

28.  From  the  last  definite  integral,  we  may  obtain  an  ex- 
pansion of  Pj  in  terms  of  cos  6  and  sin  6.  Putting  fi  =  cos  6, 
Ave  get 

1    fT  

^i  =  2Trj   [[cos^+V-lcosf  sin^)}* 

+  [cos  ^  —  V— 1  COS  t/t  sin  6]'']  dyjr 

TT  J  Q  1  .  Z 
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+  (_!)-  i{i-'^)-'{i-^^ri  +  l)  ^^^^  ^,^,„  ^^^^  ^y.,^  ^^.^  ^^^ 
J.  •  ^  •  •  •  ^77& 

Now     J^(costr^f=vr^    2m(2m-2)...2-- 

^(^-l)...(^-2m  +  l)  (2m- 1)  (2m-3)...l 
^^^  1.2... 2m  2m  (2771 -2)... 2 

_ i{i-l)...(i-2m  +  l)  ^ 
(2.4...2mf  ' 

/.  P,=  (cos  ^)'-^^^  (cos  ey-'  (sin  ^)''+  ... 

+(_i)«i(i^^(i^^^ 


CHAPTER  III. 

APPLICATIOX    OF    ZONAL    HABMONICS    TO    THE    THEORY    OF 
ATTRACTION.  REPRESENTATION      OF      DISCONTINUOUS 

FUNCTIONS    BY    SERIES    OF    ZONAL    HARMONICS. 

1.  We  shall,  in  this  chapter,  give  some  applications  of 
Zonal  Harmonics  to  the  determination  of  the  potential  of  a 
solid  of  revolution,  symmetrical  about  an  axis.  When  the 
value  of  this  potential,  at  every  point  of  the  axis,  is  known, 
we  can  obtain,  by  means  of  these  functions,  an  expression 
for  the  potential  at  any  point  which  can  be  reached  from 
the  axis  without  passing  through  the  attracting  mass. 

The  simplest  case  of  this  kind  is  that  in  which  the 
attracting  mass  is  an  uniform  circular  wire,  of  indefinitely 
small  transverse  section. 

Let  c  be  the  radius  of  such  a  wire,  p  its  density,  k  its 
transverse  section.  Then  its  mass,  31,  will  be  equal  to  ^irpch, 
and  if  its  centre  be  taken  as  the  origin,  its  potential  at  any 

M 

point  of  its  axis,  distant  z  from  its  centre,  will  be —c. 

Now,  this  expression  may  be  developed  into  either  of  the 
following  series : 

z\      2z'^2Az'     '"^^     ^       2A...2i     ^^^-^  ■•■]■- ^-)' 

We  must  employ  the  series  (1)  if  z  be  less  than  c,  or 
if  the  attracted  point  lie  within  the  sphere  of  which  the  ring 
is  a  great  circle,  and  the  series  (2)  if  z  be  greater  than  c, 
or  if  the  attracted  point  lie  without  this  sphere. 
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Now,  take  any  point  whose  distance  from  the  centre  is  r, 
and  let  the  inclination  of  this  distance  to  the  axis  of  the 
ring  be  d.  In  accordance  with  the  notation  already  em- 
ployed, let  cos  6  ^  fx.  Then,  the  potential  at  this  point  will 
be  given  by  one  of  the  following  series  : 

1.3.5...(2i-l)       r^'  ] 

+  (-1)  — 2T4Tg:3~^^^c^^"^  ••••]•■ -^^^^ 

r  f "     2    '*r='"^2.4^r*      "• 

1.3. 5.. .(2.'^])       c«         )         , 
+  ^~^^       2.4.G...2i      ^.r-.,  +  .-|...(2). 

For  each  of  these  expressions,  when  substituted  for  F, 
satisfies  the  equation  V^V  =0,  and  they  become  respectively 
equal  to  (1)  and  (2)  when  6  is  put  =  0,  and  consequently 
7'  =  z.  The  expression  (2')  also  vanishes  when  r  is  infinitely 
great,  and  must  therefore  be  employed  for  values  of  r  greater 
than  c,  while  (1')  becomes  equal  to  (2')  when  r  =  c,  and  Avill 
therefore  denote  the  required  potential  for  all  values  of  r 
less  than  c. 

These  expressions  may  be  reduced  to  other  forms  by 
means  of  the  cxj)ressions  investigated  in  Chap.  2,  Art.  25,  viz. 

Pi  =  ^'i  (/x+V?-icos^ycZ^, 

orP.=  ^-  [  {fjL  + J /j^i  COB  ylry^'Ulf. 

Substitute  the  first  of  these  in  (1')  and  (observing  that 
fir  =  z)  we  see  that  it  assumes  the  form 


+  274  c^  -'i''^' 
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which  is  equivalent  to 

Mr  ^ 


Mr 


0  [c^  +  {s+(s*^-r^)icos^ni' 

The  substitution  of  the  last  form  of  P,  in  the  series  (2') 
brings  it  into  the  form 

M  r  f        1 1 c^ 

ttJo  1s  +  (^'-0^cos^     2{^  +  (/_^2)4cos^}' 

.hi '— \d^ 

■*"2.4{2+(3^-Oicos^}^      ••■•]' 

which  is  equivalent  to 

M  r  d'^ 


0    [{z  +  (^-7^)-C0ii'^+cf 

2.  Suppose  next  that  the  attracting  mass  is  a  hollow  shell 
of  uniform  density,  whose  exterior  and  interior  bounding 
surfaces  are  both  surfaces  of  revolution,  their  common  axis 
being  the  axis  of  z.  Let  the  origin  be  taken  within  the 
interior  bounding  surface ;  and  suppose  the  potential,  at  any 
point  of  the  axis  within  this  surface,  to  be 

Then  the  potential  at  any  point  lying  within  the  inner 
bounding  surface  will  be 

A,P,  +  A^P^r  +  A,Pf  +  . . .  +  A,P/  +  . . . 

For  this  expression,  when  substituted  for  V,  satisfies  the 
equation  V^F=0;  it  also  agrees  with  the  given  value  of 
the  potential  for  every  point  of  the  axis,  lying  within  the 
inner  bounding  surface,  and  does  not  become  infinite  at  any 
point  within  that  surface. 

Again,  suppose  the  potential  at  any  point  of  the  axis 
without  the  outer  bounding  surface  to  be 

B,     B,     B,  B, 

z       z        z  z 
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Then  the  potential  at  any  point  lying  without  the  outer 
bounding  surface  will  be 

RP,     5,P,     BJ"^  B,P, 

For  this  expression,  when  substituted  for  V,  satisfies  the 
equation  V^  V=  0  ;  it  also  agrees  with  the  given  value  of  the 
potential  for  every  point  of  the  axis,  lying  without  the  outer 
bounding  surface,  and  it  does  not  become  infinite  at  any 
point  within  that  surface. 

By  the  introduction  of  the  expressions  for  zonal  har- 
monics in  the  form  of  definite  integrals,  it  will  be  found  that 
if  the  value  of  either  of  these  potentials  for  any  point  in  the 
axis  be  denoted  by  <^  (z),  the  corresponding  value  for  any 
other  point,  which  can  be  reached  without  passing  through 
any  portion  of  the  attracting  mass,  will  be 

if  <f>[z  +  {s:'-r')^cos^]d^. 

3.  We  may  next  shew  how  to  obtain,  in  terms  of  a  series  of 
zonal  harmonics,  an  expression  for  the  solid  angle  subtended 
by  a  circle  at  any  point.  We  must  first  prove  the  following 
theorem. 

The  solid  angle,  subtended  by  a  closed  j)lane  cwve  at  any 
2)oint,  is  proportional  to  the  component  attraction  p)erpendicul(ir 
to  the  plane  of  til e  curve,  exercised  upon  the  point  by  a  lamina, 
of  uniform  density  and  thickness,  bounded  by  the  closed  plane 
curve. 

For,  if  dS  be  any  clement  of  such  a  lamina,  r  its  distance 
from  the  attracted  point,  9  the  inclination  of  r  to  the  lino 
perpendicular  to  tlie  piano  of  the  lamina,  the  elementary 
solid  angle  subtended  by  dS  at  the  point  will  be 

d^  cos  e 

r'        • 
And  the  component   attraction  of  the   element  of  the 
lamhia  corresponding  to  dS  in  the  direction  perpendicular 
to  its  plane  will  be 

pIc    .^  cos  a, 
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p  being  tlie  density  of  tlie  lamina,  Tc  its  thickness.  Hence, 
for  this  element,  the  component  attraction  is  to  the  solid 
angle  as  pk  to  1,  and  the  same  relation  holding  for  every 
element  of  the  lamina,  we  see  that  the  component  attraction 
of  the  whole  lamina  is  to  the  solid  angle  subtended  by  the 
whole  curve  as  pA;  to  1. 

Now,  if  the  plane  of  xy  be  taken  parallel  to  the  plane 

of  the  lamina,  and    V  be  the  potential  of  the  lamina,  its 

component   attraction  perpendicular    to    its  plane    will   be 

dV 
—  -J— .      Now  since  F  is  a  potential  we  have  V^  F=  0,  whence 

— -V^F=0,  or  VM--T-)=0.     Hence  -j-  is  itself  a  potential, 

and  satisfies  all  the  analytical  conditions  to  which  a  potential 
is  subject.  It  follows  that,  if  the  solid  angle  subtended  by 
a  closed  plane  curve  at  any  point  {x,  y,  z)  be  denoted  by 
&),  ft)  will  be  a  function  of  x,  y,  z,  satisfying  the  equation 
V^ft)  =  0.  Hence,  if  the  closed  plane  curve  be  a  circle  it 
follows  that  the  magnitude  of  the  solid  angle  which  it  sub- 
tends at  any  point  may  be  obtained  by  first  determining 
the  solid  angle  which  it  subtends  at  any  point  of  a  line 
drawn  through  its  centre  perpendicular  to  its  plane,  and 
then  deducing  the  general  expression  by  the  employment 
of  zonal  harmonics. 

Now  let  0  be  the  centre  of  the  circle,  Q  any  point  on  the 
line  drawn  through  0  perpendicular  to  the  plane  of  the 
circle,  E  any  point  in  the  circumference  of  the  circle.  With 
centre  Q,  and  radius  QO,  describe  a  circle,  cutting  QE  in  L. 
From  L  draw  LN,  perpendicular  to  Q  0, 

Let  OE=c,  OQ  =  z. 

Then^Z  =  (c^  +  /)i-^,  ON^^-^—^^ii^  +  z^f^-z] 


{c'  +  z')^' 
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And  the  solid  angle  subtended  by  the  circle  at  0 
A     ON 

=  27r  ll  - 


To  obtain  the  general  expression  for  the  solid  angle  sub- 
tended at  any  point,  distant  r  from  the  centre,  we  first 
develope  this  expression  in  a  converging  series,  proceeding 
by  powers  of  z.     This  will  be 

if  2  be  less  than  c,  and 

"""llz'     2Az*^"'     ^     ^^        2A...2i      5--'+-j 
if  z  be  greater  than  c. 

Hence,  by  similar  reasoning  to  that  already  employed, 
we  get,  for  the  solid  angle  subtended  at  a  point  distant  r 
from  the  centre. 


n       c    ^2    c'       2.4    c''   ^ 


_  1.3. ..(2/- 1)^,7'-^ 

^     ^        2.4...2i  c-"^^     "*" 

if  r  be  less  than  c,  and 

[2  r'  2.4  1-'  ■^••"  ^  ^  2.4...21  r''  '^•••j 
if  ?*  be  greater  than  c. 

4.  We  may  deduce  from  this,  expressions  for  the  potential 
of  a  circular  lamina,  of  uniform  thickness  and  density,  at  an 
external  point.  For  we  sec  that,  if  V  be  the  potential  of 
such  a  lamina,  Ic  its  thickness,  and  p  its  density,  we  have  for 
a  point  on  the  axis. 


F.  H. 
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whence  F=  ^irph  [{c"  +  z'f^  -  z] 

if  21  be  the  mass  of  the  lamina. 

Expanding  this  in  a  converging  series,  we  get 

^~c'f     ^"^2  0      2.40^  +  2.4760^     •••" 

^       ^       2.4.G...21      c''-''^" 
if  z  be  less  than  c,  and 

F=    ,   s--^-T-,+ 


cM2s      2.4  s'  '  2.4.  G^' 


1.1.3...(2t-.S)   c''' 
^     ^       2.4. 6... 2*      s^'-'"^ 


if  s  be  cfreater  than  c. 


Hence  we  obtain  the  following  expressions  for  the  po- 
tential of  an  uniform  circular  lamina  at  a  point  distant  r 
from  the  centre  of  the  lamina  : 

^^  M(j,     „    ,  IP/   i.ipy 

_  r-  IV  i-i-S-.-C^^-s)  P,£-       1 

^     ^       2.4.6...2i        c'*-'  J 

if  r  be  less  than  c,  and 

Jf  n  P/  _  1  ^1  p/    1  j.^3  P,c^  _ 

c"M2    r       2.4  r'    "^2.4.6    r'       '" 


_  1.1.3...(2/-3)  ^5^' 

^     ^       2.4. 6. ..2i         r^'"^   "^ 


if  r  be  greater  tban  c. 
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It  may  be  shewn  that  the  solid  angle  may  be  expressed 
in  the  form 

and  the  potential  of  the  lamina  in  the  form 


C     IT 


r  [c^  +[z  +  (s^ - rf  cos 6>}^]i d9-^. 
Jo  c 


5.  As  another  example,  let  it  be  required  to  determine 
the  potential  of  a  solid  sphere,  whose  density  varies  inversely 
as  the  fifth  power  of  the  distance  from  a  given  external  point 
0  at  any  point  of  its  mass. 

It  is  proved  by  the  method  of  inversion  (see  Thomson 
and  Tait's  Natural  Philosophy,  Vol.  1,  Art.  518)  that  the 

potential  at  any  external  point  P'  will  be  equal  to  — ppr,  0' 

being  the  image  of  0  in  the  surface  of  the  sphere,  and  M 
the  mass  of  the  sphere.  We  shall  avail  ourselves  of  this 
result  to  determine  the  potential  at  a  given  internal  point. 

Let  G  be  the  centre  of  the  sphere,  0  the  given  external 
point.  Join  CO,  and  let  it  cut  the  surface  of  the  sphere  in  A, 
and  in  GA  take  a  point  0', such  that  CO.GO  =  GA\  Then 
0'  is  the  image  of  0. 

Let  P  be  any  point  in  the  body  of  the  sphere,  then  we 
wish  to  find  the  potential  of  the  sphere  at  P. 

Take  0  as  pole,  and  DC  as  prime  radius,  let  OP  =  r, 
POC  =  0.    Also  let  CA  =  a,  G0  =  c. 

Let  the  density  of  the  sphere  at  its  centre  be  p,  then  its 
density  at  P  will  be  p  -^ .     Hence 


i/=27r[[ppr'sin^cZr(Z^, 


4—2 
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the  limits  of  r  being  the  two  values  of  r  which  satisfy  the 
equation  of  the  surface  of  the  sphere,  viz. 

7^  -\-c^  —  2cr  cos  6  =  a^, 

and  those  of  6  being  0  and  sin~^  - . 

c 

Hence,  if  r^ ,  r^  be  the  two  limiting  values  of  r,  we  have 


Now 


r'^      r'^      c'  —  d'   W„     r 


, ,  1       1      2c  cos  0 

Also  -  +  -  =  -^~    ,  , 

r.      r,       c  —  a 


1  _l\^^[c'co^^e-{e-a')Y 


r„      rj  c  —a 

(a' -  c' sin' 6)' 


=  2 


c  —  a 


n 

.:  M  =  -!, ., .  -vr      -.  cos  ^  sni  6  (a  -  c  sm' 0)-  iW 

I      c"  —  a'   c  —  a"  J  ^ 

a 

=  7-2"^  2r2  cos  6  sm  ^  (a"  -  c"  sm''  ^)  ^  c?^ 

(c  —  a  )  j  0 

_  4         TT/OC*  3 

Now,  if  V  be  the  potential  at  P,  Ave  have  (see  Chap.  i. 
Art.  1) 

d'irV)         1      d    f  .    ^dV\         47rpc' 
cZr         sm  6  dO  \  dO  J  r^ 
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This  is  satisfied  by  V=  —  ^  — ^ . 

Assume  then,  as  the  complete  solution  of  the  equation, 

+  {a/+^)p,  +  .... 

It  remains  to  determine  the  coefficients  A^^  A^...A^...B^, 
B^...B^,  so  that  this  expression  may  not  become  infinite  lor 
any  value  of  r  corresponding  to  a  point  within  the  sjDhere, 
and  that  at  any  point  P  on  the  surface  of  the  sphere  it  may 

be  equal  to  j-j-.j,  where  O'P :  OP  ::  a  :  c,  and  therefore,  at 

the  surface, 

p._  Mc    1    _  4     7rpc''a^ 

a    UP      3  (c^—  a'j'  r' 


And,  at  the  surface,  we  have 
'/■*  -  2cr  cos  d  +  c"  —  a'  =  0; 

1  _      1_    /     1  .  2c  cos  ^N 

;•'      0"  —  a' 


/_1      2c  cos  ^\ 


^f-^  +  2o5); 


c  —  a  \      r  r 


+  (A^r+   'Ap^+  ...  identicalbj. 
Hence  ^  ^^'"^  =  C^   '^^''-  +b\p 

V     3  c'-  a'         7     ' 
and  i?,,  i?3,...5....1,,  J^....l.  all  =  0. 
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Hence  since  P(,  =  1, 
and  R  =  ^  "^^^ 


whence  we  obtain,  as  the  expression  for  the  potential  at  any 
internal  point, 

j^_  2     irpd'      Sa^  —  c^      4   Trpc'^    cos  6      2  irpc^ 
3  (c^  -  ay       V~  "^  3  c^^^'  'V'~  ~  3  'r' ' ' 

6.  We  shall  next  proceed  to  establish  the  proposition  that 
if  the  density  of  a  spherical  shell,  of  indefinitely  small  thick- 
ness, he  a  zonal  surface  harmonic,  its  potential  at  any  internal 
jwint  ivill  he  proportional  to  the  corresponding  solid  har- 
monic of  positive  degree,  and  its  potential  at  any  external 
point  tuill  he  proportional  to  the  corresponding  solid  liarmonic 
of  negative  degree. 

Take  the  centre  of  the  sphere  as  origin,  and  the  axis  of 
the  system  of  zonal  harmonics  as  the  axis  of  z.  Let  h  be  the 
radius  of  the  sphere,  hh  its  thickness,  U  its  volume,  so  that 
IJ=  ^Trb^Sh.  Let  (7Pf  be  the  density  of  the  sphere,  P^  being 
the  zonal  surface  harmonic  of  the  degree  i,  and  C  any  con- 
stant. 

Draw  two  planes  cutting  the  sphere  perpendicular  to  the 
axis  of  z,  at  distances  from  the  centre  equal  to  ^,  ^-\-  d^ 
respectively.     The  volume  of  the  strip  of  the  sphere  inter- 

cepted  between  these  planes  will  be  ,-y.   JJ,  and  its  mass  will  be 

CP,  If  U. 

*  26 

Now  ^=  &/A,  hence  d^=hd^,  and  this  mass  becomes 
CTT 
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Hence  the  potential  of  this  strip  at  a  point  on  the  axis  of  e, 
distant  z  from  the  centre,  will  be 

which  may  be  expanded  into 

and       ^5(p„  +  P,^+...  +  P,^;+....)cZ/.if.>J. 

To  obtain  the  potential  of  the  whole  shell,  we  must  inte- 
grate these  expressions  with  respect  to  fx  between  the  limits 
—  1  and  + 1.  Hence  by  the  fundamental  property  of  Zonal 
Harmonics,  proved  in  Chap.  II.  Art.  10,  we  get  for  the  po- 
tential of  the  whole  shell 

CU    s* 
TT^ — -  jjir,  at  an  internal  point, 

„  .  — r  -TXT  at  an  external  point. 

From  these  expressions  for  the  potential  at  a  point  on 
the  axis  we  deduce,  by  the  method  of  Art.  1  of  the  present 
Chapter,  the  following  expressions  for  the  potential  at  any 
point  whatever : 

CU  Pr^ 

V.  =  T^ — q  T^  at  an  internal  point, 

„        CU  Ph'    ^  ,        1      •  . 

y,  =  ^r^ — --t'  —rr  at  an  external  point. 
^     2^  +  1  r^  ^ 

From  hence  we  deduce  the  following  expressions  for  the 
normal  component  of  the  attraction  on  the  point. 

Normal  component  of  the  attraction  on  an  internal  point, 
measured  towards  the  centre  of  the  sphere, 

__fZT/;_   _i_      pr^ 
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Normal  component  of  the  attraction  on  an  external  point, 
measured  towards  the  sphere, 

In  the  immediate  neighbourhood  of  the  sphere,  where  r  is 
indefinitely  nearly  equal  to  b,  these  normal  component  at- 
tractions become  respectively 

and  their  difference  is  therefore 

P. 

And  writing  for  U  its  value,  47r7/S3,  this  expression  be- 
comes 

47rS& .  CP, 

Or,  the  density  may  be  obtained  by  dividing  the  alge- 
braic sum  of  the  normal  component  attractions  on  two  points, 
one  external  and  the  other  intej^nal,  indefinitely  near  the 
sphere,  and  situated  on  the  same  normal,  by  47r  x  thickness 
of  the  shell. 

7.  It  follows  from  this  that  if  the  density  of  a  spherical 
shell  be  expressed  by  the  series 

c,p,  +  c,p,+  a,p,  +  ...  +  cA+-, 

C\,  (7j,  C^  ...  (7; ...  being  any  constants,  its  potential  (FJ  at 
an  internal  point  will  be 

\  b  ^3  b'  +5    b'   ^■"^2i+i  r^  ^--v 

and  its  potential  (  FJ  at  an  external  point  will  be 

/(7„P     10,P,^     1(7„P,6^  1      C.P^h^         \ 

In  the  last  two  Articles,  by  the  word  "density"  is  meant 
"volume  density,"  i.e.  the  mass  of  an  indefinitely  small 
clement  of  tht3  attracting  sphere,  divided  by  the  volume  of 
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the  same  element.  The  product  of  the  volume  density  of 
any  element  of  the  shell,  into  the  thickness  of  the  shell  in 
the  neighbourhood  of  that  element,  is  called  "surface  den- 
sity." We  see  from  the  above  that,  if  the  surface  density 
be  expressed  by  the  series 

the  potentials  at  an  internal  and  an  external  point  will  seve- 
rally be 


/^,P       l^£r      1^;^^  J^    a,P^ 


\   r         6     r         or  2z  + 1    r 

This  variation  in  surface  density  may  be  obtained  eitlier 
by  combining  a  variable  volume  density  with  an  uniform 
thickness,  as  we  have  supposed,  or  by  combining  a  variable 
thickness  with  a  uniform  volume  density,  or  by  varying  both 
thickness  and  density. 

8.  "We  have  seen,  in  Chap,  ii.,  that  any  positive  integral 
power  of  /Lt,  and  therefore  of  course  any  rational  integral 
function  of  fi,  may  be  expressed  by  a  finite  series  of  zonal 
harmonics.  It  follows,  therefore,  that  we  can  determine  tho 
potential  of  a  spherical  shell,  whose  density  is  any  rational 
integral  function  of  fi. 

Suppose,  for  instance,  we  have  a  shell  whoso  density 
varies  as  the  square  of  the  distance  from  a  diametral  plane. 
Taking  this  plane  as  that  of  xi/,  the  density  may  be  ex- 
pressed by  pfi',  or  p  ^ .  We  have  seen  (Chaj).  ii.  Art.  20) 
that 

;.»  =  |(1+2PJ. 

Hence,  by  the  result  of  the  last  Article,  the  jDotential 
will  be 

U/1     2Py\    ,        .  ,        ,      . 
p  o  [t  "^r    fa")  ^^  ^^  internal  point. 
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Ufl     2Ph'\    ^  .        1      .  . 

/3  o  (  -  +  V  ^- J  ^^  ^^  external  point ; 

or,  since  F^r'  =  — - —  r  =  — ^ — ,  we  obtain 

77"    /I  1     0„2  2\ 

/>  «-  (r  +  V  — p — J  for  the  potential  at  an  internal  point, 

Ptt  j-+v-(— :;3+-;5-)rfor  that  at  an  external  point. 

9.  As  an  example  of  the  case  in  which  the  density  is  re- 
presented by  an  infinite  series  of  zonal  harmonics,  suppose  we 
wish  to  investigate  the  potential  of  a  spherical  shell,  whose 
density  varies  as  the  distance  from  a  diameter.  Taking  this 
diameter  as  the  axis  of  z,  the  density  will  be  represented 
by  p  sin  6,  or  p  (1  —  fji^^y-^.  We  have  investigated  in  Chap.  II. 
Art.  21,  the  expansion  of  sin  6  in  an  infinite  series  of  zonal 
harmonics.  Emplojang  this  expansion,  we  shall  obtain  for 
the  potential 

2   6  J2     »     16^6^     ■■'    2.4...^(^■+2)■2.4...(^•-2>     *6*     "**]  ' 

or 

T    TTi^P^o_l_pl 1.3■■.(^-l)    1.3.. .^-3)        V         | 

2  ^  |2  r  IG  V  •"  2.4...^(^•+2)•2A..(^■-2)^  ^'/^^  ■••}' 
according  as  the  attracted  point  is  internal  or  external  to  the 
spherical  shell,  i  being  any  even  integer.  All  these  expres- 
sions may  be  obtained  in  terms  of  surface  density,  by  writing, 
instead  oi  pU,  47rcV. 

10.  We  may  next  proceed  to  shew  how  the  potential  of 
a  spherical  shell  of  finite  thickness,  whose  density  is  any  solid 
zonal  harmonic,  may  be  determined.  Suppose,  for  instance, 
that  we  have  a  shell  of  external  radius  a,  and  internal  radius 
a ,   whose   density,   at   the   distance  c  from  the  centre,   is 

)  i  Pfi^i  '^  being  any  line  of  constant  length. 

Dividing  the  sphere  into  concentric  thin  spherical  shells, 
of  thickness  dc,  the  potential  of  any  one  of  these  shells,  of 
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radius  c,  at  an  internal  point  distant  r  from  the  centre  "will 

DC* 

he  obtained  by  writing  c  for  6,  jr-  for  C,  ^ttcWc  for  U,  in 
tbe  first  result  of  Art.  6.     This  gives 

To  obtain  the  potential  of  the  whole  shell,  we  must  inte- 
grate this  expression,  with  respect  to  c,  between  the  limits 
a  and  a.     This  gives 

Again,  the  potential  of  the  shell  of  radius  c,  at  an  external 
point,  will  be 

p  47rcVc  P,  c«        _47r    ^      c«*''  , 
/i'  2*  +  1    r'^'   ^^  2i  +  1  h'    '  7-'"^' 

Integrating  as  before,  we  obtain  for  the  potential  of  the 
whole  shell, 

47r   _      p^       (oT^-ar^ 

Suppose  now  that  we  wish  to  obtain  the  potential  of  the 
whole  shell  at  a  point  forming  a  part  of  its  mass,  distant  r 
from  the  centre.  We  shall  obtain  this  by  considering  sepa- 
rately the  two  shells  into  which  it  may  be  divided,  the 
external  radius  of  the  one,  and  the  internal  radius  of  the 
«jther,  beiug  each  r.  Writing  r  for  a,  in  the  first  of  the  fore- 
going results,  we  obtain 

o  • 7  S^i    (*    —  '>  )  '>  ' 

2l  +  l     h'    ^  ^ 

And  writing  r  for  a  in  the  other  result,  we  obtain 

47r pP,  r^-a'^ 

(2i+l)(2i  +  3)  h'         r'^'       ' 

Adding  these,  we  get  for  the  potential  of  the  whole 
spliere 

47r     pPi((^-jr^   t       r^'^s  _  ^2.«  I 

2c-\- 1  h'   Y2~  *'  "^  (2iT3y7"'j  • 
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It  is  hardly  necessary  to  observe  that  the  corresponding 
results  for  a  solid  sphere  may  be  obtained  from  the  foregoing, 
by  putting  a'  =  0. 

If  the  density,  instead  of  being  ~  P.  c\  be  ^,  T^  c™,  similar 

reasoning  will  give  us,  for  the  potential  of  the  thin  shell  of 
radius  c  and  thickness  dc  at  an  internal  and  external  point 
respectively. 

And,  integrating  as  before,  we  obtain  for  the  potential  of 
the  whole  shell, 

.p.y»-+2_^'m-it2j  ,.i  at  an  internal  point, 


{2l+l){m-i  +  2)k' 


4<7r  p   ^  a 


m+i+3 


7y>-  I  T\r"  I   •  ,  o\T^Pi~ i+i- ^t  ^^  external  point. 

(2z  + 1)  {m  +  1  +  3)  h      *  9''^  ^ 

And,  at  a  point  forming  a  part  of  the  mass, 


21  +  1  /r  \   m-i+2  m  +  i  +  '^    r" 

11.  Suppose,  for  example,  that  we  wish  to  determine,  in 
each  of  the  three  cases,  the  potential  of  a  spherical  shell 
whose  external  and  internal  radii  are  a,  a,  respectively,  and 
Avhose  density  varies  as  the  square  of  the  distance  from  a 
diametral  plane. 

Taking  this  plane  as  that  of  xij,  the  density  may  be  ex- 
pressed by  pz^,  or  p^c'/ji^     Now  yu,^  =  ^^-^- — .      Hence  the 
density  of  this  sphere  may  be  expressed  as 
^^Pc'^  +  llPr' 

The  several  potentials  due  to  the  former  term  will  be, 

2 
writing  2  for  i  and  multiplying  by  - , 

o 


TO   THE   THEORY   OF  ATTRACTION.  Gl 

And  for  the  latter  term,  writing  0  for  i,  and  2  for  m,  and 
multiplying  by -g, 

.47r  p      ,       „      477  p  a^-  a^     jir  p_  fa*-r*     /  -  a"\ 
12  ^^^  ^"^  ~  ^'  ^'    15  ^-^       r      '      3  An    4      "^     5r    J ' 

And,  since  Py  =  — -^ — ,  we  get  for  the  potential  at  an 
internal  point 

at  an  external  point 

p  f  47r  fJ  -  a'  47r  a'  -a"]  , 

^  po  ~7'>  ~  ^'^^  "'  ''  +  15  "r     J  ' 

at  a  point  forming  a  part  of  the  mass 


115"  l-2~  +  -1^'^)  (3— »  ^  +  -3-  1,-4"  +  "5r^;j  • 


12.  We  may  now  prove  that  by  means  of  an  infinite  series 
of  zonal  harmonics  w^e  may  express  any  function  of  //.  what- 
ever, even  a  discontinuous  function.  Suppose,  for  instance, 
that  we  wish  to  express  a  function  which  shall  be  equal  to 
A  from  /i-  =  l  to  fi  =  X,  and  to  B  from  fi=\  to  /m=  —  1. 
Consider  what  will  be  the  potential  of  a  spherical  shell, 
radius  c,  of  uniform  thickness,  whose  density  is  equal  to  A 
for  the  part  corresponding  to  values  of  fj,  between  1  and  X, 
and  to  B  for  the  part  corresponding  to  values  of  p,  between  \ 
and  —  1. 

Divide  the  shell,  as  before,  into  indefinitely  narrow  strips 
by  parallel  planes,  the  distance  between  any  two  successive 
planes  being  cdp.. 
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We  have  then,  for  the  potential  of  such  a  sphere  at  any 
point  of  the  axis,  distant  z  from  the  centre, 

for  the  first  part  of  the  sphere 


2'7rAc''Sc 


and  for  the  latter  part 

■T 


27rBc'Bc  '  ^^ 


These  are  respectively  equal  to 
27r^c'Sc 


at  an  internal  point ;  and  to 
^ttAc'Bc 


/;(p.+p.£+...+p,^,+....)<?^. 


z 


z 

at  an  external  point. 

Now  it  follows  from  Chap.  II.  (Art.  23)  that  if  i  be  any 
positive  integer, 

£  ^^^^  =  -  2^  +  1  ^^'-^^^  "  ^-^^^^^' 
whence,  since  1   F^d^i.  =  0,  it  follows  that 
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Also  f  P^dii  =  l-\  I  F,dfj,  =  l+X. 

Hence  the  above  expressions  severally  become  : 
For  the  potential  at  an  internal  point  on  the  axis 

-?^"  [^  (1  -  X)  +  5  (1  +  X)  -  ^  {P,(X)  -p.w  j  J 
and  for  the  potential  at  an  external  point  on  the  axis 


^[P^.{^)-P.-M]-^.--'] 


2i  + 

Hence  the  potentials  at  a  point  situated  anywhere  are 
respectively 


27rc^ 
c 


[M(1-X)  +  J?(1  +  X))P» 


at  an  internal  point; 


and 
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at  an  external  point. 

Now,  if  we  inquire  what  will  be  the  potential  for  the 
following  distribution  of  density, 

^,[.4(1  -  X)  +  B(l  +  X)  -  iA-B)[P^{X)  -  Po(X)lP,(/x) 

-  (.4- P)[P3(X)-P,(X)1P,  (/.)-..  . 

we  see  by  Art.  C  that  it  will  be  exactly  the  same,  both  at 
an  internal  and  for  an  external  point,  as  that  above  in- 
vestigated for  the  shell  made  up  of  tvfo  parts,  whose  densities 
are  A  and  B  respectively. 

But  it  is  known  that  there  is  one,  and  only  one,  dis- 
tribution of  attracting  matter  over  a  given  surface,  which 
Avill  produce  a  specified  potential  at  every  point,  both  ex- 
ternal and  internal.  Hence  the  above  expression  must 
represent  exactly  the  same  distribution  of  density.  That  is, 
writing  the  above  series  in  a  slightly  different  form,  the 
expression 
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is  equal  to  A,  for  all  values  of  fi  from  1  to  \  and  to  B  for  all 
values  of  /*  from  \  to  —  1. 

13.  By  a  similar  process,  any  other  discontinuous  function, 
whose  values  are  given  for  all  values  of  f/,  from  1  to  —  1,  may 
be  expressed.  Suppose,  for  instance,  we  wish  to  express  a 
function  which  is  equal  to  A  from  fj,  =  1  to  /u.  =  X^,  to  B  from 
fjb  =  \  to  /M  =  X^,  and  to  C  from  fji.=X^to  {i  =  —  l.  This  will 
be  obtained  by  adding  the  two  series 

For  the  former  is  equal  to  A  —  B  from  /a  =  1  to  /i  =  Xj, 
and  to  0  from  //.  =  X^  to  /i  =  —  1 ;  and  the  latter  is  equal  to 
B  from  ^  =  1  to  /x  =  Xg,  and  to  C  from  /u,  =  Xj,  to  /tt  =—  1. 

By  supposing  A  and  G  each  =  0,  and  5  =  1,  we  deduce  a 
series  which  is  equal  to  1  for  all  values  of  /x  from  /a  =  X^  to 
fi  =  Xj,  and  zero  for  all  other  values.     This  will  be 

1 


U\-\+{P^^^^-PM  -  Po{\)-Po{\)]P.iH')  +  '- 

+ {p,.,(xj-p,,,(\)  -  p,.,(xj-pz:(xj}p,  (^) +. . .]. 

This  may  be  verified  by  direct  investigation  of  the 
potential  of  the  portion  of  a  homogeneous  spherical  shell, 
of  density  unity,  comprised  between  two  parallel  planes, 
distant  respectively  cX^  and  cX^  from  the  centre  of  the 
spherical  shell. 

14.  In  the  case  in  which  X^  and  X^  are  indefinitely  nearly 
equal  to  each  other,  let  X^^  =  X,  and  X^  =  X  +  cZX.  We  then 
have,  ultimately, 

F.  H.  5 
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Hence   P,,A)  "^mW  ""^i-iW -^i-x(^) 

\    d\  dX     ) 

=  (2^•  +  l)P.(\)c^\. 
Hence  the  series 

'III  +  2,PX\)PM  +  5P.wn(iu)  + ... 

+  (2^  +  l)P,(X)PM+...} 
is  equal  to  1  when  /i  =  X  (or,  more  strictly,  when  jm  has  a«iy 
value  from  X  to  X  +  dX)  and  is  equal  to  0  for  all  other  values 
of /i. 

We  hence  infer  that 

1  +  3P,(X)P,(/.)  +  ...  +  (2;  +  l)P,(X)P,(yc.)  +  ... 
is  infinite  when  /i  =  X,  and  zero  for  all  other  values  of  /x. 
15.    Representing  the  series 
Ml  +  3P,(X)P,0a)  +  ...  +  (2^4-  1)P,(X)P,(/.)  +  ...}:' 

by  ^(X)  for  the  moment,  we  see  that  pj>QC)d\  is  equal  to  p 
when  fM  =  \,  and  to  zero  for  all  other  values.  Hence  the 
expression 

is  equal  to  p^  when  /*  =  X^,  to  p^  when  fi  =  \...  Supposing 
now  that  X^,  X^...  are  a  series  of  values  varying  continuously 
from  1  to  —  1,  we  see  that  this  expression  becomes 


/: 


p0(X)(ZX,- 


p  being  any  function   of  X,  continuotis"  or   discontinuous. 
Hence,  writing  ^(X)  at  length,  we  see  that 

l^f^pdX  +  SP^{fi)j'  pP^{X)d\+... 

+  {2i  +  l)P,{fi)j'  pP/X)JX  +  ...l 

is  equal,  for  all  values  of  /x  from  -  1  to  + 1,  to  the  same 
function  of  /j,  that  p  is  of  X. 
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16.    The  same  conclusion  may  be  arrived  at  as  follows : 

The  potential  of  a  spherical  shell,  whose  density  is  p, 
and  volume  U,  at  any  point  on  the  axis  of  z,  is  * 

which  is  equal  to ^-j  I    pd\ -\ —  I   pPiQ^) d\+  ... 

for  an  internal  point, 

and  to   ^\~  f  pd\  +  4  f  pP,{X)dX  +  ... 
^K^J-i  ^  J  _i 

for  an  external  point. 

It  hence  follows  that  the  potential,  at  a  point  situated 
anywhere,  is 

for  an  internal  point, 

andto    j[lj[pd\+?^'f^pP,{X)dX  +  ... 

for  an  external  point. 

And  these  expressions  are  respectively  equal  to  those 
for  the  potentials,  at  an  internal  and  external  point  re- 
spectively, for  matter  distributed  according  to  the  following 
law  of  density : 

5—2 
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+  {2i  +  l)P,{/x)j'  pP,{\)d\+  ..j. 

It  will  be  observed,  in  applying  tbis  formula,  tbat  if  p  be 
a  discontinuous  function  of  X,  eacb  of  tbe  expressions  of  tbe 

form  I  pPi{X)d\  will  be  the  sum  of  tbe  results  of  a  series  of 

integrations,  eacb  integration  being  taken  through  a  series  of 
values  of  X,  for  which  p  varies  continuously. 


CHAPTER  IV. 

SPHEEICAL  HARMONICS  IN  GENERAL.  TESSERAL  AND  SEC- 
TORIAL HARMONICS.  ZONAL  HARMONICS  WITH  THEIR 
AXIS  IN  ANY  POSITION.  POTENTIAL  OF  A  SOLID  NEARLY 
SPHERICAL  IN   FORM. 

1.  We  have  hitherto  discussed  those  solutions  of  the 
equation  VT=  0  which  are  symmetrical  about  the  axis  of  z, 
or  in  other  words,  those  solutions  of  the  equivalent  equation  in 
polar  co-ordinates  which  are  independent  of  ^.  We  propose, 
in  the  present  Chapter,  to  consider  the  forms  of  spherical 
harmonics  in  general,  understanding  by  a  Solid  Spherical 
Harmonic  of  the  i**  degree  a  rational  integral  homogeneous 
function  of  x,  y,  z,  of  the  i^^  degree  which  satisfies  the  equa- 
tion V*  V=  0,  and  by  a  Surface  Spherical  Harmonic  of  the 
i^  degree  the  quotient  obtained  by  dividing  a  Solid  Sphe- 

rical  Harmonic  by  («*  -{-y^-^-  z^'\  Such  an  expression,  as  we 
see  by  writing  a;  =  r  sin  ^  cos  <^,  y  =  r  sin  ^  sin  ^,  z  =  r  cos  6, 
will  be  of  the  i^  degree  in  sin  6  cos  ^,  sin  d  sin  <^,  cos  9 ;  and 
will  satisfy  the  differential  equation  in  Y^, 

1     d  f  .    ^dY\         1     ffZ      .,.     ,,^^     ^ 

sm  6  dd\  dd  J     &m  6  d(f>^        ^        /     i       ' 

or,  writing  fi  for  cos  0, 

It  will  be  convenient,  before  proceeding  to  investigate  the 
algebraical  forms  of  these  expressions,  to  discuss  some  of 
their  simpler  physical  properties. 

2.  We  will  then  proceed  to  shew  how  spherical  har- 
monics may  be  employed  to  determine  the  potential,   and 
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consequently  the  attraction,  of  a  spherical  shell  of  indefinitely 
small  thickness. 

"We  will  first  establish  an  important  theorem,  connecting 
the  potential  of  such  a  shell  on  an  external  point  with  that 
on  a  corresponding  internal  point.     The  theorem  is  as  follows : 

If  0  he  the  centre  of  such  a  shell,  c  its  radius,  P  any  in- 
ternal point,  P'  an  external  point,  so  situated  that  P'  lies  on 
OP  produced,  and  that  OP .  OP'  =  c\  and  if  OP  =  r,  OP'  =  r', 
then  the  potential  of  the  shell  at  P  is  to  its  potential  at  P' 
as  c  to  r,  or  {luhich  is  the  same  thing)  as  r'  to  c. 

For,  let  A  be  the  point  where  OP'  meets  the  surface  of 
the  sphere,  Q  any  other  point  of  its  surface.  Then,  by  a 
known  geometrical  theorem, 

QP:  QF  v.AP'.AF  .:c-r:r'~c. 

c  —  r      cr  —  r"^      cr  —  r^     r      c 
r  —  c     rr  —  cr     c''  —  cr     c      r  ' 

Again,  considering  the  element  of  the  shell  in  the  im- 
mediate neighbourhood  of  Q,  its  potential  at  P  is  to  its 
potential  at  P'  as  QP  is  to  QP,  that  is,  as  c  to  r,  or  (which 
is  the  same  thing)  as  r'  to  c,  which  ratio,  being  independent 
of  the  position  of  Q,  must  be  true  for  every  element  of  the 
spherical  shell,  and  therefore  for  the  whole  shell.  Hence 
the  proposition  is  proved. 

3.     Now,  suppose  the  law  of  density  of  the  shell  to  be 

such  that  its  potential  at  any  internal  point  is  F  {jx,  (p)  -^ . 

c 

Then  FQu,,  ^)  —.  must  be  a  solid  harmonic  of  the  degree  i. 
c 

Hence  F  (jm,  ^)  must  be  a  surface  harmonic  of  the  degree  t. 
Let  us  represent  it  by  Y^. 

By  the  proposition  just  proved,  the  potential  at  any 
external  point,  distant  r  from  the  centre,  must  be 

Y  ~  ^- 
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Hence,  the  component  of  tlie  attraction  of  the  sphere  on 
the  internal  point  measured  in  the  direction  from  the  point 
inwards,  i.e.  towards  the  centre  of  the  sphere,  is 

c* 

And  the  component  in  the  same  direption  pf  ,the  attraction 
on  the  external  point,  measured  inwards,  is 

Now  suppose  the  two  points  to  lie  on  the  same  line 
passing  through  the  centre  of  the  sphere,  and  to  be  both 
indefinitely  clo^e  to  the  surface  of  the  sphere,  so  that  r  and  r 
are  each  indefinitely  nearly  equal  to  c. 

And  the  attraction  on  the  external  point  exceeds  the 
attraction  on  the  internal  point  by 

0 

Now,  supposing  the  shell  to  be  divided  into  two  parts, 
by  a  plane  passing  through  the  internal  point  perpendicular 
to  the  line  joining  it  with  the  centre,  we  see  that  the  at- 
traction of  the  larger  part  of  the  shell  on  the  two  points  will 
be  ultimately  the  same,  while  the  component  attractions  of 
the  smaller  portions,  in  the  direction  above  considered,  will 
be  equal  in  magnitude  and  opposite  in  direction.     Hence  the 

Y 

difference  between  these  components,  viz.  (2i  +  V) — -,  will  be 

c 

equal   to   twice   the   component   attraction   of  the   smaller 

portion  in  the  direction  of  the  line  joining  the  two  points. 

But  if  Pj  be  the  density  of  the  shell,  6c  its  thickness,  this 

component  attraction  is  'lirp^hc. 

Y 

Hence  (2i  + 1)  —  =  4!^rp^Bc, 
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And,  if  flTj  be  the  corresponding  surface  density, 

2i+l  „ 
^*  =  "^^-^*- 
It  hence  follows  that  if  the  potential  of  a  spherical  shell, 
of  indefinitely  small  thickness,  he  a  surface  harmonic,  its 
potential  at  any  internal  point  will  he  proportional  to  the 
corresponding  solid  harmonic  of  positive  degree,  and  its  po- 
tential at  any  external  point  will  he  proportional  to  the 
corresponding  solid  harmonic  of  negative  degree. 

That  is,  the  proposition  proved  for  zonal  harmonics  in 
Chap.  III.  Art.  6,  is  now  extended  to  spherical  harmonics  in 
general. 

4.  The  spherical  harmonic  of  the  degree  i  will  involve 
2i  +  1  arbitrary  constants. 

For  the  solid  spherical  harmonic,  r^Y.,  being  a  rational 
integral  function  of  x,  y,  z  of  the  i**"  degree,  will  consist  of 

2 terms.      Now  the    expression   V^F,   being  a 

rational  integral  function  of  x,  y,  z  of  the  degree  i  -  2,  will 

f/y  111 

consist  of  -^^ — —-  terms ;  and  the  condition  that  it  must  be 
=  0  for  all  values  of  x,  y,  z,  will  give  rise  to  - — ~-  relations 

among  the  ^— — ^^ L  coefficients  of  these  terms,  leaving 

(i  +  l)(z  +  2)      {i-V)i  .  .    T 

2 9 ,  or  2i  + 1,  mdependent  coefficients. 

5.  We  proceed  to  shew  how  the  spherical  harmonic  of  the 
degree  i  may  be  arranged  in  a  series  of  terms,  each  of  which 
may  be  deduced  by  differentiation  from  the  Zonal  Harmonic 
symmetrical  about  the  axis  of  z.  The  solid  zonal  harmonic, 
which,  in  accordance  with  the  notation  already  employed,  is 
represented  by  r*Pj  {yi),  is  a  function  of  z  and  r  of  the  degree  i, 

satisfying  the  equation  V'F=  0,  or  ^+  ^+  ^k=  0. 

ax        dy        dz 

Now,  if  we  denote  this  expression  by  P^  (z),  we  see  that 
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since  it  is  a  function  of  z  and  r,  it  is  a  function  of  the  dis- 
tance (s)  from  a  certain  plane  passing  through  the  origin,  and 
of  the  distance  (r)  from  the  origin.  Further,  if  we  write  for  z 
the  distance  from  any  other  plane  passing  through  the  origin, 

dW    d^V     d^V 
leaving  r  unaltered,  the  equation  ^-^  +  ^fT'^'^'i  —^  ^^^^ 

continue  to  be  satisfied. 

Now  z  +  a(x  +  -J—lii/),  a  being  any  quantity  whatever, 
represents  the  distance  from  a  certain  plane  passing  through 
the  origin,  since  in  this  expression,  the  sum  of  the  squares 
of  the  coefficients  of  z,  x,  y  is  equal  to  unity.  Hence 
P,  {3  + a  («  +  >/— 1^)}    is   a  solid    zonal  harmonic    of    the 

degree  %  its  axis  being  the  imaginary  line  -  =  —  =  z. 
Therefore  the  equation 

rfx*  ■*"  df  ^  dz^        ' 

is  satisfied  by  V=Pi  [z-\-a  {x  +  'J—ly)],  that  is,  expanding 
by  Taylor's  Theorem,  it  is  satisfied  by 

a*(x  +  ^^^yydT,(z) 
■^        1.2.. .i  dz'     ' 

for  all  values  of  a. 

Hence,  since  the  equation  in  V  is  linear,  it  follows  that 
it  is  satisfied  by  each  term  separately,  or  that,  besides  Pj  (z) 
itself,  each  of  the  i  expressions, 

(.+VrT,)'^^.(.+VrT,).<?:gf),  ...(.W=r,).^). 

satisfies  the  equation  F=  0. 

By  similar  reasoning  we  may  shew  that  each  of  the  ^  ex- 
pressions, 

(Wri,)^),(.W=I,)'™,...(.-VrT,).'?^). 

satisfies  the  same  equation. 
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Now  eacli  of  the  2i  solutions,  thus  obtained,  is  imaginary. 
But  the  sum  of  any  two  or  more  of  them,  or  the  result 
obtained  by  multiplying  any  two  or  more  by  any  arbitrary 
quantities,  and  adding  the  results  together,  will  also  be  a 
solution  of  the  equation.  Hence,  adding  each  term  of  the 
first  series  to  the  corresponding  term  of  the  second,  we  ob- 
tain a  series  of  i  real  solutions  of  the  equation.  Another 
such  series  may  be  obtained  by  subtracting  each  term  of  the 
second  series  from  the  corresponding  term  of  the  first,  and 
dividing  by  V— 1.  We  have  thus  obtained  (including  the 
original  term  Pi{z))  a  series  of  2t  +  l  independent  solutions 
of  the  given  equation,  'which  will  be  the  %i  +  1  independent 
solid  harmonics  of  the  degree  i. 

6.  We  may  deduce  the  surface  harmonics  from  these  by 
writing  r  sin  6  cos  ^  for  x,  r  sin  6  sin  ^  for  y,  r  cos  Q  for  0, 
and  dividing  by  r*.     Then,  putting  cop  ^  =  /*,  and  observing 

that  P,  [z)  =  rT,  {il),  ^4^  =  ^*  ^-^^ ...  we  obtain  the  fol- 

lowing  series  of  2z  + 1  solutions  : 

cos</>sin^^|^\  cos2^sin^^'^'^^^^,  ...  cos^>sin*6>^4^^ 
sm  <^  sm  6      '       ,  sm  2^  sm^  6  — ,  ^  -  ,  . . .  sm  i0  sm'  d      ,'  ^     . 


Expressions  of  the  form 

j>  sin^^ 


(7  cos  0-6  sm^^  — y^^^  , 


or  /S'sino-^sin^^ 


dfi" 


or  their  equivalents, 

^  ^       ^  ^      dfji'' 
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(G  and  8  denoting  any  quantities  independent  of  0  and  <j>) 
are  called  Tesseral  Surface  Harmonics  of  the  degree  i  and 
order  a:  The  particular  forms  assumed  by  them  when 
a-  =  z  are  called  Sectorial  Surface  Harmonics  of  the  degree  i. 

It  will  be  observed  that,  since  — f-f-^  is  a  numerical  constant, 

a/x. 

Sectorial  Harmonics  only  involve  6  in  the  form 

The  product  obtained  by  multiplying  a  Tesseral  or 
Sectorial  Surface  Harmonic  of  the  degree  i  by  r'  (that  is, 
the  expression  directly  obtained  in  Art.  5)  is  called  a  Tesseral 
or  Sectorial  Solid  Harmonic  of  the  degree  i. 

7.    We  shall  denote  the  factor  of  a  Tesseral  or  Sectorial 

Harmonic  which  is  a  function  of  6,  that  is  sin"^^ — r^-^ ,  or 

dfjb'^ 

-  d^P  (a) 
{l—fj?y      1  S^  ,  by  the  symbol  Ty^,  or,  when  it  is  necessary 

to  particularize  the  quantity  of  which  it  is  a  function,  by 
T/«^)(yii)or  J,(-)(cos^). 

It  will  be  convenient,  for  the  purpose  of  comparison  with 
the  forms  of  Tesseral  Harmonics  given  in  the  Mecanique 
Celeste,  and  elsewhere,  to  obtain  Tl-'')  in  a  completely  de- 
veloped form. 

JNow,  smce  P^Ox)  =   ,  i    t>   o — •  — ^"i — ->  we  see  that 
dfi'    ~2M.2.3...t      dfx,'-^'' 

-      ^      ^'^'  L«   *  ^«i-2  4. « (*  -  ^^)  u?^-*    I 
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_*(2^_2)(2^-S)...(*-o■-l)/t^-^-2 

J.  •  ^ 


^-o--l)(^•-o^-2)  (^' 
2.4.(2i-l)(2^-3) 


i-<7--2 


^  (,•_^)(^-^-l)(^-o■-2)(^-<7-3)  ) 


And  therefore 

■^'  2M.2.3...*      ^^    ^^    1^  2(24-1)     ^ 

(^-^)  (^-o-l)(^-(T-2)(^-g-3)  _ 

"^  2.4(2i-l)(2i-3)  ^ 

The  form  given  by  Laplace  for  a  Tesseral  Surface  Har- 
monic of  the  degree  t  and  order  cr  is  (see  Mecanique  Celeste,. 
Liv.  3,  Chap.  2,  pp.  40—47) 

^  (1  -  f.^)^  |/.-<^ -  ^'"^^(ll^)"^^  H''---'  +  -]  ^o^<^<l>' 

A  being  a  quantity  independent  of  6  and  </>.  The  factor  of 
this,  involving  yu,,  is.  denoted  by  Thomson  and  Tait  [Natural 
Philosophy,  Vol.  1,  p.  149)  by  the  symbol  ©/""^  Thomson 
and  Tait  also  employ  a  symbol  ^y>,  adopted  by  Maxwell  in 
his  Treatise  on  Electricity  and  Magnetism,  Vol.  1,  p.  164, 
which  is  equal  to 

1     2      cr 
or  2" tJ^--Z ^w 

(^+a-)(^+c^-l)...(^-o■  +  l)     '    * 
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Heine  represents  the  expression 

■   (^-o^)(^-o— l)(^-a^-2)(^•-o^-3)     ._^_^_ 
2.4.(2t-l)(2«-3}  ^ 

9 

or  (-1)*  ©/'),  by  the  symbol  PJ{ii),  and  calls  these  expres- 
sions by  the  name  Zugeordnete  Functionen  Erster  Art  [Hand- 
buck  der  Kugelfunctionen,  pp.  117,  118)  which  Todhunter 
translates  by  the  term  "Associated  Functions  of  the  First 
Kind,"  which  we  shall  adopt, 

Heine  also  represents  the  series 
(«  — cr)  (/— 0-— 1)     .       „ 

'^  2(2i-l)  '^ 

"^  2.4(2i;-l)(2t-3)  ^ 

by  the  symbol  ^^(/t),  (p.  117). 

The  several  expressions,  Ti''\  Bl'\  ^(<'),  PI,  ^^,  are  con- 
nected together  as  follows : 

T^^  =  0  ('' 


-.^H=(-i)^p^=(i-/.r^!,. 


2^-'^i(i-l)...(q-+l) 
■~(t  +  <r  +  l)  {i  +  <r  +  2)...2i 

8.     It  has  been  already  remarked  that  the  roots  of  the 

equation  P^  =  0  are  all  real.    It  follows  also  that  those  of  the 

dP.  d^P 

equations  -f^  —  ^i  -r-i  =0...  are  real  also.     Hence  we  may 

arrive  at  the  following  conclusions,  concerning  the  curves, 
traced  on  a  sphere,  which  result  from  our  jiutting  any  one 
of  these  series  of  spherical  harmonics  =  0. 

By  putting  a  zonal  harmonic =0,  we  obtain  t  small  circles, 
whose  planes  are  parallel  to  one  another,  perpendicular  to 
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the  axis  of  the  zonal  harmonic,  and  symmetrically  situated 
with  respect  to  the  diametral  plane,  perpendicular  to  this 
axis.  If  i  be  an  odd  number  this  diametral  plane  itself 
becomes  one  of  the  series. 

By  putting  the  tesseral  harmonic  of  the  order  o-=0,  we 
obtain  i  —  <T  small  circles,  situated  as  before,  and  o-  great 
circles,  determined  by  the  equation  cos  a^  =  0,  or  sin  <r^  =  0, 
as  the  case  may  be,  their  planes  all  intersecting  in  the  axis 
of  the  system  of  harmonics,  the  angle  between  the  planes  of 

any  two  consecutive  great  circles  being  -  . 

By  putting  the  sectorial  harmonic  =  0,  we  obtain  i 
great  circles,  whose  planes  all  intersect  in  the  axis  of  the 
system,  the  angle  between  any  two  consecutive  planes  being 

IT 


9.  The  tesseral  harmonic  may  be  regarded  from  another 
point  of  view.  Suppose  it  is  required  to  determine  a  solid 
harmonic  of  the  degree  i,  and  of  the  form  Yf^  such  that  Yj 
shall  be  the  product  of  a  function  of  ix,  and  of  a  function  of  ^, 
which  functions  we  will  denote  by  the  S3mabols  Jf^,  <l>^,  respec- 
tively.    The  differential  equation,  to  which  this  will  lead,  is 


^•(^•+l)^^^^^  +  |{(l-/^^)fj^^  +  T^^^-^^  =  «- 


Now  this  will  be  satisfied,  if  we  make  M.  and  ^^  satisfy 
the  following  two  equations : 

d  (,         2\<^K] 


dfji[  ■^  d/x )      1—  fJ'^ 

d^ 

df 

The  latter  equation  gives 

^i  =  (7  cos  o-</)  +  C  sin  a(f). 

And,  taking  a  as  an  integer,  positive  or  negative,  the 
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former  is  satisfied  by  3f,=  Ty\  i.e.  {I- iJi^fl^Y\l- (ij, 
as  we  proceed  to  prove. 

We  know  that 

Differentiate  a  times,  and  we  get 


{(l-^')f(  +  .-(«  +  l)g  =  0; 


whence,  by  Leibnitz's  Theorem, 


or 


and,  multiplying  by  (1  —  fj^f, 

+  (*-<.)  (t  +  cr+l)(l-;.f^'  =  0...(l). 

Now,  putting      (l-/.^/p  =  r,n 
we  get 
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,<r+a 


i)^(i-/)'^i*-{0 


dfjb 

=(i-/^'r5^'-2(-+i)Mi-^')' -rf^,.x 

(  SI  --0  (7<^P 


./.' 


Hence  the  equation  above  given  for  Jij  is  satisfied  by 
M.  =  T/"'^  and  the  equation  in  I^,  is  satisfied  by 

Y,  =  Cr/'^)  COS  o-(/)  +  0'  r/<^>  sin  o-<^. 

10.  In  Chap.  II.  Art.  10  we  have  established  the  fundamental 
property  of  Zonal  Harmonics,  that  if  i  and  m  be  two  unequal 

positive  integers,   I     P.P^^d/j,  =  0.     This  is  a  particular  case 

of  the  general  theorem  that  if  Y^,  Y,„  be  two  surface  har- 
monics of  the  degrees  ^  and  m  respectively, 

r  rY,Yjf.d4>=^o. 

J  -l.'O 
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For,  let  Vi,  V^  be  the  corresponding  solid  hannonics,  so 
that  V,=r'Yi,  V^  =  r"'Y^.  Then,  by  the  fundamental  pro- 
perty of  potential  functions,  we  have  at  every  point  at  which 
no  attracting  matter  is  situated, 

d^^  dy''^  dz'       *      da?'^  df  ^  dz"  ~   ' 
and  therefore 

'^'[dcc'  '^  df  "^  dz'J      ^'^{dx''^  df'^  dz'J~    ' 

or,  in  accordance  with  our  notation,  F^y''  F^  —  F,„v'  V^  =  0. 

Now,  integrate  this  expression  throughout  the  whole 
space  comprised  within  a  sphere  whose  centre  is  the  origin 
and  radius  a,  a  being  so  chosen  that  this  sphere  contains  no 
attracting  matter.     We  then  have 

jjj{V,v"V^-V^^^'V,)  dxdydz  =  0. 

But  also,  when  the  integration  extends  over  all  space 
comprised  within  any  closed  surface,  we  have 

dS  denoting  an  element  of  the  bounding  surface,  and  -r- 
°  °  '  dn 

differentiation  in  the  direction  of  the  normal  at  any  point. 

Now,  in  the  present  case,  the  bounding  surface  being  a 
sphere  of  radius  a,  and  F],  V^  homogeneous  functions  of  the 
degrees  t,  m,  respectively, 

dS  =  a^df.d4>,    ^  =  ia}-'Y.     ^  =  ma--^Y^, 

and,  the  integration  being  extended  all  over  the  surface  of 
the  sphere,  the  limits  of  fi  are  —  1  and  1,  those  of  ^,  0  and  27r. 
Hence 

r.H.  6 
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whence,  ifmi  —  ibe  not  =  0, 

/I  r^ir 
I     Yidfid(j)  will  be  investigated  here- 
after. 

11.  We  may  hence  prove  that  if  a  function  of  fi  and  <j) 
can  be  developed  in  a  series  of  surface  harmonics,  such  de- 
velopment is  possible  in  only  one  way. 

For  suppose,  if  possible,  that  there  are  two  such  develop- 
ments, so  that 

F{i^,j>)=Y,+  Y,  +  . ..  +  ¥,+  ,.. 

and  also 

Then  subtracting,  we  have 
0=  r,-  F,'+  Y,-  r/+  ...  +  Y,-  Y:  +  ...  identically. 

Now,  each  of  the  expressions  Y^  —  Y/,  I^^  —  Y"/...  F^—  Y"/ 
being  the  difference  of  two  surface  harmonics  of  the  degree 
0,  1,  ...i...  is  itself  a  surface  harmonic  of  the  degree 
0,  1,  ...i....  Denote  these  expressions  for  shortness  by 
Z^,Z^...Zi...  so  that 

0  =  i^o  +  ^  + ...  +  ^j  + ...  identically. 

Then,  multiplying  by  Z^  and  integrating  all  over  the 
surface  of  the  sphere,  we  have 


0=f     p^,V/^<^. 


That  is,  the  sum  of  an  infinite  number  of  essentially 
positive  quantities  is  =  0.  This  can  only  take  place  when 
each  of  the  quantities  is  separately  =  0.  Hence  Z^  is  identi- 
cally =  0,  or  Y/  =  Yf,  and  therefore  the  two  developments 
are  identical. 

"We  have  not  assumed  here  that  such  a  development  is 
always  possible.     That  it  is  so,  will  be  shewn  hereafter. 
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12.  By  referring  to  the  expression  for  a  surface  har- 
monic given  in  Art.  4,  we  see  that  each  of  the  Tesseral  and 
Sectorial  Harmonics  involves  (1  —  fi*)\  or  some  power  of 
(1  _  fj,^^h^  as  a  factor,  and  therefore  is  equal  to  0  when  /x  =  + 1. 
From  this  it  follows  that  when  fi  =  ±l,  the  value  of  the 
Surface  Harmonic  is  independent  of  ^,  or  that  if  Y  {fi,  ^)  repre- 
sent a  general  surface  harmonic,  Y  (±  1,  ^)  is  independent  of 
^,  and  may  therefore  be  written  as  F  (+  1).  Or  Y  (1)  is  the 
value  of  Y{fJ^,  ^)  at  the  pole  of  the  zonal  harmonic  P<  (/a), 
F(— 1)  at  the  other  extremity  of  the  axis  of  P,  (jx). 

We  may  now  prove  that 


r2ir 


Y,d<j>  =  i'wY,{l)F,M. 


For,  recurring  to  the  fundamental  equation, 

Now,  if  we  integrate  this  equation  with  respect  to  <p, 
between  the  limits  0  and  27r,  we  see  that,  since 


/ 


d'Y         dY, 


and  the  value  of  F<  only  involves  (j)  under  the  form  of  cosines 

or  sines  of  <b  and  its  multiples,  and  therefore  the  values  of 

dY 

-jT*  are  the  same  at  both  limits,  it  follows  that 

Hence 


/: 


Hence    I     F/?^  is  a  function  of  /m  which  satisfies  the 
Jo 
fundamental  equation  for  a  zonal  harmonic,  and  we  therefore 
have 

6—2 
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rSTT 


f 

Jo 


Y,dcf>=CPM> 

0 

G  being  a  constant,  as  yet  unknown. 

To  determine  C,  put  fi=l,  then  by  the  remark  just  made, 
Yf  becomes  ^"^(1),  and  is  independent  of  <}>.    Hence,  when 

^=1,  pY;tZ<^  =  27rr.(l).     Also  P, (/Li)  =  1.     We  have  there- 
fore      "  27rr,(l)  =  C: 

Jo 
It  follows  from  this  that 

13.     We  may  now  enquire  what  will  be  the  value  of 

jj^   Y.Z.d/Mdcj,, 

Yf,  Z.  being  two  general  surface  harmonics  of  the  degree  i. 
Suppose  each  to  be  arranged  in  a  series  consisting  of  the 
zonal  harmonic  P,  whose  axis  is  the  axis  of  z,  and  the  system 
of  tesseral  and  sectorial  harmonics  deduced  from  it.  Let  us 
represent  them  as  follows : 
Y,=  AP, 

+  C,Tp cos  (^  +  (7,T,(2)  cos  2^  +  ...  +  a^i^'^  cos o-<^  +  ... 

+  C;T©cos«<^ 

+  S,TPsm(f>  +  S,Tp  sin  2(/)  +  ...  +  S^T,^''^  sm<T<f>  +  ... 

^.=  aP, 

+  Cj  r/i)  cos  (ji  +  c^  7^(2)  COS  2</)  +  . . .  f  c.,^/"^)  cos  o-<^  +  . . . 

c^T^^"^  COS  i^ 
+  s, T.(i)  sin  <!>  +  s^Tp  sin  2c^  + . . .  +  s'^Ti''^  sin  o-</)  +  . . . 

Hence  the  product  Y^Z.  will  consist  of  a  series  of  terms, 
in  which  ^  will  enter  under  the  form  cos  a-cp  cos  a<J3,  or 
cos  a^  sin  c'<f).     This  expression  when  integrated  between 
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the  limits  0  and  27r  vanishes  in  all  eases,  except  when 
(7=  a-  and  the  expression  consequently  becomes  equal  to 
cos"  o-^,  or  sin^^cr^.  In  these  cases  we  know  that,  a  being  any 
positive  integer, 

/'2ir  rErr 

I     cos^cr^cZ^sr       sbo?  o-^d<f>  =  IT. 
Jo  Jo 

Hence  the  question  is  reduced  to  the  determination  of  the 
value  of 


j[m*df^. 


Now  TW  =  (1-/^T"^ 


'     d-^p^^^-^^' 


But,  by  the  theorem  of  Rodrigues,  proved  in  Chap.  II. 
Art.  8,  we  know  that 


d^i+-    ~     ^  lizs  ^^'~'    ' 

Hence  ^^/"^^  may  also  be  expressed  under  the  form 
whence  it  follows  that 

m<r)Y-(.iy  (       1       V\^±^d^^'^{^.'-lydi-^{^.^-ly 

K^i    }  -K    ^)    1^2M.2.3...i7  \i-a-       d/M'+'^  dfi^-'^      ' 

Now,  putting  {fjb^  —  iy  =  M  for  the  moment,  and   inte- 
grating by  parts, 

Cd^-^^Md^-'M  ,    _d^^''-'^ M d'-'' M 
}  iLfJiP^  "djiF^   ^  ~  "diiF^^  d^jf-"' 
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The  factor  —r-- —  vanishes  at  both  limits,  hence 

by  a  repetition  of  the  same  process. 

And  by  repeating  this  process  <r  times,  we  see  that 

2 


Hence  /     (T/-^))' cZya  =  ^ 


1 1  +  q-      2 

^ "^27+1' 

and  therefore 

{T,(''^co&a<}>ydfidc{>=  {T('^'>  sin  a(j>ydfid(f> 

J  -I  J  0  J -iJ  Q 

_['i  +  o-     27r 
~[i-o-  2i  +  l* 

It  will  be  observed  that  this  result  does  not  hold  when 
cr  =  0,  in  which  case  we  have 


Hence  I     rY,Z,d/xd<l> 


47r      , 


2i  + 1 
*  In  this  case  J   '^cos-<r^cZ0  =  f  ^ sia?(r<pd^=2w. 
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"      Stt     fh'  +  l  l»+2 

+  27Ti  1[^  (^'">  +  ^'*')  +  \m  ^^'''+  ^^^  +  •••• 

.  14.  We  have  hitherto  considered  the  Zonal  Harmonic 
under  its  simplest  form,  that  of  a  "  Legendre's  Coefficient "  in 
which  the  axis  of  z,  i.  e.  the  line  from  which  0  is  measured,  is 
the  axis  of  the  system.  We  shall  now  proceed  to  consider  it 
under  the  more  general  form  of  a  "Laplace's  Coefficient," 
in  which  the  axis  of  the  system  of  zonal  harmonics  is  in  any 
position  whatever,  and  shall  shew  how  this  general  form  may 
be  expressed  in  terms  of  F^  {fj,)  and  of  the  system  of  Tesseral 
and  Sectorial  Harmonics  deduced  from  it. 

Suppose  that  0',  ^'  are  the  angular  co-ordinates  of  the 
axis  of  the  Zonal  Harmonic,  i.e.  that  the  angle  between  this 
axis  and  the  axis  of  z  is  ff,  and  that  the  plane  containing 
these  two  axes  is  inclined  to  a  fixed  plane  through  the  axis 
of  z  which  we  may  consider  as  that  of  zx,  at  the  angle  0'. 
In  accordance  with  the  notation  already  employed,  we  shall 
represent  cos  0'  by  /x'. 

The  rectangular  equations  of  the  axis  of  this  system 
will  be 

^        ^        y       ^^    -g 

sin  &  cos  <^'      sin  6'  sin  ^'     cos  ff  ' 

Hence  the  Solid  Zonal  Harmonic  of  which  this  is  the  axis 
is  deduced  from  the  ordinary  form  of  the  solid  zonal  har- 
monic expressed  as  a  function  of  z  and  r  by  writiug,  in  place 
of  z,  X  sin  0'  cos  <f>'  +  y  sin  0'  sin  <f>'  +  z  cos  0'. 

To  deduce  the  Surface  Zonal  Harmonic,  transform  the  solid 
zonal  harmonic  to  polar  co-ordinates,  by  writing  rsin^cos^ 
for  X,  r  sin  0  sin  <f)  for  y,  r  cos  0  for  z,  and  divide  by  r*. 

The  transformation  from  the  special  to  the  general 
form  of  surface  zonal  harmonic  may  be  at  once  effected, 
by  substituting  for  fi,  or  cos  0,  cos^cos^'+sin^sin^cos(^— ^'). 

Now,  in  order  to  develope 

P,  {cos  ^  cos  ^  +  sin  ^  sin  ff  cos  {<}>  —  (f>')] 
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in  the  manner  already  pointed  out,  assume 

P,  {cos  9  cos  6'  +  sin  6  sin  6'  cos  (<^  -</>') } 

=  AP,  (ji)  +  ((7(1)  cos  </)  +  >§(!)  sin  <^)  TP 

+  (C(2)  cos  2(^  +  /S(2)  sin  2c/))  2!(2)  + ... 

+  ( C^'^)  cos  a<j>  +  /S^-^)  sin  <t<}>)  T W  +  . . . 

+  ( (7<^)  cos  i<^  +  /S^^)  sin  i(f>)  Tp, 

the  letters  A,  ...  (7^^\  S^'^K..  denoting  functions  of  fx!  and 
^',  to  be  determined. 

To  determine  C^''\  multiply  both  sides  of  this  equation 
hycoso-cpT}"^  and  integrate  all  over  the  surface  of  the  sphere, 
i.e.  between  the  limits  —  1  and  1  of  fi,  and  0  and  Stt  of  (p. 
We  then  get 

/•I       rZTr 

P,  {cos  e cos 0'  +  sin  6 sin d' cos {<J3 - (ji)] cos o-<^ 2^(<^>  t7/icZ</) 

=  CWr    r  {cob acf>T(-^ydfid({> 
h'  +  q-    27r 


I  -  a-  2^  +  1 


Of''). 


It  remains  to  find  the  value  of  the  left-hand  member  of 
this  equation. 

Now  cos  cr(j)  Ti'^'>  is  a  surface  harmonic  of  the  degree  i,  and 
therefore  a  function  of  the  kind  denoted  by  y^  in  Art.  12. 

And  we  have  shewn,  in  that  Article,  that 


a 


2t  Ai-TT 


that  is,  that  if  any  surface  harmonic  of  the  degree  i  he  multi- 
plied hy  the  zonal  harmonic  of  the  same  degree,  and  the  product 
integrated  all  over  the  surface  of  the  sphere,  the  integral  is 

equal  to  ^. — r-  into   the  value  which  the  surface  harmonic 
2i  -f  1  -^ 

assumes  at  the  pole  of  the  zonal  harmonic. 
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Hence 

•1      r2tr 


ri    r2tr 

I     J    P,  {cos  e  cos  e'  +  sin  d  sin  ff  cos  (</>  -  </>')]  1\  {fi,  <^)  (?/t  J(^ 


+ 
and  therefore 


=  2^r.(^',^0, 


ri   rzir 

I     I     P,  [cos  ^  cos  ^  +  sin  6  sin  6>'  cos  (<^  -  </>')}  cos  <T(f>T^''^d//.d(f) 

=  2-^  cos  .f  21'./.  . 

Hence 

or  0(')  =  2  H^^  cos  a<f>  W  (>')• 

Similarly  ,S(<^)  =  2  Li~5  sin  o-^'  i;^'^)  (/). 

And  to  determine  A,  we  have 

/I    f2ir 
I    P,  {cos  e  cos  ^'  +  sin  ^  sin  ^  cos  (<^  -  <^') }  P,  (/^)  dfidcfi 

=  Af_J^{PAH')YdfMdcf>', 

or    ^  =  P,(/). 
Hence,  P^  [cos  6  cos  ^'  +  sin  0  sin  ^  cos  (^  —  <f)')] 

=  P,  (/.')  P,  (/.)  +  2  ^cos  (0  -  </,')  7^(1)  (/)  T/^>  (/.) 

+  2  t^i  cos  2  ((^  -  <^')  T(2)  (/.')  P/^)  Ga)  +  ... 

t  "t"  ^ 
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Li± 


cr 


+  2  i  COS  z  (<^  -  f )  ^/n/)  ^>M/^). 

15.  We  have  already  seen  (Chap.  II.  Art.  20)  how  any 
rational  integral  function  of  /m  can  be  expressed  by  a  finite 
series  of  zonal  harmonics.  "We  shall  now  shew  how  any 
rational  integral  function  of  cos  d,  sin  6  cos  ^,  sin  6  sin  ^, 
can  be  expressed  by  a  finite  series  of  zonal,  tesseral,  and 
sectorial  harmonics. 

For  any  power  of  cos  </>  or  sin  ^,  or  any  product  of  such 
powers,  may  be  expressed  as  the  sum  of  a  series  of  terms  of 
the  form  cos  a^,  or  sin  o-^,  the  greatest  value  of  <t  being  the 
sum  of  the  indices  of  cos  (f>  and  sin  (jj,  and  the  other  values 
diminishing  by  2  in  each  successive  term.  Hence  any 
rational  integral  function  of  cos  6,  sin  6  cos  (}>,  sin  0  sin  <f),  will 
consist  of  a  series  of  terms  of  the  form 

cos"*  6  sin"  0  cos  a^  or  cos'"  0  sin"  0  sin  a^, 

where  n  is  not  less  than  a. 

If  n  be  greater  than  cr,  w  —  a-  must  be  an  even  integer.  Let 
n  —  or  =  2s,  then  writing  sin"^  under  the  form  (1  —  cos^^)*  sin"^^, 
we  reduce  cos™  0  sin"  0  cos  cr^  to  the  sum  of  a  series  of  terms 
of  the  form  cos^  0  sin"'  0  cos  (7(f),  or,  writing  cos  0  =  fi,  oi  the 

form  fi"  (1  —  /i^)  *  cos  o-^. 

Similarly    cos"*  0  sin"  ^  sin  cr^  is  reduced  to  a  series  of 

terms  of  the  form  fi^  (1  —  fi^y  sin  (r<p. 

Now        LLP  =  7 r-7 r^ > :r^  -7—  /*^      > 

and  fi^'^'^  can  be  developed  in  a  series  of  terms  of  the  form 
of  multiples  of  Fp+<T}  Pp+a-2  •••  •     (Chap.  II.  Art.  17.) 

Hence  fiP  can  be  expressed  in  a  series  of  the  form 

Jo- 
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Aq,  A^  representing  known  numerical  constants,  and  therefore 

fjP  i}.  —  jj^y  assumes  the  form 

(-^0  Tp+v  +  A^  Tp+^-2  +  ...), 
consequently  multiplying  these  series  by  cos  o-^  or  sin  o-^,  we 
obtain  the  developments  of 

/^  (1  —  /i')  ^  COS  o-^  and  ^  (1  —  /u,')  "^  sin  a<^ 
in  series  of  tesseral  harmonics. 

16.    We  will  give  two  illustrations  of  this  transformation. 

First,  suppose  it  is  required  to  express  cos'  6  sin*^  sin  ^  cos^ 
in  a  series  of  Spherical  Harmonics. 

Here  we  have    sin  ^  cos  ^  =  ^  sin  20. 

Hence  cos*  6  sin'  B  sin  ^  cos  ^  =  ^  cos'  0  sin'  B  sin  20. 

Comparing  this  with  cos"*  ^  sin"  B  sin  <r0,  we  see  that  n  is 
not  greater  than  <r. 

Hence  cos'  B  sin'  ^  sin  0  cos  0  =  -  /x'  (1  —  /i,')  sin  20. 

and  /i*  =  g^P,  +  ^P,  +  -P„, 

.  J.  /_8^^     4^\ 

••  '^      12V35  dyi^'^l  d/^V 
^    2   cf  P       1  <f  P, 
105  dii:      21  dyu," ' 
.*.  cos'  B  sin'  ^  sin  0  cos  0 
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Next,  let  it  be  required  to  trac^form  cos^  6  sin^  d  sin  <f  cos^  (f> 
into  a  series  of  Spherical  Harmonics. 

1  1 

Here  sin  <f)  cos^  <^  =  ^  sin  2</)  cos  <^  =  ^  (sin  3<^  +  sm  <f)). 

Now  cos'  e  sin'  6  sin  d(}>  =  //-'(!-  /i')^  sin  3<^ 

Also  cos'  ^  sin'  ^  sin  (^  =  /u-'  (1  -  /a')  (1  -  ft')*  sin  (^ 
=  (/.=^-/.^)(l-)a^)^sin«^ 


=  (i|^^-||^0^1-'^^)^-^- 


Also  (Chap.  IL  Art.  17) 


^^=3-5^^+^^^n^o' 


^e_16  p  .24  10  1 

^  ~23l    «'^77    ^^Sl^^'^T    "• 
Hence  cos'  6  sin'  ^  sin  Scji 

1/16  c^'P,     24J'PA  ,T        2^^   •    o^ 

=  120(231^ +r7 5?)  (^-^^'^^^^^ 

Andcos'^sm  ^sm<^  =  - (^^^3-^ +^^.  +  ^3^ 

2   dP.       1  JPA  ,,  srsi    •      . 

So  dfj,      7  dfj,J^       ^^         ^ 

[(i9'3dfM      SH5  dfx,      (iSdfi)^       ^;   sm^ 

=  -(  ^^  T<y>  -  -^  T  (1)  -  —  Tj^A  sin  6  ; 
^693    -^        385^        63    W        ^' 

.-.  cos'^  sin'^  sinc^  ^^^'^  =  {69^  ^«^'^  +  1^  r;3>|sin3^ 

_  I  A.  r  (1)  -  -!-  T  (1) -  r  ^^'1  sinrf). 

[693^^       770-^*       63    M 
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17.  The  process  above  investigated  is  probably  the  most 
convenient  one  when  the  object  is  to  transform  any  finite 
algebraical  function  of  cos^,  sin^iCos<^,  and  sin  ^  sin  0,  into 
a  series  of  spherical  harmonics.  For  general  forms  of  a 
function  of  //.  and  ^,  however,  this  method  is  inapplicable, 
and  we  proceed  to  investigate  a  process  which  will  apply 
universally,  even  if  the  function  to  be  transformed  be  discon- 
tinuous. 

We  must  first  discuss  the  following  problem. 

To  determine  the  potential  of  a  spherical  shell  whose 
surface  density  is  F(ji,^),  -F' denoting  any  function  whatever 
of  finite  magnitude,  at  an  external  or  internal  point. 

Let  c  be  the  radius  of  the  sphere,  /  the  distance  of  the 
point  from  its  centre,  6',  ^'  its  angular  co-ordinates,  Fthe 
potential.     Then  fjb  being  equal  to  cos  0 


F{fi,  ^)  c^dtid(f> 


0 '  [/'-  2cr  {cos  0  cos  0'  +  sin  0  sin  0'  cos  {<}> -<!>')]  +  c']^  ' 

The  denominator,  when  expanded  in  a  series  of  general 
zonal  harmonics,  or  Laplace's  coefficients,  becomes 

i|l+P,(/.,^)^4-P,(^,^)^V...  +  P,0.,<^)^'  +  ...|, 

for  an  internal  and  an  external  point  respectively,  P^  {fi,  ^) " 
being  written  for 

P,  {cos  0  cos  0'  +  sin  0  sin  0'  cos  (<f>  -  ^')}. 

Hence,  V^  denoting  the  potential  at  an  internal,  V^  at  an 
external,  point, 

^1  =  0  [f_fyiH'>  <!>)  dtid<\>  +  ^  £^  £"P,  (/.,  (^)  F{^l,  j>)  dfMdc}> 
+ ...  +7  f  fjPiiM',  ^)  F{f^>  </>)  dfMdcf>  +  ...|, 
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It  will  be  observed  that  the  expression  P^  (/u-,  ^)  involves 
/A  and  /a'  symmetrically,  and  also  <}>  and  (}>.  Hence  it  satisfies 
the  equation 


^•{(1-"") 


And,  since  fi  and  ^  are  independent  of  fi  and  ^',  this 
differential  equation  will  continue  to  be  satisfied  after  P,  has 
been  multiplied  by  any  function  of  /u,  and  cf),  and  integrated 
with  respect  to  fj,  and  ^.  That  is,  every  expression  of  the 
form 


fl     C2rt 


is  a  Spherical  Surface  Harmonic,  or  "Laplace's  Fuuction" 
with  respect  to  //  and  ^'  of  the  degree  i.  And  the  several 
terms  of  the  developments  of  V^  are  solid  harmonics  of  the 
degree  0,  1,  2...i...  while  those  of  V^  are  the  corresponding 
functions  of  the  degrees  —  1,  —2,  —  3...  —  (i+  1), ...  And 
these  are  the  expressions  for  the  potential  at  a  point  (r,  fi,  (fi) 
of  the  distribution  of  density  F{/jf,  (f>')  at  a  point  (c,  fi,  ^'). 

Now,  the  expressions  for  the  potentials,  both  external 
and  internal,  given  in  the  last  Article,  are  precisely  the  same 
as  those  for  the  distribution  of  matter  whose  surface  density  is 


~|[' J'V(>,  cj>)  dixd<\>  +  S/'JJV.O,  ^)  Fi/J,,  <^)  dfj,dcj,+  ... 


4i7r 

+ 


or,  as  it  may  now  be  better  expressed, 
47r 


\[jy{l^,  <!>)  diMd4> 
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Pj  {cos^  cos^'+  sin^  sin^  cos  (^-«^')  Fiji,  <f>)  dfid<f) 
+  ... 

+  (2i  +  l)l    /  P,{cos^cos^'+sm^sm^'cos(^-<^')l-^(/*'^)^/^^^+---    • 

And,  since  there  is  only  one  distribution  of  density  which 
will  produce  a  given  potential  at  every  point  both  external 
and  internal,  it  follows  that  this  series  must  be  identical 
with  F(jjb,  (f>).  We  have  thus,  therefore,  investigated  the 
development  of  F  {//,',  <j>)  in  a  series  of  spherical  surface 
harmonics*. 

The  only  limitation  on  the  generality  of  the  function 
F(fjo',  0')  is  that  it  should  not  become  infinite  for  any  pair  of 
values  comprised  between  the  limits  —1  and  1  of  /it,  and  0 
and  2'ir  of  <p. 

18.  Ex.  To  express  cos  2</>'  in  a  series  of  spherical  har- 
monics. 

For  this  purpose,  it  is  necessary  to  determine  the  value  of 
(2i  +  1)1    I    Pj  [cos^  cos ^+ sin ^  sin ^'  cos  (^-^')|  cos  2j>dfj,d(f). 

Now  Pi  [cos  9  cos  6'  +  sin  0  sin  ^  cos  ((/>  -  <^') } 

=  P,(cos^)P,(cos^) 

^       2        .   ^dP,(cose)   .    ^dP,  (cos 6')        .^     ^,, 

'^{i-l)i{i  +  l){i  +  2f^''^ 

d'P, (cos 6)   .  a^,6fP,(oos(90        „,,       ,,. 
dM,'       "^^^        df.-     '^os2(<f>-<i>)  +  ... 

Now  I     cos  cr{j>  —  (j)')  cos2<pd<l>  =  0, 

for  all  values  of  a  except  2. 

*  In  connection  with  the  subject  of  this  Article,  see  a  paper  by  Mr  G.  H. 
Darwin  in  the  Messenger  oj  Mathematics  for  March,  1877. 
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/•2ir  ^ 

And         I    cos  2(^  — ^')  cos2(^(i'^  =  7rcos2^. 

Jo 
Also 

j^^^^^^^  =  -2M.2.3...J^  -1)       ^^..       d^. 
And 

Now  when  /-t  =  1, 

^     ^^       di^^'  '    ^  *      '  dfju'-'  "• 

And  when  /a  =  —  1, 

Hence 

/\in'«g<?^  =  ^-3-.2'.1.2.3...qi-(-in 

=  4  or  0,  as  i  is  even  or  odd ; 

.  • .  I  sm'^  0  — ^~ — -  cos  2  (^  -  ^ )  cos  2(p  dfid^ 

=  47r  cos  2(^'  or  0,  as  i  is  even  or  odd ; 

.*.  cos 2^' 

.    .  2 ^,  cZ'P, (cos ^)  _,, 

4  sm''  ^  ^^Si — ^  -Tr  cos  2<f> 


47r  I    1 . 2 .  3 .  4  dfi 

.    n         2         .    .   ,^cfP,(cos0')  „,, 
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* 1 

/    K    7T(2)  q    7^(2)  1Q    7^(2)  \ 

Hence  the  potential  of  a  spherical  shell,  of  radius  c  and 
surface  density  cos  2(f)',  will  be 

/       gr^(2)         ^'2  rp{2)  ^'4  y(«         ^'6  V 

and 

/        y(2)  g3  ^(S)  ^5  /y7(2)  ^7  N 

87rC0s2(/>'      — ^r^-r    -73  +  »    /.    ^    -75+.     ^V    q        ,+   ...., 

^  Vl-2. 3. 4<r       3. 4.0. or        !i.b.7.8r^  J 

at  an  internal  and  external  point  respectively. 

19.  We  will  now  explain  the  application  of  Spherical 
Harmonics  to  the  determination  of  tlie  potential  of  a  homo- 
geneous solid,  nearly  spherical  in  form.  The  following 
investigation  is  taken  from  the  Mecanique  Celeste,  Liv.  in. 
Chap.  II. 

Let  r  be  the  radius  vector  of  such  a  solid,  and  let 
r  =  a->ra{a^Y^-ira^Y^-^ ...  +aiF, +  ...), 
a  being  a  small  quantity,  whose  square  and  higher  powers 
may  be  neglected,  a^,  a.^,...a^...  lines  of  arbitrary  length,  and 
1\,  Fjj,...yj...  surface  harmonics  of  the  order  1,  2,...i...  re- 
spectively. 

The  volume  of  the  solid  will  be  _  Tra". 

o 

For  it  is  equal  to 

I    I  r^drdjxdj) 

=  Q7ra^  since         I      Y^d^d(\)=^Q, 

o  J  -1  J  0 

for  all  values  of  i. 

F.  H.  7 
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Again,  if  the  centre  of  gravity  of  the  solid  be  taken  as 
origin,  o^  =  0. 

For  if  z  be  the  distance  of  the  centre  of  gravity  from  the 
plane  of  xi/, 


-1.'0 

'27r 


=  \  i    ('''a*  +  4!a'x  {a^J\  +  aJ\+,..+  c,Y,+  ...)dfid(l> 

**  J  -1  J  0 

=  4a'j.(Xj  fjil\d/.id(}>. 

J  -1  .'  0 
Similarly 

,,  izd'x  =  4a^  a  . «!  (1  —  fi^)^  cos  <^  J\  dfid(p, 

o  J  -iJ  0 

-  Tra^y  =  4r/ (x.  a^  I  (1  —  /x^)i sin  (f)  J\  dfid^. 

^  '  J  -1 J  0 

Now  Y^  is  an  expression  of  the  form 

Afjb  +  B{1—  fi^)^  cos  0  +  C(l  —  //,^)^  sin  (^, 

a!id  therefore  all  the  expressions  x,  y,  z  cannot  be  equal  to  0, 
unless  a^  =  0. 

"We  may  therefore,  taking  the  centre  of  gravity  as  origin, 
write 

r  =  a+a{a^Y^+  ...  +  a.Yi+  ...), 

as  the  equation  of  the  bounding  surface  of  the  solid. 

Now  this  solid  may  be  considered  as  made  up  of  a  homo- 
geneous S2)here,  radius  a,  and  of  a  shell,  whose  thickness  is 

a(a,F,+  ...  +  a,r.-t-  ...). 

The  potential  of  this  shell,  at  least  at  points  whose  least 
distance  from  it  is  consideral)le  compared  with  its  thickness, 
will  be  the  same  as  that  of  a  shell  whose  thickness  is  aa,  and 
d'jnsity 

r^a  ^ '+■■■  + a  ^ '  +  ■■■■}■ 
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Pq  being  the  density  of  the  solid.  Therefore  the  potential, 
for  any  external  point,  distant  R  from  the  centre,  will  be 

^■^P^M  +  ^W^^  C'^l'^  +  ••••  +^  J^  +  ••••)  • 

The  potential  at  any  internal  point,  distant  R  from  the 
centre,  will  be  made  up  of  the  two  portions 

I  irp.B'^  +  27rp„  {a-^  -  R^)  or  27rp  (a^  - 1^)  . 

for  the  homogeneous  sphere, 

^^       \   o     or  '2i  +  l  a  J 

for  the  shell,  and  will  therefore  be  equal  to 

20.     If  the  solid,  instead  of  being  homogeneous,  be  mide 

up  of  strata  of  different  densities,  the  strata  being  concentric, 

and  similar  to  the  bounding  surface  of  the  solid,  we  may 

c 
deduce  an  expros>ion  for  its  potential  as  follows.     Let  -  r  be 

the  radius  vector  of  any  stratum,  p  its  density,  r  having  the 
same  value  as  in  the  last  Article,  and  p  being  a  function 
of  c  only.  Then,  8c  being  the  mean  thickness  of  the  stratum, 
that  is  the  difference  between  the  values  of  c  for  its  inner 
and  outer  surfaces,  the  potential  of  tlie  stratum  at  an  ex- 
ternal point  will  be 

^irpc'Sc      ,         c'BcfaYc'      a,Yn' 
R  ^     a    \   '0     R'        1    R^ 

To  obtain  the  potential  of  the  wliole  solid  at  an  external 
point  we  must  integrate  this  expression  with  respect  to  r, 
between  the  limits  0  and  a,  remembering  that  p  is  a  func- 
tion of  c. 

7—2 
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Again,  the  potential  of  the  stratum,  above  considered, 
at  an  internal  point  will  be 

-^f--) (^)- 

To  obtain  the  potential  of  the  whole  solid  at  an  internal 
point  we  must  integrate  the  expression  (1)  with  respect  to  c 
between  the  limits  0  and  R,  and  the  expression  (2)  with 
respect  to  c  between  the  limits  R  and  a,  remembering  in 
l)oth  cases  that  />  is  a  function  of  c,  and  add  the  results 
tosrether. 


CHAPTER  V. 

SPHERICAL  HARMONICS  OF  THE  SECOND  KIND.     ' 

1.  We  have  already  seen  (Chap.  ii.  Art.  2)  that  the 
differential  equation  of  which  P^  is  one  solution,  being  of 
the  second  order,  admits  of  another  solution,  viz. 

Now  if  /i  between  the  limits  of  integration  be  equal 
to  +  1,  or  to  any  roots  of  the  equation  P,  =  0  (all  of  which 
roots  lie  between  1  and  —1),  the  expression  under  the 
integral  sign  becomes  infinite  between  the  limits  of  inte- 
gration. We  can  therefore  only  assign  an  intelligible 
meaning  to  this  integral,  by  supposing  /t  to  be  always  be- 
tween 1  and  00  ,  or  between  —  1  and  —  oo  .  We  will  adopt 
the  former  supposition,  and  if  we   then  put   (7=  — 1,   the 

expression  ~pri-rzi~^\  V'^'  p^  (  '_  iv  ^^^  ^^  always  posi- 
tive.    We  may  therefore  define  the  expression 

dfjb 

as  the  zonal  harmonic  of  the  second  kind,  which  we  shall 
denote  by  (2,-,  or  Q.  (/a),  when  it  is  necessary  to  specify  the 
variables  of  which  it  is  a  function. 


r 

J  ft. 


It  will  be  observed  that,  if  fx  be  greater  than  1,  P^  is 
always  positive.  Hence,  on  the  same  supposition,  Qi  is 
always  positive. 

We  see  that  Q^  =  I     —, — r  =  ^  log , 


102        SPHERICAL  HARMONICS   OF  THE   SECOND   KIND, 


1  fM  +  1 

And,  in  a  similar  manner,  the  values  of  Q.^,  Q^,...  may 
be  calculated. 

2.  But  there  is  another  manner  of  arriving  at  these 
functions,  which  will  enable  us  to  express  them,  when  the 
variable  is  greater  than  unity,  in  a  converging  series,  with- 
out the  necessity  of  integration. 

This  we  shall  do  in  the  following  manner. 

Let  U= ,  V  being  not  less,  and  /a  not  greater,  than 

unity. 

Then  ^=__J__  ^^^       1 

dv         {v  —  fiy  '  d/x~  {v  —  fif ' 

..        ,.dU      v-^-l        ..        ,.dU      1-fi' 


du         {l^-fl)"     ^         ^^   dfl~  {v-fif 

Aig^,?)^^  -    ^'-1   f^ 2    \_^l-tiv 

1 

Now,  let be  expanded  in  a  series  of  zonal  harmonics 

V  —  IJ,  '■ 


r=   1 


■P.W,  -P,W...P,(m),  sothat 
by  the  definition  of  P  (/x) . 
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And  also  ,4(a-.)fj  =  ....|,{a-„.)^*^}p,M,... 

And  these  two  expressions  are  equal.  Hence,  equating 
the  coefficients  of  P^  {/j,), 

Hence  <f>i{^)  satisfies  the  same  differential  equation  as  /* 
and  Q,.  But  since  U=  0  when  j/'=  oo  ,  it  follows  that  ^j(i/)  =  () 
when  v  =  <x> .  Hence  (f).{v)  is  some  multiple  of  Qi{v)=AQi{if) 
suppose.     It  remains  to  determine  A. 

Now,  <f>i{v)  may  be  developed  in  a  series  proceeding  by 
ascending  powers  of  -  ,  as  follows. 

V 

We  have      =  - +  ^i  + ...+  ^,  + .... 

V  —  fJb        V        V  V 

and  also        =  <^,{v)  PM+cf>Xp)  PM  +...  +  <^,M  PM  +  ... 

Now,  by  Chap.  II.  Art.  17,  we  see  that,  if  m  be  any 
integer  greater  than  i,  the  coefficient  of  P<  in  /a"'  is 

/o-  .  i\         (m-«+ 2)  (w-t  +  4)...(m-l)  .,  . ,        ,, 

(2i  + 1)  -7 — > . — ^,^-7-  ^—  . — TT — y — --■--  ,        ^,  if  I  be  odd, 
^  ^  (m  +  i  +  1)  (m  +  i-l)...(w  +  4)  (w  +  2) 

and  (2i  4- 1)  7^- — .—m ■    r^ — 7 — ~.^'-^ — -ttii  *  be  even, 

m  —  {  being  always  even. 

Hence,  writing  for  m  successively  i,  i  +  2,  «  +  4,  ...  we  get 
.  ,  ^      .«•     -.N  f  2.4...(i-l)  1 

4.G...0'+1^  1 


+ 


(2t+3j  (2i  +  l)...(i  +  4)  V 


4+3 


6.8...(t  +  3)  1  ,  ..  .,        ,, 

+  TT-^— -x-r^-.     .., — -r. — TT  -;+5  +  ••••r  II  «  be  odd, 
(2i  +  o)(2i  +  3)...(i  +  6)  i/"^' 


104        SPHERICAL   HARMONICS   OF  THE  SECOND  KIND. 
/«.       niN   f  2.  4... I  1 

^^^  =  ^2^+ ^)  Wi)(2i-i)...(i+i)  ];7^ 

4. 6. ..(1  +  2)  J_ 

'^(2i  +  3j(2i  +  l)...(t+8)'i.'^' 

6.8...(t  +  4)  1    ,       ]  .„  ., 

+  yK-- — ^,  ,,-. .     .-,. — n — ^  -71T+  ....^  it  t  be  even. 

Now,  recurring  to  the  equation 


.p,(^)i-(^-^-i)' 

we  see  that,  if  Qi(v)  be  developed  in  a  series  of  ascending 

powers  of  -,  the  first  term  will  be  ~7^^:—~t\ — ^,  where   G 

is  the  coefficient  of  jJ  in  the  development  of  jPi(^) ; 

(i  +  2)(z  +  4)...(2/-l) 
that  IS    6  =  — — s^— 7,-  S .-   V  r  —  if  «  he  odd, 
2.  4.b...(i— 1)  ' 

and  =  (i+ lHi_+  3)  («:+ 3)..,(2,-l)  .^  .  ^^  ^^^^^ 

Hence  the  first  term  in  the  development  of  Q^  {v)  is 
2.4.6...(t-l) 


(•i  +  2)(i  +  4)...(2i-lj  (2i+l) 


if  I  be  odd, 


,  2.4.6...^  .„., 

and         =  ..  ,  _,.  ,. — „, — ,  - — ■-.  ■-,. — z-r  II  *  be  even, 
(i+l)  (i  +  8)...(2i-l)  (2^  +  l) 

Avhich  is  the  same  as  the  first  term  of  the  development  of 

Pi{v),  divided  by -r-^ — =-. 

Hence  A  =  2i+  1,  and  we  have 

3.     The  expression  for  Q.  may  be  thrown  into    a    more 
convenient  form,  by  introducing  into  the  numerator  and  de- 
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nominator  of  the  coefficient  of  each  term,  the  factor  neces- 
sary to  make  the  numerator  the  product  of  i  consecutive 
integers.  We  shall  thus  make  the  denominator  the  product 
of  i  consecutive  odd  integers,  and  may  write 

^,,  1.2.3...t        1        3.4.5...(i  +  2)     1 


1.3.5...(2i+l)  y'^'     3.5.7...[2i  +  S)v' 

5.6.7...(t  +  4)     1 
5.7.9...(2t+5)  v' 


+  ^-.^t;— 7k .  — ^V  7m  +  . .  . 


(2k +1)  {2k  +2)...{i+2k) 1_ 

"^  {2k  +  1)  (2^^  +  3) . . . (27  +  2k'+l)  p'^"^-' ■*"■*" 
whether  i  be  odd  or  even. 

4,  We  shall  not  enter  into  a  full  discussion  of  the  pro- 
perties of  Zonal  Harmonics  of  the  Second  Kind.  They  will  be 
found  very  completely  treated  by  Heine,  in  his  Handhuch  der 
Kugelfunctionen.  We  will  however,  as  an  example,  investi- 
gate the  expression  for  -r-^  in  terms  of  Qi^^,  Qi^.^...     , 

Recurring  to  the  equation 

+  {2i  +  l)oQ,{v)PM  +  ... 
we  see  that 

+  (2i  +  l)5^,(.)^^^  +  (2.'+3)Q,,W^^^ 
Now  we  have  seen  (Chap.  ii.  Art,  22)  that 
^^  =  (2^-^)  P,.,(/.)+  (2.-5)  P,_M  +  ... 

Hence     ^-^^  =  (2/+  1)  P^  +  (2i-  3)  P^ij.)  +... 
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^^^  =  (2/ +  9)  P,,M  +  (2^  +  5)  P,,M 

+  {2i+l)FM  +  ... 
And  therefore  the  coefficient  of  -PiO".)  in  the  expansion 

n     d  1  . 

01  -J IS 

a/jb  V  —  fi 

{2i+l)  {(2^+3)  Q,,,W  +  (2z  +  7)  ^,,3(.)  +  (2^  +  ll)  <?,.,,(.)+...}. 
Again, 


11  d      1  d      1         ^ 

And  , +  , =  0. 

(iv  V  —  jM      a/j,  V  —  fi 


Hence,  comparing  coefficients  of  P.  (/i), 


^t~  ==  -  ^^^'  +  ^^  ^-^  ^'^^  -  '^^'  +  ^^  ^'•-  ^' 


Hence  it  follows  that 
and  thei'efore  that 

/\. W  du  =  ^^  [QAv)-Q^M]' 

5.  By  similar  reasoning  to  that  by  which  the  existence  of 
Tcsseral  Harmonics  was  established,  we  may  prove  that  there 
is  a  system  of  functions,  which  may  be  called  Tesseral  Har- 
monics of  the  Second  Kind,  derived  from  T/"^)  in  the  same 
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manner  as  Q^  is  derived  from  P..    The  general  type  of  such 
expressions  will  be 


and  this  when  multiplied  by  cos  a(f)  or  sin  a^,  will  give  an 
expression  satisfying  the  differential  equation 

|(1  -  /.«)  I^y  U+  [i  {i  +  1)  (1  -  H'")  -  <r']  U=  0, 

and  which  may   be   called   the  Tesseral  Harmonic  of  the 
second  kind,  of  the  degree  i  and  order  a. 


CHAPTER  yi. 

ELLIPSOIDAL  AND   SPHEROIDAL   HARMONICS. 


1.  The  characteristic  property  of  Spherical  Harmonics 
is  thus  stated  by  Thomson  and  Tait  (p.  400,  Art.  587). 

"A  spherical  harmonic  distribution  of  density  on  a  spheri- 
cal surface  produces  a  similar  and  similarly  placed  spherical 
harmonic  distribution  of  potential  over  every  concentric 
spherical  surface  through  space,  external  and  internal." 

The  object  of  the  present  chapter  is  to  establish  the  ex- 
istence of  certain  functions  which  possess  an  analogous  pro- 
perty for  an  ellipsoid.  They  have  been  treated  of  by  Lamd, 
in  his  Legons  sur  les  fonctions  inverses  des  transcendantes  et 
les  fonctions  isothermes,  and  were  virtually  introduced  by 
Green,  in  his  memoir  On  the  Determination  of  the  Exterior 
and  Interior  Attractions  of  Ellipsoids  of  Variable  Densities, 
(Transactions  of  the  Cambridge  Philosophical  Society,  1835). 
We  shall  consider  them  both  as  functions  of  the  elliptic  co- 
ordinates (as  Lam^  has  done)  and  also  as  functions  of  the 
ordinary  rectangular  co-ordinates ;  and  after  investigating 
some  of  their  more  important  general  properties,  shall  pro- 
ceed to  a  more  detailed  discussion  of  the  forms  which  they 
assume,  wheii  the  ellipsoid  is  a  surface  of  revolution. 

2.  For  this  purpose,  it  will  be  necessary  to  transform 
the  equation 

dW     d'V    d'V     ^  ,^,    ^ 

into  its  equivalent,  when  the  elliptic  co-ordinates  e,  v,  v  are 
taken  as  independent  variables.  If  a,  b,  c  be  the  semiaxes 
of  the  ellipsoid,  the  two  sets  of  independent  variables  are 
connected  by  the  relations 
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d'-^v   '  ¥  +  v   '  c'  +  v 

Thus  a'  +6,  &"  +  6,  c*  +  6  are  the  squares  on  the  semiaxes 
of  the  confocal  ellipsoid  passing  through  the  point  x,  y,  z. 

o*  +  v,  ¥  +  V,  &  +  V,  the  squares  on  the  semiaxes  of  the 
confocal  hyperboloid  of  one  sheet. 

a*  +  v,  W  +  V ,  &  +  V ,  the  squares  on  the  semiaxes  of  the 
confocal  hyperboloid  of  two  sheets. 

Thus,  6  is  positive  if  the  point  a;,  y,  z  be  external  to  the 
given  ellipsoid,  negative  if  it  be  internal. 

And,  if  o^  be  the  greatest,  c*  the  least,  of  the  quantities 
a\  h\  c\ 

e  will  lie  between  —  c*  and  oo  , 

V  „         „  -^'    »    -c', 

V  „         „         -a     „    -6. 

dJ^V     d^V    d'V 
3.     Now   -j-^  +  -Ts  +  zn  =  0  is  the  condition  that 
ax        dy        cLz 

taken  throughout  a  certain  region  of  space,  should  be  a  mini- 
mum. In  the  memoir  by  Green,  above  referred  to,  this 
expression  is  transformed  into  its  equivalent  in  terms  of  a 
new  system  of  independent  variables,  and  the  methods  of  the 
Calculus  of  Variations  are  then  applied  to  make  the  resulting 
expression  a  minimum.  We  shall  adopt  a  direct  mode  of 
transformation,  as  follows : 

Suppose  a,  yS,  7  to  be  three  functions  of  x,  y,  z,  such  that 

V'^a  =  0,  v')S  =  0,  vV  =  0 (1), 

such  also  that  the  three  families  of  surfaces  represented  by 
tlie  equations  a  =  constant,  /3  =  constant,  7  =  constant,  inter- 
,  sect  each  other  everywhere  at  right  angles,  i.e.  such  that 
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<7y5  J7      d^  d'y     dl3  dy  _  ^      dy  dx     dy  dx      dydx  _ 
dx  dx      dy  dy     dz  dz        '     dx  dx     dy  dy      dz  dz       ' 

dxdjS     d^djS     did^     ^^ 

cl  _| L-  _j JZ  —  0 (2). 

dx  dx     dy  dy     dz  dz        ^  '' 


Then 


dV_dVda.     dVd§     dVdy 
dx       dx  dx      djS  dx      dy  dx ' 

dn^_  d'Vfd^     ^f^]'     d'Vfdf^ 


dx'      dr  \dxj  "^  d^'  [dxj  "^  df  \dx) 

^d^Vdl3 dy         dn^ dydx     ^  dW  da d_^ 
"^  d/:idy  dx  dx        dydx  dx  dx     ^  dxdji  dx  dx 

dVd^x     dVd'^      dVd'y 
da.  dx''      djS  dx''      dy  dx'' ' 

J'F         d'^V     . 
Y  .f  and  -7^  being  similarly  formed,  we  see  that,  when  the 

three  expressions  are  added  together,  the  terms  involving 

dV    dV     dV      .„    ,.  ,         ,  ,.   •  /nN  n     1       ' 

,-  ,  -777 ,  ~Y~  ^^'^^  disappear  by  the  conditions  (1),  and  those 

^^'^^^'^"-  my '  U".'  Jacl  ^'^  *^'  conditions  (2).     Hence 

_.     d'V(/dy.V      (dxV      uh\' 

^^^MU-'Uj-'idz) 

'^  d/3'\[dx)  '^[dyj  '^[dzj 


d_n'{{dy\\(cIyV     fdy\n^ 
"*■  dy'  \\dx)  "^  \dyj  "*"  [dz 


4.     Now,  let 


a  = 


d^jr 


[{a:  +  f)(b'  +  \lr)(/^^]^' 
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r  -c2 

J  V 


All  these  expressions  satisfy  the  conditions  (1),  for  a  is 
the  potential  of  a  homogeneous  ellipsoidal  shell,  of  proper 
density,  at  an  external  point,  and  yS  and  7  possess  the  same 
analytical  properties. 

Again,  a  is  independent  of  v  and  v,  and  is  therefore  con- 
stant when  6  is  constant.  Similarly  /9  is  constant  when  v  is 
constant,  and  7  is  constant  when  v  is  constant.  Hence  a,  /3, 
7  satisfy  the  conditions  (2). 


Now 


©"-(I)'-©" 


{sy+(iy+(ij 


{^a'  +  e)  [(}'  +  €)  (c'  +  €)[  \dxj       \dy  J       \dz 
222 
And  -,— -  +  , ./--  +  -T— -  =  1. 

(a^  +  ef  "^  (6*  +  ep  "*"  (c*  +  e)'^j  dx  ~  d'  +  e  ' 

with  similar  expressions  for    ,    and  -y- .      Hence,  squaring 
and  adding, 


^— -.--^-iier-dy-©]-- 


x'  y 

But  from  the  equations 

9  2  f2  9  9  * 

+  7;/-rT  +  :7  ,  _=-i'  t.^  ,  t +7:2-— -  +  ...,  ..-  ^' 


a'  -{■  e      b'  +  e     v  +  e       '  a^  +  v  '  b^  -^  v      c''  +  v 

*  jS  is  a  purely  imagiuary  quantity.    We  inay,  if  we  please,  wrile  .^  -  ip' 
for/3. 
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'  "T"    I.*      1     ..'    "'      Ji     ,     ..'   ~    ^f 


tt'+u       6*4-1;'       c'  +  v 
we  deduce 


1- 


x^  y^  z^      _   (ft)  —  e)  (<B  —  y)  («  —  v') 


o)  being  any  quantity  whatever.  For  this  expression  is  of 
0  dimensions  in  &>,  e,  f,  v,  it  vanishes  when  co  =e,  u,  or  u', 
and  for  those  values  of  w  only,  it  becomes  infinite  when 
<u  =  —  a^,  —  h\  or  —  c\  and  for  those  values  of  co  only,  and  it  is 
=  1  when  0)  =  00  . 

From    this,    multiplying  by   a^  +  w,   and   then    putting 
o)  =  —  a^  we  deduce 

a  result  which  will  be  useful  hereafter. 

Again,  differentiating  with  respect  to  co,  and  then  putting 
o)  =  e, 

of'  y^  ^       _         (e  —  v)  (e  —  v) 


(y.r j  "^  Uyj  ^  U^  y  ""^  (e  -v)  (e  -  u')       "  ' 

/dzV      fdotV      (dA^^ 4 

\dx)       \dy/        \dz)       (e  —  v)  {e  —  v)  ' 

(u- u)  (u  -  e)  (e-u)  (^  ^  aa        ^  a/3"' 

The  equation  V^F=  0  is  thus  transformed  into 
,,       ^d'V     ,        ,.d'V     .         .d'V     ^ 
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V 


+  (6  -  1.')  [[(v  +  a^)  (v  +  J^)  (u  +  c^)]4  J^J  V 

+  {v-e)  [{(.'  +  a^)  (v'+  h')  iv'  +  c')]i  ^,]  V. 

5.  A  class  of  integrals  of  this  equation,  presenting  a  close 
analogy  to  spherical  harmonic  functions,  may  be  investigated 
in  the  following  manner.  Suppose  ^  to  be  a  function  of  e, 
satisfying  the  equation 

\[{e  +  a'){e  +  h')(e  +  c')]^  ^ 


de 


V=  {me  +  r)  E, 


m  and  r  being  any  constants. 

Then,  if  H  and  H'  be  the  forms  which  this  function 
assumes  when  v  and  v  are  respectively  substituted  for  e, 
the  equation  V'^F=  0  will  be  satisfied  by  V=EHH'. 

6.  We  will  first  investigate  the  form  of  the  function 
denoted  by  E,  on  the  supposition  that  ^  is  a  rational  integral 
function  of  €  of  the  degree  n,  represented  by 

„_,      n(n  —  \)      __, 
e»  +  np/'  '  +  ^^^p,/  ^  +  ...  +2\- 

We  see  that 


+ 


—  n 


{n  -  1)  (e  +  a')  (e  +  h')  (e  +  c')  le""''  +  (n  -  2)  p/'- 


•+  P,j 


{n-2){n-S)  ] 

H :^; Y    +{)i~\)p^e 


1.2 


F.  H. 
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Hence  writing 


(6  +  a')  (e  +  ¥)  {e  +  c')  =e'  +  Sf/  +  Sf,e  4/3, 
we  see  that 

^  - 1)  (e«  +  Sf/  +  3/,6  +/3)  je"-  +{n-2)  ^^e""^ 


{n-2)(n-  8) 
1.2 


+^— ^^^:^-^^^v"-+...+p„. 


=  (me  +  r)  je"  +  np^e'^-'  +  ^^^^^  p.e""^  +  . . .  +  ^1 . 

Hence,  equating  coefficients  of  like  powers  of  e,  we  get 


=  wwpj  +  r, 


,3((^)(«_^)^,  +  2(«-l)/.,.+/j^ 


n  [  71  +  - )  =  m, 


.  - 1)  {{n  -  2)^,  +  3/J  + 1  {(,^  -  1)  ^,^  +  2/J 


1.2 


n  (w  —  1) 


(  8  1 

or,  as  they  may  be  more  simply  written, 


w  ( w  +  2  )  =  m, 
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(re  -  1)  (n  -  2  j  p^^  +  Sn/A  =  nnvp^  +  r, 


w  (ti  —  ]  ) 


It  thus  appears  that  p^  is  a  rational  function  of  r  of  the 
first  degree,  p^  of  the  second,  p^  of  the  ?i"',  and  when  tlie 
letters  2^1,  l^a-'i^fi  ^^^^  been  eliminated,  tlie  resulting  equa- 
tion for  the  determination  of  r  will  be  of  the  {n  +  Ij"*  degree. 
Each  of  the  letters  p,,  p^-.-p^  will  have  one  determinate 
value  corresponding  to  each  of  these  values  of  r;  and  we 

have  seen  that  m  =  nin  +  -^\.    There  will  therefore  be  (n  +  1 ) 

values  of  E,  each  of  Avhich  is  a  rational  integral  expression 
of  the  n^  degree,  n  being  any  positive  integer. 

7,  But  there  will  also  be  values  of  E,  of  the  tf"  degree, 
of  the  form 


(6+6^)H6+04e»-'  +  (n-l)(?/-^+^-^^^2./"+-+^„-, 


We  thus  obtain 

=  (e  +  a^i  (6  +  V) (e  +  c') {n-\)  |e"-^  +  {n  -  2)  q^e''-' 

(n  -  2)  (n  -  3)       „.,  ] 

8—2 
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d 


[{e-^a'){e-^h'){e  +  c')]^ 


E 


=  {(e  +  a?)  [e  +  ¥) (e  +  c'^)]l\\\^^±I}^  +  (e  +  a^)l (e  +  ^ 
+  (6  +  a^)i  (6  4-  6^)|  (n  -  1)  |e"-='  +  (ri  -  2)  ^^^e"-' 


(n-2)(7i-3)^  „-. 


1.2 


^,e" '  +  ...  +  ?„ 


+  (e+a^)'  (e+Z>^)  (e+c^)  (w-1)  (w-2)  L''-«+(n-8)  ^/ 


+ 


(,i_3)(^_4)^   „. 


1.2 


^/'°  +  ...  +  ^„_, 


Hence 


(,,  _  1)  |e-  +  (^  -  2)  ^/-  +  (^-^^)^(^--^)  ^^e"-  + . . .  +  ^„.  j 
+  (e  +  a')  (e  +  6')  (e  +  c'')  (n  - 1)  (w  -  2)  |€"-='  +  (w  -  3)  q/"-' 


+ 


(n_3)(n-4)_„ 


(?7ie  +  r)  |e""'  +  (7i  -  1)  q/ 


1.2 


^/'^  +  ...+^„. 


+ 


{n-r){n-^)       „. 


1.2 


^/   +--  +  ^n-,  ; 


(71  -  1)  ^1  +  w  -  2J  =  m, 


{n-l)\2a^  +  \{}f  +  e)  +  ^^{n~2)q^ 


+  (.„,_  1)  (7i-2)(a'-'  +  J^  +  c^  +  (n-3)gJ  =  (n-l)77i^,  +  r, 
(^  -  1)  |(^'  +  aV  +  aT)  <7„_,  +  (n  -  2)  aW  <7„_,j  =  r^„.,. 
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By  a  similar  process  to  that  applied  above,  we  shall  find 
that  r  is  determined  by  an  equation  of  the  n^^  degree,  and 

that  m=(w  — 1)  (  w— ^j,   and  that  each  of  the  letters  g-,, 

5'2"-?'n-i  is  ^  rational  function  of  r.  Thus,  there  will  be  n 
solutions  of  the  form 

(6  +  V)^  (e  +  &f  {e"-  +{n-l)  q^e"^  +  . . .  +  q^,). 

There  will  also  be  n  solutions  of  a  similar  form,  in  whicli 
the  factors  (e  +  c')^  (e  +  a^)^,  (e  +  a^)^  (e  +  b^)^  are  respectively 
involved.  Hence,  the  total  number  of  solutions  of  the  n^^ 
degree  will  be  4n  +  1. 

8.     We  may  now  investigate  the  number  of  solutions  of 

the  degree  n-+  -  ,  n  being  any  positive  integer.     These  will 

be  of  the  following  forms  :  three  obtained  by  multiplying  a 
rational  integral  function  of  e  of  the  degree  n  hy  (e  +  a^)s 
(e  +  b^)^,  (e  +  c")^,  respectively,  and  one  by  multiplying  a 
rational  integral  function  of  e  of  the  degree  n  —  1  by  the 
product 

{{e  +  a'){€  +  b'){e  +  c')]K 

An  exactly  similar  process  to  that  applied  above  will 
shew  us  that  there  will  be  ?n- 1  solutions  of  each  of  the 
first  three  kinds,  and  n  of  the  fourth.  Hence  the  total  number 
of  such   solutions  will  be  3  (ji  +  1)  +  n,  or  4/i  +  3,  that  is 


hi) 


+1. 


To  sum  up  these  results,  we  may  say  that  the  total 
number  of  solutions  of  the  n^^  degree  is  4>i  + 1,  n  denoting 
either  a  positive  integer,  or  a  fraction  with  an  odd  numerator, 
and  denominator  2. 

Similar  forms  being  obtained  for  U,  H' ,  we  may  proceed 
to  transform  the  expression  EHH'  into  a  function  of  x,  y,  z. 

9.     Consider  first  the  case  in  which 

E=  e"  +  np^e^-'  +  — '^i>,6"-  +  ...  +p„. 
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Write  this  under  the  form 

E={€-(o^  (e-ft)J...(6-(»J. 
Then  /T  =  (u  -  wj  (u  -  wj  . . .  (7^  -  «J, 

H'  =-{v' -  (O^  [V  -  W^  ...  {V    -  ft)J . 

Hence 
EHE'  =  (e  -  wj  {v  -  &) J  {v  -  ft) J ...  (e  -  « J  {v  -  «„)  (v  -  ft)„). 

Now  we  have  shewn  (see  Art.  4  of  the  present  Chapter) 
that  (e  —  ft) J  {y  —  &)J  {v  —  w^ 

=  (ft),  +  a')  (ft),  +  b')  (ft),  +  O  f-.r^^  +  jX-  +  -/—  -  1]  . 

^  '      ^^  '      ^^  '      'W+co^    6+ ft),    c^+ft),      y 

Each  of  the  factors  of  EHH'  being  similarly  transformed, 
we  see  that  EHH'  is  equal  to  the  continued  product  of  all 
expressions  of  the  form 


(ft)  +  a^)  (ft)  +  i^)  (ft)  +  c^)  f  ^  +  ^^— + -^ 


+  ft) 


the  several  values  of  w  being  the  roots  of  the  equation 

_  „_,     w  (w  —  1)       „  „ 

ft)"  +  np,a>^  '  +  ~IT2~" i'^"     +  •  •  •  +^"  =  ^• 

As  this  equation  has  been  already  shewn  to  have  {n  +  1) 
distinct  forms,  we  obtain  {n  + 1)  distinct  solutions  of  the 
equation  V^F=0,  each  solution  being  the  product  of  n 
expressions  of  the  form 

^^  if'  z^ 


O?  +  ft)       6''  +  ft)        c''  +  ft) 

That  is,  there  will  he  n  +  1  independent  solutions  of  the 
degree  2w  in  x,  y,  z,  each  involving  only  even  powers  of  the 
variables. 

10.  To  complete  the  investigation  of  the  number  of  solu- 
tions of  the  degree  2w,  let  us  next  consider  the  case  in  which  E 

=  (e+J^)i(e+c^)4|e"-V(ri-l)i),e-+^|:::^V.^»-+...+^„_,|. 
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The  object  here  will  be  to  transform  the  product 

since  the  other  factors  will,  as  already  shewn,  give  rise  to  the 
product  of  71  —  1  expressions  of  the  form 


o?  +  a>      W  +  a>      c*  +  w 

Now,  by  comparison  of  the  value  of  a?  given  in  Art.  4, 
we  see  that 

(e  +  If)  {v  +  If)  {v  +  ¥)  {e+c')  {v  +  c')  {v'  +  c') 

=  Qf  -  c')  {V  -  a')  {c'  -  a')  (c'  -  If)  fz\ 

Hence,  we  obtain  a  system  of  solutions  of  the  form  of 
the  product  of  (n  —  1)  expressions  of  the  form 

x^  f  z^         . 


a'+w     }f-\-(o     c^  +  o> 

multiplied  by  yz.  Of  these  there  will  be  w,  and  an  equal 
number  of  solutions  in  which  zx,  xy,  respectively,  take  the 
place  of  yz. 

Thus,  there  will  be  4?i  + 1  solutions  of  the  degree  2??  in 
the  variables  of  which  n  +  1  are  each  the  product  of  n 
expressions  of  the  form 

x^     ,     y' 


.     „^ L-Jr 1 


a"'*  +  ft)      If  +  (o      <f  +  ay 

n  are  each  the  product  of  (n— 1)  such  expressions,  multiplied 

by  2/2, 

n  ...  ...  ...  zx, 

n  ...  ...  ...  xy. 

11.  We  may  next  proceed  to  consider  the  solutions  of  the 
degree  2?i  + 1  in  the  variables  x,  y,  z. 

Consider  first  the  case  in  which 
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Here  the  product  (e  +  a^)^  (v  +  a^)^  (v  +  a^)'^  will,  as  just 
shewn,  give  rise  to  a  factor  x  in  the  product  EHH'. 

Hence  we  obtain  a  system  of  solutions  each  of  which  is 
the  product  of  n  expressions  of  the  form 

multiplied  by  x.  Of  these  there  will  be  n  + 1,  and  an  equal 
number  of  solutions  in  which  y,  z,  respectively  take  the 
place  of  the  factor  x. 

Lastly,  in  the  case  in  which 

E={e  +  a')^  (€  +  h'f  (e  +  c')^  je"-^  +  (w  -  1)  p^e""' 

{n  -  1)  (w  -  2)       „_„  ] 

+  - f72 ^^'     +  -  +^"-|' 

Ave  see  that  in  EHH'  the  product 

will  give  rise  to  a  factor  xyz. 

Hence  we  obtain  a  system  of  solutions  each  of  which  is 
the  product  of  [n  —  1)  expressions  of  the  form 

x^  it  ^ 

a^'  +  ft)     6'  +  G)     c'  +  ft) 

multiplied  by  xyz.     Of  these  there  will  be  n. 

Thus  there  will  be  472  +  3  solutions  of  the  degree  2n  +  1 
in  the  variables,  of  which 

{n  + 1)  are  each  the  product  of  /i  expressions  of  the  form 

;/  z" 


-^ V-rf^ V  -k 1  multiplied  by  x, 

{n  +  1)  are  each  the  product  of  n  such  expressions,  multiplied 
(n  +  1)  ...  ....  ...  ...  z, 
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n  are  each  the  product  of  {n  —  1)  such  expressions,  multi- 
plied by  xyz. 

12.    Now  an  expression  of  the  form  C .  EHIT,  C  being 
any  arbitrary  constant,  is  an  admissible  value  of  the  potential 

a^     V^     z^ 
at  any  point  within  the  shell  '-^  +  jj+  -§  =  1.      But  it   is 

not  admissible  for  the  space  without  the  shell,  since  it 
becomes  infinite  at  an  infinite  distance.  The  factor  which 
becomes  infinite  is  clearly  E,  and  we  have  therefore  to 
enquire  whether  any  form,  free  from  this  objection,  can  be 
found  for  this  factor.  We  shall  find  that  forms  exist,  bearing 
the  same  relation  to  JS  that  zonal  harmonics  of  the  second 
kind  bear  to  those  of  the  first. 


Now  considering  the  equation 


[{(^+«^)(^ 


d 


V'=(me-\-r)  U, 


(6  +  &^)(e  +  0}i^ 

which  we  suppose  to  be  satisfied  by  putting  U=  E,  we  see 
that,  since  it  is  of  the  second  order,  it  must  admit  of  another 

particular  integral.     To  find  this,  substitute  for   U,  E  ivde, 

we  then  have 


[ 


((6+a^)(6  +  i'^)(e  +  c^)|i^' 


U 


=  \^[{e+o?)  {e  +  F)  {e  +  c')]^  |1  E , ^vde 


{(6  +  a^)(6  +  i^)(e  +  c^)}4-f 


+  [[e  +  a'){e+ly){e^-c')]^Ev; 
U 


[[e-\-a'){e  +  y'){e  +  c')]^ 


ae 


I  d 
de 


■!■■ 


vde 


+  2  1(^+  ^')  C^  +  c'}  +  {e  +  c')  (e  +  a')  +  {e  +  a')  {e  +  h')]  Ev 
+  ie  +  a^{e+¥)(e  +  c')(^^.v  +  E^y 
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Now,  since  by  supposition,  the  equation  for  the  determi- 
nation of  Z7  is  satisfied  by  putting   U=E,  it  follows  that 

when  Elvde  is  substituted  for  U,  the  terms  involving  \vde 

will  cancel  each  other,  and  the  equation  for  the  determina- 
tion of  V  will  be  reduced  to 

Idv  .  2  dE     1  /^_         1  1    \  _  0 

whence    log  v  +  2  log  ^  +  log  |  (e  +  aj")  (e  +  h'^)  (e  +  c')}* 

=  log  v^+2  log  E^+  log  ahc, 
Vq  and  Eq  being  the  values  of  v  and  E,  corresponding  to  6  =  0. 
E.^  '  ahc 


V  de      E  ae 


Hence 


v  =  v^ 


•*•  ^/^^^  =  ^«^»^  "^^  •  ^ Ie^  Ke  +  a^)  (fn-  b^)  {e  +  c^)}^  ' 

We  may  therefore  take,  as  a  value  of  the  potential  at 
any  external  point, 

F=  v,E'  ahc  ETIH'  f -r . 

For  this  obviously  vanishes  when  e=  go  .  It  remains  so 
to  determine  v^  that  this  value  shall,  at  the  surface  of  the 
ellipsoid,  be  equal  to  the  value  C.EIIH',  already  assumed 
for  an  internal  point.     This  gives 

Hence,  putting  v^ .  E^ .  ahc  =  V^,  we  see  that  to  the  value 
of  the  potential 

V.ERH'I -^ i, 
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for  any  internal  point,  corresponds  the  value 


1 


■y    Tp'TTTJ 

for  any  external  point. 

13.  We  proceed  to  investigate  the  law  of  distribution  of 
density  of  attracting  matter  over  the  surface  of  the  ellipsoid, 
corresponding  to  such  a  distribution  of  potential. 

Now,  generally,  if  hn  be  the  thickness  of  a  shell,  p  its 
volume  density,  the  difference  between  the  normal  compo- 
nents of  the  attraction  of  the  shell  on  two  particles,  situated 
close  to  the  shell,  on  the  same  normal,  one  within  and  the 
other  without  will  be  Aiirphn.  This  is  the  attraction  of  the 
shell  on  the  outer  particle,  minus  the  attraction  on  the  inner 
particle. 

But  the  normal  component  of  the  attraction  on  the  outer 

particle  estimated  inwards  is  — j-  . 

And,  if  V  denote  the  potential  of  the  shell  on  an  in- 
ternal particle,  the  normal  component  of  the  attraction  on 

,    .       dV 
it  estimated  inwards  is  —    ,    . 

an 

XT  A       ^  d^        "^^ 

Hence  ^-n-pon  =  -r ,  - . 

'  an       dn 

dV_dVdx     dVdy     dV  ds 
dn      dx  dn      dy  dn      dz  dn' 

And  -Y-  is  the  cosine  of  the  inclination  of  the  normal  at 
dn 

the  point  oc,  y,  z  to  the  axis  of  x,  and  is  therefore  generally 

equal  to  e  ~,  — ,  e  denoting   the   perpendicular   from   the 
^  a  +€ 

centre  on  the  tangent  plane  to  the  surface 

o  9  Q 

a'+e^h'+e^c'  +  e 


whence 


or 
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And  we  have  shewn  that 

(a''-6^)(a^-c^        ' 

2dx_      1 
X  de      a^  +  e ' 

X     _^dx 
a'  +  e  ~    ^  ' 

•  ^  _  9    ^'^ 
'  '  dn  de' 

Similarly  ^=2^^,   $=2.1^ 

dn  de      dn  de^ 

'  '  dn  [dx  de  '^  dy  de'^  dz  de)  ~       de  ' 

Similarly  ^'  =  2e^'. 

dn  de 


Now     V'=V^.EHH 


/«QO 

Jo  ^ 


de 


de  ''  de  Jo  K'{{e  +  a')  (e  +  ^'0  (e  +  c')]'' ' 


And       V=V^.EHE''  ^^ 


W{{<i'-\-e){h'  +  e){&-re)f' 
therefore,  generally, 

-K.ERH' — i r. 

But,  when  the  attracted  particle  is  in   the   immediate 
neighbourhood  of  the  surface,  e  =  0.     Hence,  the  first  line 
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dV  ...  dV 

of  the  value  of  —  becomes  identical  with  the  value  of  — ,— , 
de  ae 

and  we  have 

dV^_dV^  yHH'   1 

de       de  ~~    "  E,   abc' 

Eq  denoting  the  value  which  E  assumes,  when  e  =  0. 

TTTJ'    1 
Hence,  ^Trphn  =  2e  F^  — ^ r- . 

-/->A     CLOG 

But  ^n,  being  the  thickness  of  the  shell,  is  proportional  to 
e,  and  we  may  therefore  write  f"  =  f-  »  ^ci  being  the  thick- 
ness of  the  shell  at  the  extremity  of  the  greatest  axis ; 

^  F,  g  1  mr 

"  ^     27r  ha  abc    E,    ' 

and  this  is  proportional  to  the  value  of  F  corresponding  to 
any  specified  value  of  e,  since  HH'  is  the  only  variable 
factor  in  either. 

Hence  functions  of  the  kind  which  we  are  now  considering 
possess  a  property  analogous  to  that  of  Spherical  Harmonics 
quoted  at  the  beginning  of  this  Chapter.  On  account  of 
this  property,  we  propose  to  call  them  Ellipsoidal  Harmonics, 
and  shall  distinguish  them,  when  necessary,  into  surface  and 
solid  harmonics,  in  the  same  manner  as  spherical  harmonics 
are  distinguished.  They  are  commonly  known  as  Lamp's 
Functions,  having  been  fully  discussed  by  him  in  his  Legons. 
The  equivalent  expressions  in  terms  of  x,  y,  z  have  been  con- 
sidered by  Green  in  his  Memoir  mentioned  at  the  beginning 
of  this  chapter,  and  for  this  reason  Professor  Cayley  in  his 
"  Memoir  on  Prepotentials,"  read  before  the  Royal  Society 
on  June  10,  1875,  calls  them  "  Greeniaus." 

We  may  observe  that  the  factor 

J_  a      1_ 

47r  ha  abc 

is  equal  to  -. — r-?^ ,  and  therefore  also  to  -, ?:r  or  ^r-  7  ^-  - 

^TTOcoa  4:ircaob       ^irabbc 
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Hence,  it  is  equal  to 


4<7r 

-^  {bc8a  +  caBb  +  ahSc) 
o 


or  to 


volume  of  shell ' 
and  the  potential  at  any  internal  point 


=  I  volume  of  shell  x  EE^  .pi    — 


and  the  potential  at  any  external  point 

r                        de 
=  i  volume  of  shell  x  EE^ .  p      ' • 

where  for  p  must  be  substituted  its  value  in  terms  of  v and  v. 

14.    We  will  next  prove  that  if  V^,   V^  be  two  different 
ellipsoidal  harmonics,  dS  an  element  of  the  surface  of  the 

ellipsoid,   lie  V^V^d8  =  0,  the  integration  being  extended  all 

over  the  surface. 

We  have  generally 

And  throughout  the  space  comprised   within  the  limits  of 
integration,  V  F,  =  0,  V  V^  =  0.     Hence 


IK^-'T-y^'I^-o- 


Now  it  has  been  shewn  already  that  V^,  V^  are  each  of 
the  form  EHH' ,  where  E  is  a  function  of  e  only,  H  the  same 
function  of  v,  H'  of  v.     We  may  therefore  write 

v.=fAe)A{y)A(-'), 

and  similarly  F,  =/^  (e)/^ (v) /^ {v). 
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Hence  F;  ^^  =  F,  F,  ^tt  ' 

"'^  de    '-  de    '^'n/;(^)  'fM' 

Now,  all  over  the  surface,  e  =  0.     Hence 

Hence,  unless  - 1    „,  —  4hW-  =  0,  which  cannot  happen 

unless  the  functions  denoted  by  j^  and  /,  are  identical*,  or 
only  differ  by  a  numerical  factor,  we  must  have 


//' 


eFF„^^=0. 


Now  e  is  proportional  to  the  thickness  of  the  shell  at 
any  point.     Calling  this  thickness  he,  we  have  therefore 


//■ 


heV,VA8=Q. 


Hence,  adding  together  the  results  obtained  by  integrating 
successively  over  a  continuous  series  of  such  surfaces,  we  get 


jjjv,V,dxd^dz  =  0; 


F, ,  Fg  now  representing  solid  ellipsoidal  harmonics,  and  the 
integration  extending  throughout  the  whole  space  comprised 
within  the  ellipsoid, 

*  This  may  be  shewn  more  rip;oroualy  by  integratinf;  throuf^h  the 
space  bouuded  by  two  coiit'ocal  ellipsoida,  dcliued  by  the  values  \  aud/x  of  e. 
We  then  get,  as  iu  the  text, 

Now  the  factor  within  { }  cannot  vanish  for  all  values  of  X  and  fi,  unless  the 
functions  devoted  by  /^  and  /^  be  identical,  or  only  differ  by  a  uumericaJ. 
factor. 


128  ELLIPSOIDAL  AND  SPHEROIDAL   HARMONICS. 

15.     It  will  be  well  to  transform  the  expression 

(eVVdS 


to  its  equivalent,  in  terms  of  v,  v. 

For  this  purpose  we  observe  that  if  ds,  ds  be  eletaents  of 
the  two  lines  of  curvature  through  any  point  of  the  ellipsoid, 
dS=  dsds. 

Now, 
ds^  is  the  value  of  dx^  +  dy"^  +  dz^  when  e  and  v  are  constant, 
ds^  ...  ...  ...  6  and  v 

therefore  if  e  and  v  do  not  vary, 
'^dx        dv 

.'.  dx  =  -— dv. 

2  a'  +  u 

It/  1      2? 

Similarly  dij  =  -  y^-^—  dv,    dz=.-:  ~ dv ; 


Again,  differentiating  with  respect  to  co  the  expression 

X^  11^  z^ 

obtained  for  —. \-  ~- h  — 1,  we  sret 


(a'  +  (of      {b'  +  cof      (c^  +  o)f      {d'  +  co)  {b'  +  CO)  (c^  +  ft)) 

,  (v  -  <»)  (g  -  f>) (e  -w)  (i^  -  ft))  (u'  -ft)) 

"^  (a'  +  co)  (JJ'  +  co)  {c'  +  ft))  "^  (f/+  ft))^  ip-^c^  (c'  +  ft))  "^  •  ■  • ' 
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therefore,  putting  w  =  v, 

A  similar  expression  holding  for  ds'^  we  get 
J  (j2 ^  _  J^ W -vY  {e-v)  {e-v) 2j,/. 

1  cr'  v"  ^ 

Again,  -5  =  r:2-T— V.  +  //,a  ,  „n2  + 


^        (g-  ^')  (g  -  ^) 

writing  e  for  w  in  the  expression  above  ; 

•  •  ^ '''^  ~     16  {d'+v) {ly+v){(f+^/)(a;+v){l^+v}  {c'+v)  "^^ '*''  * 

It  has  been  shewn  that,  integrating  all  over  the  surface, 
the  limits  of  v  are  —  c*  and  —J',  those  of  v,  —  b*  and  —  a'. 

Hence,    F^,   V^,   denoting  two  different  ellipsoidal  har- 
monics 


J  -b^  J  -a=  [{a'+v) {b'+v)  {c'+v)  {a'+v)  {b'+v)  {c'+v')]^ 


=  0. 


The  value  of  the  expression   1 1 1  Wdxdydz,  or  its  equiva- 


lent 


a6c 


f-c'    r-6» 


-«'  /-«-' V^{v-v)dvdv' 

{{a'+v)  (b'+v)  {c'+v)  {a'+v)ib'  +  v')  (c^+u')]'- 


in   any  particular   case,  is   most  conveniently  obtained  by 
expressing  V  as  a  function  of  cc,  y,  z, 

F.  H.  9 
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16.  Before  proceeding  further  with  the  discussion  of  ellip- 
soidal harmonics  in  general,  we  will  consider  the  special  case 
in  which  the  ellipsoid  is  one  of  revolution.  We  must  enquire 
what  modification  this  will  introduce  in  the  quantities  which 
we  have  denoted  by  a,  /3,  7,  viz. 

„^r d± 

_  /""'    d^ 

and  in  the  differential  equation 

We  will  first  suppose  the  axis  of  revolution  to  be  the 
greatest  axis  of  the  ellipsoid,  which  is  equivalent  to  supposing 
b^  =  c^.  To  transform  a  and  7,  put  a'  +  -\^  =  6\  a^+  €  =  rj', 
a^  -\-  v  =  (o^ ;  we  then  obtain 

«=2|     ...     ''?     ...^-J^.W-'+t"'-^')' 


^f"         dd  1         ,      (a'-b')-^-(o 

Jo  &'-a'+  ¥      {a^ -h'f     °  {a'  -b'}'^+(o 

To  transform  /3,  we  must  proceed  as  follows. 

Put    yjr  =  -  c^  cos'^  'ST  —  b^  sin*  ot,    v  =  —  c"  cos*  (p  -  6''sin*  cf), 
wc  then  get  generally 

b'  +  f={b'-  c')  cosV,    c^  +  ^  =  (c*  -  b')  sinV  ; 
dyjr  =  2  (c"  —  6*)  COS  -57  sin  ct  dt^  ; 
.    o  ^  __2      [^      d^ 2  V  -  1  0 
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Hence,  ^  =  -1  (,=  _„.+ J-)  | . 

Also,  €  =  if  —  a^  v  =  ft)"''  —  a',  v=—  h\  and  our  differential 
equation  becomes 

+  iv'  -a'  + 1')  1(0,^  -  a'  +  b')  £|'  r 

This  equation  may  be  satisfied  in  the  following  ways. 

First,  in  a  manner  altogether  independent  of  0,  by  sup- 
posing V  to  be  the  product  of  a  function  of  rj  and  the  same 
function  of  <«,  this  function,  which  we  will  for  the  present 
denote  by/(?;)  or/(w),  being  determined  by  the  equation 

or  -f-  \(ai'-  a'  +  h')  4-]  fH  =  mf(oi}. 

ciu)  [ '  riw  j  ■ 

0-   i! 
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Secondly,  by  supposing  -  -^   a  constant  multiple  of    V, 
=  —  cr^V,  suppose. 

Our  equation  may  then  be  written 

-  [rf  -  aV  If)  [{of-  a'  +  If)  ^^X  V 

-  a'  {a'  -  V)  [{of  -a'  +  If)  -  {rf  -  a' +  }f)]  V=  0, 

which  may  be  satisfied  by  supposing  the  factor  of  V  inde- 
pendent of  (j)  to  be  of  the  form  F  (97)  ^{0)),  where 

|(^^_  a^  +  J^)  ~-Xf(v)  -  a'  {a'  -  h')  Firj)  =  m  {v'-a'+b')  F{v), 

\{co'-  a'  +  h')  ^'f{<^)  -0-'  {a'-h')  F{w)  =m {of-o?+¥)  F{co). 

The  factor  involving  ^  will  be  of  the  form 

A  cos  (T<f>  +  B  sin  acf). 
Now,  returning  to  the  equation 

A 

we  see  that,  supposing  the  index  of  the  highest  power  of  tj 
involved  \rLf{rf)  to  be  i,  we  must  have  m  =  i{i+l). 

Now,  it  will  be  observed  that  rj  may  have  any  value 
however  great,  but  that  of,  which  is  equal  to  o^  +  v,  must 
lie  between  a^  —  }f  and  0.  Hence,  putting  tu"  =  (a^  —  Z>^)  yu,^, 
where  ixf'  must  lie  between  0  and  1,  we  get 

I,  {(1  - /^^)  |,}/i(«^  -  O^ /^l +  « (^ +  !)/{(«''- Z'^)  VI  =  0. 


I 


[{v'-o^'-^^')2^f{v)  =  mf{ri), 


ELLIPSOIDAL  AND   SPHEROIDAL  HARMONICS.  133 

Hence  this  equation  is  satisfied  by  f[{a^—  h"^)^  fi]  =  CPi, 
C  being  a  constant ;  and  supposing  C  =  1  we  obtain  the 
following  series  of  values  for  /  (tu), 

i  =  0,    /(a,)  =  l, 
.       t  =  l,    /(a))  = 

i=2,    f[a>)  = 
i  =  3,    /(a,)  = 


2  id'  -  b')     ' 


Exactly  similar  expressions  may  be  obtained  iov  f(ri),  and 
these,  when  the  attraction  of  ellipsoids  is  considered,  will 
apply  to  all  points  within  the  ellipsoid.  But  they  will  bt; 
inadmissible  for  external  points,  since  t)  is  susceptible  of  in- 
definite increase. 

The  form  of  integral  to  be  adopted  in  this  case  will  bo 
obtained  by  taking  the  other  solution  of  the  ditferential 
equation  for  the  determination  of  /(17),  i.e.  the  zonal  har- 

f      "n 
monic  of  the  second  kind,  which  is  of  the  form  Qi 


[a'-b-")^)' 
where 

dd 


H{a^-b^)^~^\a^\^U 


P. 


^       \ 


{d'-a'+b') 


[{a'-b')^\ 
Or,  putting  rj'^  =  (a"  -  ¥)  v^,  &"  =  {a^  —  h"^)  \^,  we  may  write 

17.  We  may  now  consider  what  is  the  meaning  of  the 
quantities  denoted  by  rj  and  eo.  They  are  the  values  of  "^ 
which  satisfy  the  equation 
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and  are  therefore  the  semi-axes  of  revolution  of  the  surfaces 
coiifocal  with  the  given  ellipsoid,  which  pass  through  the 
point  X,  y,  z.  One  of  these  surfaces  is  an  ellipsoid,  and 
its  semi-axis  is  t).  The  other  is  an  hyperboloid  of  two  sheets 
whose  semi-axis  is  <o. 

Now,  if  6  be  the  eccentric  angle  of  the  point  x,  y,  z, 
measured  from  the  axis  of  revolution,  we  shall  have 

x^  =  Iff  cos'^  6. 

But  also,  since  rf^,  w^,  are  the  two  values  of  ^^  which 
satisfy  the  equation  of  the  surface, 

Hence  oi^  =  (a?  —  If)  cos^  6, 

and  we  have  already  put 

Avhence  the  quantity  which  we  have  already  denoted  by  /i 
is  found  to  be  the  cosine  of  the  eccentric  angle  of  the  point 
X,  y,  z  considered  with  reference  to  the  ellipsoid  confocal 
with  the  given  one,  passing  through  the  point  x,  y,  z.  We 
have  thus  a  method  of  completely  representing  the  potential 
of  an  ellipsoid  of  revolution  for  any  distribution  of  density 
symmetrical  about  its  axis,  by  means  of  the  axis  of  revo- 
lution of  the  confocal  ellipsoid  passing  through  the  point 
at  which  the  potential  is  required,  and  the  eccentric  angle 
of  the  point  with  reference  to  the  confocal  ellipsoid.  For 
any  such  distribution  can  be  expressed,  precisely  as  in  the 
case  of  a  sphere,  by  a  series  of  zonal  harmonic  functions  of 
the  eccentric  angle. 

18.  When  the  distribution  is  not  symmetrical,  we  must 
have  recourse  to  the  form  of  solution  which  involves  the  factor 
A  cos  a(f)+  B  sin  acf).  It  will  be  seen  that,  supposing  F  to 
represent  a  function  of  the  degree  i,  and  putting  m  =  «  (t-1-1), 
the  equation  which  determines  F{(o)  is  of  exactly  the  same 
form  as  that  for  a  tesseral  spherical  harmonic.  For  F{r}),  if 
the  point  be  within  the  ellipsoid,  we  adopt  the  same  form, 
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if  without,  representing  the  tesseral  spherical  harmonic  by 

Ty)  I — '^ — .1 ,  or  r>^  (z/),  we  adopt  the  form 
[{a  -  hy) 

d\ 


^/^K^)£^ 


19.  It  maybe  interesting  to  trace  the  connexion  of  sphe- 
rical harmonics  with  the  functions  just  considered.  This  may 
be  effected  by  putting  h^  =  a'.  We  see  then  that  tj  will  become 
equal  to  the  radius  of  the  concentric  sphere  passing  through 
the  point,  and  rf  —  a?  4-  }f  will  become  equal  to  rj^  Hence 
the  equation  for  the  determination  of/ (77)  will  become 

|(.'|)/W -•(.•+!)/(.)- 

which  is  satisfied  by  putting  /(t;)  =  V,  or  rf^^'^^K  The  former 
solution  is  adapted  to  the  case  of  an  internal,  the  latter  to 
that  of  an  external  point. 

With  regard  to /(to),  it  will  be  seen  that  the  confocal 
hyperboloid  becomes  a  cone,  and  therefore  on  becomes  inde- 
finitely small.     But  u,,  which  is  equal  to  — , ,  remains 

finite,  being  in  fact  equal  to  -  or  cos  9.  Hence /(/i)  becomes 
the  zonal  spherical  harmonic. 

Again,  the  tesseral  equations,  for  the  determination   of 

F{t]),  F {(o),  become 

which  are  satisfied  by  F{r])  =rf  or  77"''"^". 

And,  writing  for  w^  [a^-l^)  fi,  w^e  have,  putting  F{o))=x{P'), 

|(/x^  -  1)  |,[  X  if^)  +  <^'X  (/-)  =  **  i^  + 1)  if-"  - 1)  %  if^l 
which  gives  x  (z^)  =  TJ-'^Kl-)- 


1:36 


ELLIPSOIDAL  AND  SPHEROIDAL  HARMONICS. 


20.  We  will  next  consider  the  case  in  which  the  axis  of 
revolution  is  the  least  axis  of  the  ellipsoid,  which  is  equi- 
valent to  supposing  c^  =  6^  To  transform  a  and  yS,  put 
c"  +  yjr  =  O'^,  c^  +  €  =  Tf,  c^  +  v  =  ft)*,  we  thus  obtain 

J^a'-e'+ff'     {a'-c',i  , 


]^dr  —  c  - 


-J  tan" 


To  transform  7,  we  must  proceed  as  follows : 

Put  i/r  =  —  a^  sin*  ot  —  J*  cos*  ot,    v  =  —  a*  sin*  <^  —  J*  cos*  <^, 
we  then  get,  generally, 

a'  +  a|r  =  (a*  -  ¥)  cos*  -ir,     W^<f'  =  -{a^-  W)  sin*  ct, 

c*+'\/r=c*— a*sin*<^  — J*cos*(^,  c?'\/r  =  — 2  (a*— 6*)  sinw  cosizr  c?ct. 

Hence 


7  =  2 


2<^ 


<^  (a*  sinV  + 1^*  cos'^'ot  -  c'")^         (a*  -  c'') 
Hence, 


—  ifa*=6* 


also, 


2  2 

6  =  77  -c, 


U  =  ft)'  —  c  , 
v  =  —  a*, 
and  our  differential  equation  becomes 


(a^  _  c*  +  ft)*)  |(a*  -  c*  +  7?*) -^1 V 

_  (a^_c^  +  ^'')  |(a«_c*  +  a,*)  -Xv 
[  a(o) 

+  (,?*-ft,*)(a*-c*)jJ=0. 
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We  will  first  consider  how  this  equation  may  be  satisfied 
by  values  of  V  independent  of  <f). 

We  may  then  suppose  Fto  be  the  product  of  a  function 
of  77,  and  the  same  function  of  w,  this  function,  which  we  will 
suppose  to  be  of  the  degree  i,  being  determined  by  the 
equation 

On  comparing  this  with  the  ordinary  differential  equa- 
tion for  a  zonal  harmonic,  it  will  be  seen  that,  on  account 
of  a**  being  greater  than  c*,  the  signs  of  the  several  terms  in 
the  series  for  /(t;)  will  be  all  the  same,  instead  of  being 
alternately  positive  and  negative.     We  shall  thus  have 


(a^-Oi' 


.  ^       _  Soy*  +  30  (g'  -  c')  v'  +  3  (g'  -  c' 


and  generally 

We  will  denote  the  general  value  of /(t^)  by  p^  \ ,  , 

or,  writing  77=  {a^  —  &)-v,  by  j^A^)' 
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For  external  points,  we  must  adopt  for  /  (7;)  a  function 
which  we  will  represent  by  qA  ^  ^  I,  or  q.{v),  which  will 
be  equal  to 


Hia'^-cf 


(^0^+a'-c') 


It  is  clear  that  /(w)  may  be  expressed  in  exactly  the 
same  way.  But  it  will  be  remembered  that  r]^  and  a>^  are 
the  two  values  of  ^'^  which  satisfy  the  equation 


a'-c'  +  y  '  ^* 

Hence  77,  as  before,  is  the  semi-axis  of  revolution  of  the 
confocal  ellipsoid  passing  through  the  point  {x,  y,  z).  But 
»7^w^  —  —  {a?  —  c^)  z^,  an  essentially  negative  quantity,  since 
a'"*  is  greater  than  c^  Hence  w^  is  essentially  negative.  Now, 
if  6  be  the  eccentric  angle  of  the  point  {x,  y,  z)  measured 
from  the  axis  of  revolution,  we  have  z^  =  rj'  cos^^.     Hence 

and  therefore  w^  =  —  {a?  —  c^)  cos'^  6 

=  —  (a^  —  &)  yLt^  suppose. 

Hence  the  equation  for  the  determination  of  /(&))  assumes 
the  form 

the  ordinary  equation  for  a  zonal  spherical  harmonic.  Hence 
we  may  write 

fjb  being  the  cosine  of  the  eccentric  angle  of  the  point  x,  y,  z, 
considered  with  reference  to  the  confocal  ellipsoid  passing 
through  it. 
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21.  We  have  thus  discussed  the  form  of  the  potential, 
corresponding  to  a  distribution  of  attracting  matter,  sym- 
metrical about  the  axis.  When  the  distribution  is  not 
symmetrical,  but  involves  <^  in  the  form  A  cos  a<f>-\-  B  sin  cr^, 
we  replace,  as  before,  P^  (a*-)  by  r/'>  (/a),  and  pjfx.)  by  a 
function  tj''^^{i^)  determined  by  the  equation 


and  q,  [v]  by  */-)  {v)  \ 


d\ 


.</<^)(x)r(\^+i) 

22.     As  an  application  of  these  formulae,  consider  the  fol- 
lowing question. 

Attracting  matter  is  distributed  over  the  shell  whose 

x'^     y^  +  z^ 
surface  is  represented  by  the   equation  — ^  +     j^i —  —  1>   so 

that  its  volume  density  at  any  point  is  P,  (fi),  /x  being  the 
cosine  of  the  eccentric  angle,  measured  from  the  axis  of 
revolution ;  required  to  determine  the  potential  at  any 
point,  external  or  internal. 

The  potential  at  any  internal  point  will  be  of  the  form 

CPAh)P.{v) (1), 

and  at  an  external  point,  of  the  form 

C'PMQA^') (2), 

where  {a^  —  V^)^  v  =  ih.e  semi-axis  of  the  figure  of  the  con- 
focal  ellipsoid  of  revolution  passing  through  the  point  (/a,  v). 

Now  the  expressions  (1)  and  (2)  must  be  equal  at  the 

surface  of  the  ellipsoid,  where  v  =  — — , . 

(a  -  by 

Hence 
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But  generally 


\ 


J  V 


Hence 

6X 


(«•■:- &2)^ 


We  may  therefore,  putting  V  =  AP.  \ yl ,  write 

'  [{(I'-by) 

c=aqJ ^-X, 

and  we  thus  express  the  potentials  as  follows  : 

^Pi  {h)  Pi  (^)  Qi  Itv-7^1  f  at  an  internal  point, 

■^Pi  {h)  Qi  (^)  Pi  \  —^ ^r  at  an  external  point. 


Or,  substituting  for  Q.  its  value  in  terms  of  P., 

d\ 


F.  =  AP,  if,)  P,  (v)  P,  I — ~     \  f        = 
at  an  internal  point, 

V, = AP,  if,)  p,  {v)  p,  I — '^i  r.=^^ — - 

at  an  external  point. 

Now,  to  determine  A,  we  have,  Ba  being  the  thickness 
of  the  shell  at  the  extremity  of  the  axis  of  revolution, 
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-  JL     A_  (^  _  ^\ 
4:7r  8a  .  7)  \  dr]         drj  Jn  =  a 

~  ^7r6ad'-b^\dv       dp  ]"=,  ,"    i 


p  .        a 


"  47r  8a  a^  -  6"^  ^*  ^^^  ^^^l?  " 

Hence,  if  p  =  P^  (/x),  we  obtain 

A  =4!7r =  4!TrbBb. 

a 

And  we  thus  obtain 


(/..-) 


T\  =  47r  &S6  P,  ifi)  P,  {v)  P,  [ ^    ,1  I  = 


rfx. 


MhhP,{H)P,{v)Q,^^^-^^, 


«     )  T. , .  r       d\ 


V,  =  ^irUh  P,  (/*)  P,  I  — ^ X  P  (,.)  ["  - 


{x):o^-\) 


a 


i7rhSbPMQ,(v)P       _^ 

[{a  —  0  )2) 

If  the  shell  be  represented  by  the  equation 
it  may  be  shewn  in  a  similar  manner  that  we  shall  have 
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C 


■V,  =  4'7radaPMp,{v)q, 


V^=4<7radaP,{fi)q,{v)p 


23.  We  may  apply  this  result  to  the  discussion  of  the 
following  problem. 

If  the  potential  of  a  shell  in  the  form  of  an  ellipsoid  of 
revolution  about  the  greatest  point  he  inversely  proportional 
to  the  distance  from  one  focus,  find  the  potential  at  any 
internal  point,  and  the  density. 

If  the  potential  at  P  be  inversely  proportional  to  the 
distance  from  one  focus  S,  and  H  be  the  other  focus,  we  have, 

HP+SP  =  2rj,    HP-SP=2co, 
.'.    SP  —  rj  —  ro. 

Hence  if  M  be  the  mass  of  the  shell,  V^  the  potential  at 
any  external  point, 

'^        7]—  CO 

M         1 


[a  —  b  )  ^ 

Now,  by  what  has  just  been  seen,  the  internal  potential, 
correspundiijg  to  P^  (fj.)  (J^  {v),  is 


PMPM-   ^^'^~^'' 


Hence,  if  V^he  the  potential  at  any  internal  point, 
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a 


Q. 


y.-^,ti.i.i)-^^^PMFA^)- 


And  the  volume  density  corresponding  to  P,  (ji)  Q^  {v)  is 

Hence  the  density  corresponding  to  the  present  distri- 
bution is 


P 


[[a'-b')^ 
If  Fj  had  varied  inversely  as  HP,  we  should  have  had 

'      V  +  fo' 

and  our  results  would  have  been  obtained  from  the  foregoing 
by  changing  the  sign  of  to,  and  therefore  of  /u,. 

24.  Now,  by  adding  these  results  together,  we  obtain 
the  distributions  of  density,  and  internal  potential,  corre- 
sponding to 

Tj  —  (Ji         rj-\-  Oi  7}    —w 

or,  in  geometrical  language, 

y  _M^      M  ^      SP  +  HP 
^~  SP^'liP  JSP.  Up  ' 

=  M  multiplied  by  the  axis  of  revolution  of  the  confocal 
ellipsoid,  and  divided  by  the  square  on  tho  conjugate  semi- 
diameter.  We  may  express  this  by  saying  that  the  potential 
at  any  point  on  the  ellipsoid  is  inversely  proportional  to  the 
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square  on  the  conjugate  semi-diameter,  or  directly  as  the 
square  on  the  perpendicular  on  the  tangent  plane. 

Corresponding  to  this,  we  shall  have,   writing  2k  for  i, 
since  only  even  values  of  i  will  be  retained, 


'     {a'-b')i 


X  (4A;  +  1) -^ — — -^  P,,  (fi)  P^  (v), 


a       I 


{a'  -  b') 


^      ^TT  {a' -  b')HBb     ^  p  « 

k  being  0,  or  any  positive  integer. 

Again,  subtracting  these  results  we  get 

'a  —  1  2  2  > 

^        7}  —  CO         r]  +  (O  7)  —CO 

=  M  multiplied  by  the  distance  from  the  equatoreal  plane, 
and  divided  by  the  square  on  the  conjugate  semi- diameter. 

This  gives,  writing  2^  +  1  for  i, 

a    I ^ 


47r(a^-6'0i^»S6 


P„. 


2*+l)/„2  72,1 


[{d'-b')^} 

25.  In  attempting  to  discuss  the  problem  analogous  to 
this  for  an  ellipsoid  of  revolution  about  its  least  axis,  we  see 
that  since  its  foci  are  imaginary,  the  first  problem  would  re- 
present no  real  dist^-ibution.  But  if  we  suppose  the  external 
potential  to  be  the  sum  or  difference  of  two  expressions,  each 
inversely  proportional  to  the  distance  from  one  focus,  we 
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obtain  a  real  distribution  of  potential — in  the  first  case 
inversely  proportional  to  the  square  on  tbe  conjugate 
semi-diameter,  in  the  latter  varying  as  the  quotient  of  the 
distance  from  the  equatoreal  plane  by  the  square  on  the 
conjugate  semi-diameter. 

It  will  be  found,  by  a  process  exactly  similar  to  that  just 
adopted,  that  the  distributions  of  internal  potential,  and 
density,  respectively  corresponding  to  these  will  be : 

In  the  first  case 

Tc  being  0,  or  any  positive  integer. 
In  the  second  case 


Pik+l 


\{a'-c')^) 


^     47r(a^-c7^aSa     ^  ^ 

7c  being  0,  or  any  positive  integer. 

26.  AVe  may  now  resume  the  consideration  of  the  ellip- 
soid with  three  unequal  axes,  and  may  shew  how,  when  the 
potential  at  every  point  of  the  surface  of  an  ellipsoidal  shell 
is  known,  the  functions  which  we  are  considering  may  be 
employed  to  determine  its  value  at  any  internal  or  external 
point.     We  will  begin  by  considering  some   special   cases, 

r.  H.  10 
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by  which  the  general  principles  of  the  method  may  be  made 
more  intelligible. 

27.  First,  suppose  that  the  potential  at  every  point  of 
the  surface  of  the  ellipsoid  is  proportional  to  a?  =  -^  suppose. 

In  this  case,  since  x  when  substituted  for  V,  satisfies  the 
equation  v'^=  ^>  ^^  ^^^  *^^^  ^o|  ^'"^^  ^1^^  be  the  potential 
at  any  internal  point.  But  this  value  will  not  be  admissible 
at  external  points,  since  x  becomes  mfinite  at  an  mhnite 
distance. 

Now,  transforming  to  elliptic  co-ordinates 

And  the  expression 

p  Ja/t 

"Jo  if  +  a')  Kt  +  a')  (t  +  ^")  (t+c')}^ 

"satisfies,  as  has  already  been  seen,  the  equation  ^'V^O,  is 

equal  to  V,~  at  the  surface  of  the  ellipsoid,  and  vanishes 

at  an  infinite  distance.     This  is  therefore  the  value  of  the 
potential  at  any  external  point.    It  may  of  course  be  written 

lo^  r ^± , 

a  je  {yir  +  a')  [(t  +  a')  {^r  +  V)  if  +  c')Y 


28.     Next,  suppose  that  the  potential  at  every  point  of 
the  surface  is  proportional  to  ?/s=Fo|^-,  suppose.    In  this 
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case,  as  in  the  last,  we  see  that,  since  yz  when  substituted 
for  V,  satisfies  the  equation  v^F=  0,  the  potential  at  any 

internal  point  will  be  F^  ?- ;  while,  substituting  for  y,  z  their 

values  in  terms  of  elliptic  co-ordinates  we  obtain  for  the 
potential  at  any  external  point 

y,y^  r d± 


J  a 


29.    We  will  next  consider  the  case  in  which  the  po- 
tential,  at  every  point  of  the  surface,  varies  as  x^  =  V^  -, 

Cv 

suppose.  This  case  materially  differs  from  the  two  just  con- 
sidered, for  since  x^  does  not,  when  substituted  for  V,  satisfy 
the  equation  y"  F=  0,  the  potential  at  internal  points  cannot 
in  general  be  proportional  to  x\  We  have  therefore  first  to 
investigate  a  function  of  x,  y,  z,  or  of  e,  v,  v  Avhich  shall 
satisfy  the  equation  \7"F=0,  shall  not  become  infinite  within 
the  surface  of  the  ellipsoid,  and  shall  be  equal  to  x"^  on  its 
surface. 

Now  we  know  that,  generally 
(b'i-(o)(c''  +  (o)x'+{c'+(o)ia'  +  a>)7/+(a''+(o){b'  +  co)z^ 

-  (a'  +  co)  {b^  +  (o)  (c'  +  0))  =  (e  -  o))  (u  -  co)  (v  -  (o). 

And,  i£  Oj^,  0^  be  the  two  values  of  w  Avhich  satisfy  the 
equation 

{b'+a>)  (c'  +  co)  +  (c^+o))  (a'+oo)  +  (a^+w)  (6^+a))  =  0...(l), 

we  see  that 

V^(6-^J(u-^J(u'-^,)=0, 
■and  ^^e-e,){v-d;){v'-d,)  =  0. 

And,  by  properly  determining  the  coefficients  A^,  A^,  A^, 
it  is  possible  to  make 

A-^A,{e-e,) (y-exv-e,)+A,{e-ej{v-0:){v'-e^)...{2) , 

=  — f-  when  6'cV  +  c V/  +  a'b'z'  -  a'bV  =  0. 

10—2 
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Hence,  the  expression  (2)  when  A^,  A^,  A^  are  properly 
determined  will  satisfy  all  the  necessary  conditions  for  an 
internal  potential,  and  will  therefore  be  the  potential  for 
every  internal  point. 

Now,  we  have  in  general 

iij" + e^  (c^  +  e^  x' + (c^  +  e;)  {a? + e,)  if  +  (a'  +  e,Y  {h' + e,)  ^ 

-  {a^ + e^)  (6^  +  e,)  (c^ + e,)  =  (e  -  e,)  (u  -  e^)  (v'  -  e,) 

(b' + e,)  (c^  +  ^,)  x-^  +  (c^ + e^)  {a'  +  e,)  f  +  (a^ + e,)  {w  +  ^,)  z^ 

and,  over  the  surface 

hVx''+c'ay-ha'h'z'-a'hV  =  0. 

Hence,  ^  being  any  quantity  whatever,  we  have,  all  over 
the  surface, 

and  therefore,  putting  ^  =  —  a^ 

Hence,  the  right-hand  member  of  this  equation  possesses 
all  the  necessary  properties  of  an  internal  potential.  It 
satisfies  the  general  differential  equation  of  the  second  order, 
does  not  become  infinite  luithin  the  shell,  and  is  proportional 
to  x^  all  over  the  surface. 

"VVe  observe,  by  equation  (1),  that 
(Z/^  +  w)  (c^  +  «)  +  {^-\-w)  [en  (o)  +  {a' + co)  {h'+cc)  =  8  {d-co)  {6 -(d) 
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identically,  and  therefore,  writing  —  c^  for  «, 

Hence,  over  the  surface  of  the  shell, 


«•  = 


and  we  therefore  have,  for  the  internal  potential. 

This  is  not  admissible  for  external  points,  as  it  becomes 
infinite  at  an  infinite  distance.  We  must  therefore  substi- 
tute for  the  factor  e  —  ^^ 

df 


J  0 


(^-^J-^[(^  +  a^)(^+6-^)(^  +  c-0]" 


with  a  similar   substitution  for  e  —  0^,  thus  giving-,  for  the 
external  potential, 

i€-e,)(^^e,)(v'-0J)_  r rf^f- 


F.=^- 


•'  0 


dylr 


+  (e-^,)(u-^,)(^-^,)  r d^ 

^2  iK-^i)  (^'+^2)  ie  if-d^y  {(>^+«')  (^+z>"0  (t+Ol^ 

^  r d^ 

^.  ["_ d± ^  r d^jr __- 
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The  distribution  of  density  over  the  surface,  correspond- 
ing to  this  distribution  of  potential,  may  be  investigated  by 
means  of  the  formula 

^     2'7rda[de        de~/\^o' 
or  its  equivalent  in  Art.  13  of  this  Chapter.   We  thus  find  that 

_    1    a     V,    r        {v-0^){v'-e,) 
^      27r  da  Sabc  [    ^i'  (6,  -  0^)  {a'  +  6^) 


J  0 


dyjr 


{v-e;){v'-0;)    ^  r d± 

+i.r '-±- -i 

Jo  {{ir  +  a'){f  +  b'){-f  +  c')\U 

80.  The  investigation  just  given,  of  the  potential  at  an 
external  point  of  a  distribution  of  matter  giving  rise  to  a 
potential  proportional  to  x^  all  over  the  surface,  has  an  in- 
teresting practical  application.  For  the  Earth  may  be  re- 
garded as  an  ellipsoid  of  equilibrium  (not  necessarily  with 
two  of  its  axes  equal)  under  the  action  of  the  mutual  gravi- 
tation of  its  parts  and  of  the  centrifugal  force.  If,  then, 
V  denote  the  potential  of  the  Earth  at  any  point  on  or  with- 
out its  surface,  and  ft  the  angular  velocity  of  the  Earth's 
rotation,  we  have,  as  the  equation  of  its  surface,  regarded  as 
a  surface  of  equal  pressure, 

.*.  V  +  -^  H^  (a?  +  2/^)  =  a  constant,  IT  suppose. 

Hence,  if  a,  h,  c  denote  the  semi-axes  of  the  Earth,  we 
have,  for  the  determination  of  V,  the  following  conditions : 
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y=0  at  an  infinite  distance (2), 

2 


F=  n  -  J  O'  (^'  +  /)  when 


-!  +  C  +  ^=l (3). 

a      0      G 

The  term  IT  will,  as  we  know,  give  rise  to  an  external 
potential  represented  by 

n  r        ^^        -  r____it_____ 

1  1 

The  two  terms —^  nV,  — -fl^/',  will  give  rise  to  terras 

which  may  be  deduced  from  the  value  of  V^  just  given  by 
successively  writing  for  F^,  — -^flV,  and  —^Qj^W,  and  (in 
the  latter  case)  putting  6"  for  a^  throughout.     We  thus  get 

je-e;)  (v-e,)  (v'-e;)  r d^ 

^r d± a' /  a'       b^\ 

(e-e,)(v-9;)(v'-e,)  r df 

_^  r d^p^ 
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31.  Any  rational  integral  function  V  of  x,  y,  z,  which 
satisfies  the  equation  y^F=0,  can  be  expressed  in  a  series 
of  Ellipsoidal  Harmonics  of  the  degrees  0,  1,  2...i  in  x,  y,  z. 
For  if  V  be  of  the  degree  i,  the  number  of  terms  in  V  will 
Ije    (^•+l)(^'+2)(^  +  3)^     ^^^  ^^^  condition  v'F=0  is 

equivalent  to  the  condition,  that  a  certain  function  of  x,  y,  z 

of  the  degree  i  —  2,  vanishes  identically,  and  this  gives  rise 

(i  —  1)  i  ii  +  1) 
to  -^^ ^^7^^ conditions.     Hence  the  number  of  inde- 

0 

pendent  constants  in  V  is 

(^+  1)  (^  +  2)  {i+  3)      (^  -  1)  t  {i  +  1) 
6  6 

or  (i  +  l)^    And  the  number  of  ellipsoidal  harmonics  of  the 

1         3     i 

degrees  0,  1,  2..d  in  cc,  y,  z  or  of  the  degrees  0,  -,1,  ^...-^ 

in  e,  V,  v,  is,  as  shewn  in  Arts.  6  to  10  of  this  Chapter, 

1  +  3  +  5+.. .  +  2i  +  l, 

or  {i  + 1)\  Hence  all  the  necessary  conditions  can  be  satis- 
fied. 

32.  Again,  suppose  that  attracting  matter  is  distributed 
over  the  surface  of  an  ellipsoidal  shell  according  to  a  law  of 
density  expressed  by  any  rational  integral  function  of  the 
co-ordinates.  Let  the  dimensions  of  the  highest  term  in  this 
expression  be  i,  then  by  multiplying  every  term,  except  those 
of  the  dimensions  i  and  ^  —  1  by  a  suitable  power  of 

^     f     ^ 

we  shall  express  the  density  by  the  sum  of  two  rational  inte- 
gral functions  of  x,  y,  z  of  the  degrees  i,  i—  1,  respectively. 
The  number  of  terms  in  these  will  be 

(i±I)Ji±^+i(i+l)  or  (,•+!)'. 
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And  any  ellipsoidal  surface  harmonic  of  the  degree  i,  i  —  2... 
in  X,  y,  Zy  may,  by  suitably  introducing  the  factor 

a      0      c 

be  expressed  as  a  homogeneous  function  of  a;,  y,  z  of  the 
degree  i  ;  also  any  such  harmonics  of  the  degree  «— 1,  i— 3... 
in  X,  y,  z  may  be  similarly  expressed  as  a  homogeneous 
function  of  x,  y,  z  of  the  degree  i  —  1.  And  the  total  number 
of  these  expressions  will,  as  just  shewn,  be  {i+  l)^  hence  by 
assigning  to  them  suitable  coefficients,  any  distribution  of 
density  according  to  a  rational  integral  function  of  x,  y,  z 
may  be  expressed  by  a  series  of  surface  ellipsoidal  harmonics, 
and  the  potential  at  any  internal  or  external  point  by  the 
corresponding  series  of  solid  ellipsoidal  harmonics. 

83.  Since  any  function  of  the  co-ordinates  of  a  point  on 
the  surface  of  a  sphere  may  be  expressed  by  means  of  a  series 
of  surface  spherical  harmonics,  we  may  anticipate  that  any 
function  of  the  elliptic  co-ordinates  v,  v  may  be  expressed  by 
a  series  of  surface  ellipsoidal  harmonics.  No  general  proof, 
however,  appears  yet  to  have  been  given  of  this  proposition. 
But,  assuming  such  a  development  to  be  possible  at  all,  it 
may  be  shewn,  by  the  aid  of  the  proposition  proved  in 
Art.  15  of  this  Chapter,  that  it  is  possible  in  only  one  way, 
in  exactly  the  same  way  as  the  corresponding  proposition 
for  a  spherical  surface  is  proved  in  Chap.  IV.  Art,  11. 

The  development  may  then  be  effected  as  follows.  De- 
noting the  several  surface  harmonics  of  the  degree  i  in  x,  y,  z, 

or  \  in  u,  v\  by  the  symbols  V,<-^\   VP,...\T^^\  and  by 

F{v,  v)  the  expression  to  be  developed,  assume 

+  (7^(1)7/1)4....  + (7W7W  4.... 

Then  multiplying  by  eF/"^^  and  integrating  all  over  the 
surface,  we  have 

jeF {v,  v)  17-> dS  =  C^^le (F/-))' dS. 
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The  values  of  jeFiv,  v')F/<^)  dS,  and  of  Je  {V(''^y dS must 

be  ascertained  by  introducing  the  rectangular  co-ordinates 
X,  y,  z,  or  in  any  other  way  which  may  be  suitable  for  the 
particular  case.  The  coefficients  denoted  by  C  are  thus 
determined,  and  the  development  effected. 


EXAMPLES. 


o  -I  a  g 

1.  Prove  that   {&m  6)' = -^  P,-^F,  + :^P^. 

Why  cannot  (sin  Of  be  expanded  in  a  finite  series  of  spherical 
harmonics  1 

1  1  1  l  +  sins 

2.  Prove  that  1  +  2  ^i+  3  ^3  +  4  ^3+  -  =log —  . 


sm^- 


3.     Establish  the  equations 


oiP„  =  {2n-l),.P„_,-in-l)P„_,. 

4.  If  /u.  =  cos  6,  prove  that 

^,(^)  =  l-^(^  +  l)sin--4-.■■4■(-l^^-^^--=— (sin-^)-f... 

and  also  that 

J\{h-)  =  (- 1)'  +  (-  1)'''  i{i  + 1)  cos''  I  +  ... 

+  (-l)'-^'».pT:r7--(cos^^)  +... 

^      '       (\m)  \i  —  m  \        2/ 

5.  Prove    that,  if    a    be  greater    than    c,   and  i   any   odd 
integer  greater  than  m, 

6.  Prove  that  f  (^p\  dfj.  =  i  (i  +  1). 


156 


EXAMPLES. 


] 


7.     Prove  that,  when  w  =  ±  1,  -r-^  =  ^^=  \^/     . 


8.     Prove  that 


IP  P 

F    F  F 

P    P  P 


is  a  numerical  multiple  of 


9.     Prove  the  following  equation,  giving   any   Laplace's   co- 
efficient in  terms  of  the  preceding  one : 


P+j  =  ^P„  +  w  j  FJp  +  C, 


where  Cp  =  it-ix!  +  J^  - /x^  ^1  -  /x'^  cos  (w  -  w')  and  C  is  zero  if  n  be 
even,  and 

«+i  I  w  +  1 


(-1) 


2""Mi^(^+i)r 


,  if  w  be  odd. 


10,     If  i,  j,  k  be  three  positive  integers  whose  sum  is  even, 
pi'ove  that 


i 


\.^...{j+h-i-\)   1.3...(/c  +  *-i-l)   l.^  ...{l+j-h-l) 


2A...{j  +  k-i) 


2.4:...{k  +  i-j)  2A...{i+j-k) 

2A...{{+j  +  k)  1 


l.'S  ...  {i  + j +  k-l)     i+j+k+l' 
Hence   deduce   the  expansion  of  P^  in    a   series  of  zonal 
harmonics. 

11.  Express  a?y -v 'tf  ■\- yz -vy  +  z  as  a  sum  of  spherical 
harmonics, 

12.  Find  all  the  independent  symmeti-ical  complete  harmonics 
of  the  thii'd  degree  and  of  the  fifth  negative  degree. 

13.  Matter  is  distributed  in  an  indefinitely  thin  stratum  over 
the  surface  of  a  sphere  whose  radius  is  unity,  in  such  a  manner 
that  the  quantity  of  matter  laid  on  an  element  i^S)  of  the  surface 
is  l>S{\^ax+ly-{-cz-\-fx^-\-gy'^\hz'')t 
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where  a;,  y,  z  are  rectangular  co-ordinates  of  the  element  ^S  re- 
ferred to  the  centre  as  origin,  and  a,  h,  c,  f,  g,  h  are  constants. 
Find  the  value  of  the  potential  at  any  point,  whether  internal 
or  external. 

14.  If  the  radius  of  a  sphere  be  r,  and  its  law  of  density  be 
p  =  ax  +  hy  +  C2,  where  the  origin  is  at  the  centre,  prove  that  its 

potential  at  an  external  point  (^,  -q,  ^)  is  (a^  ■\-  hr]  +  c^)  where 

ii  is  the  distance  of  (^,  rj,  ^)  from  the  origin. 

15.  Let  a  spherical  portion  of  an  infinite  quiescent  liquid  be 
separated  from  the  liquid  round  it  by  an  infinitely  thin  flexible 
membrane,  and  let  this  membrane  be  suddenly  set  in  motion, 
every  part  of  it  in  the  direction  of  the  radius  and  with  velocity 
equal  to  S^,  a  harmonic  function  of  position  on  the  surface.  Find 
the  velocity  produced  at  any  external  or  internal  point  of  the 
liquid.  State  the  corresponding  proposition  in  the  theory  of 
Attraction. 

16.  Two  circidar  rings  of  fine  wire,  whose  masses  are  Jf  and 
M,  and  radii  a  and  a',  are  placed  with  their  centres  at  distances 
6,  h',  from  the  origin.  The  lines  joining  the  origin  with  the 
centres  are  perpendicular  to  the  planes  of  the  rings,  and  are  in- 
clined to  one  another  at  an  angle  6.  Shew  that  the  potential  of 
the  one  rin"  on  the  other  is 


MM'^"^-Q~,BF^Q^, 


where  B  ^^n.^^C^  5>^v^^(^- 1)(^- 2)  (^- 3)  ^^  ,_ 
2.2  2.2.4.4 

and  B\  and  Q^  are  the  same  functions  of  h'  and  a'  and  of  cos  6 
and  sin  6  respectively,  and  c  is  the  greater  of  the  two  quantities 
\/a'  +  6'^and  Ja!^TV\ 

17.  A  uniform  circular  wire,  of  radius  a,  charged  with 
electricity  of  line-density  e,  surrounds  an  uninsulated  concentric 
spherical  conductor  of  radius  c;  prove  that  the  electrical  density 
at  any  point  of  the  surface  of  the  conductor  is 

2cV  •2.r^-+'^-2.4a^^^^"^-2.4.6a«^  +-;' 

the  pole  of  the  plane  of  the  wire  being  the  pole  of  the  harmonics. 
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18.  Of  two  splierical  conductors,  one  entirely  surrounds  the 
otter.  The  inner  has  a  given  potential,  the  outer  is  at  the 
potential  zero.  The  distance  between  their  centres  being  so 
small  that  its  square  may  be  neglected,  shew  how  to  find  the 
potential  at  any  point  between  the  spheres, 

19.  If  the  equation  of  the  bounding  surface  of  a  homo- 
geneous spheroid  of  ellipticity  e  be  of  the  form 


r  =  <.(l-|cP.), 


prove  that  the  potential  at  any  external  point  will  be 

M     C-A 
r  r^ 


P.^ 


where  C  and  A  are  the  equatoreal  and  polar  moments  of  inertia 
pf  the  body. 

Hence  prove  that  Y  will  have  the  same  value  if  the  spheroid 
be  heterogeneous,  the  surfaces  of  equal  density  differing  from 
spheres  by  a  harmonic  of  the  second  ord«r. 

20.  The  equation  E  =  a{\  -^-ay)  is  that  of  the  bounding 
surface  of  a  homogeneous  body,  density  xmity,  differing  slightly 
in  form  and  magnitude  from  a  sphere  of  radius  a;  a  is  a 
small  quantity  the  powers  of  which  above  the  second  may  be 
neglected;  and  y  is  a  function  of  two  co-ordinate  angles,  such 
that 

where  F^,  Y^..,  Z^,  Z^  ...  are  Laplace's  'functions.  Prove  that 
the  potential  of  the  body's  attraction  on  an  external  particle, 
the  distance  of  which  from  the  origin  of  co-ordinates  is  r,  is 
given  by  the  equation 

^     47rV      47raaV  a  «"         xr  "I 
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21.  If  J/  be  the  mass  of  a  uniform  hemispherical  shell  of 
radius  c,  prove  that  its  potential,  at  any  point  distant  r  from 
the  centre,  will  be 

M     \M  (\        _    1     j,t^ 
2(;'^2?V2     '^    2.4  ^c^ 


.pr^_^44-.pX.. 


2.4.6    'c*     2,4.6.8    '  c« 


•)' 


""'  2r 


2r     2      V2^r*     2.4^r* 


c^_      3.5        ^c 


2.4.6'r''     2.4.6.8 


A^-..> 


according  as  r  is  less  or  greater  than  c;  the  vertex  of  the  hemi- 
sphere being  at  the  point  at  which  fi  —  l, 

22.  A  solid  is  bounded  by  the  plane  of  xy,  and  extends  to 
infinity  in  all  directions  on  the  positive  side  of  that  plane. 
Every  point  within  the  circle  v?  -k-if  =  «*,  2;  =  0  is  maintained  at 
the  uniform  temperature  unity,  and  every  point  of  the  plane  xy 
without  this  circle  at  the  uniform  temperature  0.  Prove  that, 
when  the  temperature  of  the  solid  has  become  permanent,  its 
value  at  a  point  distant  r  from  the  origin,  and  the  line  joining 
which  to  the  origin  is  inclined  at  an  angle  Q  to  the  axis  of  z  will 
be 

r      1      r'     1 . 3       r" 

_  1.3...(2^-1)  r- 

^      '       2.4...2t     ^-°'+>tt^'+i  ^••• 
if  r  <  a,  and 

2    '  r'     2.4     V*  ^      '       2.4...2t       "^  2'+i  ,.-^*  •  "  • 

if  r  >  a. 

23.  Prove  that  the  potential  of  a  cii'cular  ring  of  radius  c, 
whose  density  at  any  point  is  cos  m\^,  c\^  being  the  distance  of  the 
point  measured  along  the  ring  from  some  fixed  point,  is 
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2n-  cos  md>  (sm  6)   \  - — -. — -^ — ^ r-    -i:nri 

^  ^        ''(2.4.  6. ..2m  £//a"*  r"'*' 


+  - 


^"i'... 


2.4.6,..(2m  +  2)      c^/x"     r'""^*  '  '" 

1.3.5...(2A;-1)    (^'"P^,^^    c"-^^'*  ) 

"^  2.4.6. ..2(m  +  /fc      dfjT      r"+«+i"'" -j 

where  r  is  gi-eater  than  c.     If  r  be  less  than-  c,  r  and  c  must  h 
interchanged. 

24.  A  solid  is  bounded  by  two  confocal  ellipsoidal  surfaces,  aiK 
its  density  at  any  point  F  varies  as  the  square  on  the  perpendicula 
from,  the  centre  on  the  tangent  plane  to  the  confocal  ellipsoic 
passing  through  P.  Prove  that  the  resultant  attraction  of  sue] 
a  solid  on  any  point  external  to  it  or  forming  a  part  of  its  mas 
is  in  the  direction  of  the  normal  to  the  confocal  ellipsoid  passing 
through  that  point,  and  that  the  solid  exercises  no  attraction  on  ; 
point  "vvithin  its  inner  surface. 


S'/ 
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